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Abstract 

We present a review of parton saturation/Color Glass Condensate physics in the context of 
deuteron-gold (d+Au) colhsions at RHIC. Color Glass Condensate physics is a universal description 
of all high energy hadronic and nuclear interactions. It comprises classical (McLerran-Venugopalan 
model and Glauber-Mueller rescatterings) and quantum evolution (JIMWLK and BK equations) 
effects both in small-x hadronic and nuclear wave functions and in the high energy scattering 
processes. Proton-nucleus (or d + A) collisions present a unique opportunity to study Color Glass 
Condensate predictions, since many relevant observables in proton- nucleus collisions are reasonably 
well- understood theoretically in the Color Glass Condensate approach. In this article we review the 
basics of saturation/Color Glass Condensate physics and reproduce derivations of many important 
observables in proton (deuteron)-nucleus collisions. We compare the predictions of Color Glass 
physics to the data generated by d + Au experiments at RHIC and observe an agreement between 
the data and the theory, indicating that Color Glass Condensate has probably been discovered at 
RHIC. We point out further experimental measurements which need to be carried out to test the 
discovery. 
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1 Introduction 



Saturation/Color Glass Condensate physics is a rapidly developing field of strong interactions at 

high energy. Color Glass Condensate (CGC) physics describes high parton densities inside the hadronic 
and nuclear wave functions at small values of Bjorken x variable. It demonstrates how the gluon fields 
in the hadronic and nuclear wave functions reach their maximum allowable values in quantum chromo- 
dynamics (QCD) corresponding to ~ 1/g The CGC formahsm is also successfully applied 

to calculation of total, elastic and diffractive cross sections of high energy hadronic and nuclear scatter- 
ing. There, by resumming the strong gluon field dynamics, it resolves such long-standing questions as 
unitarity of the scattering S'-matrix |25, and infrared (IR) safety [221; which are known to be violated 
by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation P7| corresponding to the weak field 
limit of CGC. Saturation/CGC approach allows one to calculate particle production in hadronic and 
nuclear scattering [2H]-|11]- The resulting inclusive particle production cross sections are infrared-safe, 
which is a significant theoretical improvement over the IR-divergent perturbative QCD results j451 
making the small coupling approach to particle production self-consistent. One of the interesting appli- 
cation of particle production in CGC framework is understanding the initial conditions for the evolution 
of the quark-gluon system produced in heavy ion collisions toward the possible thermalization leading 
to formation of quark-gluon plasma (QGP) jUl US 1121 EH • 

Perturbative QCD (pQCD) has been extremely successful in describing the particle production 
data over a large kinematic window Applications of pQCD to particle production in high energy 
hadronic or nuclear collisions are based on the use of coUinear factorization theorems. The essence 
of a collinearly factorized cross section is the idea of incoherence. In other words, a hadronic cross 
section can be written as a convolution of parton distributions and fragmentation functions, which are 
universal non-perturbative objects that are also subject to perturbative evolution (DGLAP), with a hard 
scattering cross section involving partons, which is perturbatively calculable but is process dependent. 
Parton distribution functions are typically measured in Deep Inelastic Scattering experiments such as 
the ones performed at HERA, where it has been observed that the gluon and sea quark distributions 
grow very fast with decreasing Bjorken x. This fast growth can be understood in pQCD as driven 
by radiation of gluons with small Bjorken x via DGLAP evolution equations. Collinear factorization 
theorems are not exact and are violated by effects that are typically suppressed by inverse of the hard 
momentum transfer but can be enhanced by energy (or In or A dependent factors, which may be 
large at high energy and/or for large nuclei. This necessitates construction of a new formalism that 
does not rely on collinear factorization and can include these potentially large effects. The hint for this 
new formalism comes from pQCD itself, noticing that the rise of parton distribution functions can not 
continue for ever since it would lead to growth of hadronic cross sections at a rate which would violate 
unitarity. A weak coupling mechanism which can tame this fast growth is gluon recombination and 
saturation. Color Glass Condensate formalism is the natural generalization of pQCD in order to make 
it applicable to dense partonic systems. 

The extensive theoretical progress in the field of saturation/Color Glass has been summarized in sev- 
eral review articles, mostly concentrating on the issues of non-linear small-x evolution ^\ H^l ISO] • Our 
article here deals with saturation/ Color Glass Condensate physics putting more emphasis on particle 
production in proton (deuteron)-nucleus collisions {p{d)A) and in deep inelastic scattering (DIS). Many 
of the relevant particle production observables in pA and DIS have been well-understood theoretically 
in the saturation/Color Glass approach, at least at the partonic level [2Hl-|in]- It is therefore, very im- 
portant to be able to verify our theoretical understanding by comparing the predictions of CGC physics 
for particle production in p{d)A to the experimental data produced by deuteron-gold {d + Au) scat- 
tering program at Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) 
jST]-[nnj- Below we will review both the CGC predictions |61j-|65j and experimental data reported by 
RHIC experiments jS3j-[nnj. We will point out the apparent agreement between the two indicating a 
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possible discovery of Color Glass Condensate at RHIC jSI]. We will also review future experimental 
test which can be carried out to test the CGC discovery both hy d + Au program at RHIC and by pA 
scattering program at Large Hadron Collider (LHC) at CERN. 

The paper is structured as follows. We begin in Sect. |21with a general review of saturation/CGC 
physics. This review is by no means all-inclusive: we will concentrate on the material that we will 
need later in our discussion of particle production. We refer the interested reader who wants to learn 
more about various aspects of small-x evolution to the dedicated reviews in |1H1 EHl EOl IHE] • In Sect. 
12. II we discuss the BFKL evolution equation j27j along with its problems, such as violation of unitarity 
PH] and diffusion into infrared [20! • We also review the Gribov-Levin-Ryskin-Mueller-Qiu (GLR-MQ) 
evolution equation |T1I2]- We proceed in Sect. I2.2l bv discussing quasi-classical approximation in small-x 
physics. We review Glauber-Mueller multiple rescatterings in DIS |^ IHH] and McLerran-Venugopalan 
model of small-x wave functions of large nuclei |Sl IH El EH] • Quasi-classical regime at small-x is valid 
when X is small enough, so that coherent interactions of nucleons in the nucleus with the projectile are 
possible. This translates into the requirement for parton coherence length l^oh to be larger than the 
nuclear radius R 170] , 

Icoh = > (1) 

2 rriN XBj 

with rriN the nucleon mass and xbj the Bjorken x variable of a parton. Defining the rapidity variable 
Y = \nl/ Xbj we recast the condition as y > In A, with A the atomic number of the nucleus. 
On the other hand, when xbj becomes too small, BFKL evolution effects become important, breaking 
down the quasi-classical approximation. BFKL evolution brings in powers of \n{l/xBj) ~ ^ • 
Requiring for such effects to be small, OgY < 1, we obtain an upper bound on the allowable rapidity 
range. The applicability window for the quasi-classical approximation is then 

InA < Y < — . (2) 

as 

We also discuss in Sect. 12.21 how the saturation scale Qs arises in the quasi-classical limit. 

We continue our review of saturation/ Color Glass physics by re-deriving the Jalilian-Marian-Iancu- 
McLerran-Weigert-Leonidov-Kovner (JIMWLK) |71-|in|, [12]- [El and Bahtsky-Kovchegov (BK) [23 
[I^ non-linear evolution equations in Sect. 12.31 Quantum small- a; evolution corrections become impor- 
tant when asY > 1 [27|, such that 

Y>-. (3) 

as 

Eq. ((HI) gives a lower bound on the region of applicability of JIMWLK and BK evolution equations. 

We conclude the review of CGC by solving the non-linear evolution equations in Sect. 12.41 There we 
discuss the solution of linear (BFKL) evolution equation outside of the saturation region, demonstrate 
an interesting property of the solution of JIMWLK and BK known as geometric scaling inside the 
saturation region [TH [721 ^^^1 reproduce the derivation of extended geometric scaling outside of that 
region 173]. We demonstrate how saturation scale Qs grows with energy once the quantum evolution 
effects are included [ZH ESI IZOl- We explain how JIMWLK and BK evolution equations resolve the 
problems of the BFKL evolution by unitarizing the corresponding total cross sections of DIS and by 
prohibiting diffusion into the infrared, making the small-coupling approximation self-consistent f77'-f81j. 

We continue in Sect. [HI by deriving expressions for a number of particle production observables 
in pA collisions in the saturation/CGC framework. We start in Sect. 13.11 bv calculating inclusive 
gluon production cross section in pA in the quasi-classical approximation [S2l IH21 [HHl (HHI • We show 
that quasi-classical multiple rescatterings lead to Cronin enhancement [Hll of gluon production in pA 
[nH [n2l [HS|-[Hni; as can be seen from Fig. [30 We then proceed in Sect. 13.21 bv including the effects 
of quantum BK evolution in the expression for gluon production cross section in pA |3ll [^01 [Si- As 
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one can see from Fig. |^ the effect of small-a; evolution is to flatten the Cronin maximum introducing 
suppression of particle production at all transverse momenta (see Fig. I¥T|) jni]-|nil IHD ES- In Sect. 
13.31 we calculate valence quark production cross section in pA both in the quasi-classical limit and 
including small-x evolution |SZI-|ini- We move on to electromagnetic probes in Sect. 13. 4[ where we 
calculate prompt photon and dilepton production cross sections in pA collisions |93j- |lUUj . Finally, in 
Sect ]3. 51 we analyze two-particle correlations |101j - [TTU] . We rederive production cross section in pA for 
two gluons at mid-rapidity |in3j . for a quark and a gluon at forward rapidity |in3j and for qq pair at 
mid-rapidity [105]- [110]. 

We review some of the data generated by d + Au scattering program at RHIC |SI]-[nn| in Sect. El We 
show that the data reported by BRAHMS [HZI, PHENIX [HI], PHOBOS |S2| and STAR jMI experiments 
show Cronin-like enhancement of particle production at mid-rapidity. We discuss how this result, 
combined with suppression of produced hadrons in Au + Au collisions |lllj - |114j . serves as a control 
experiment for the signals of quark-gluon plasma formation in heavy ion collisions at RHIC |115j - 
|12Uj . We then review BRAHMS data at forward rapidity [23 IS3 ; indicating suppression predicted by 
saturation/CGC approach IHHIHSESl- This data is confirmed by preliminary results from PHENIX 
PHOBOS 58j and STAR jHOI, demonstrating that Color Glass Condensate has probably been 
discovered at RHIC |l2(Jj . We conclude Sect. ID by listing future experimental tests |12H 11221 H^^! 
which need to be carried out to verify the discovery of CGC at RHIC |12Uj . 

We conclude in Sect. Elby discussing some of the issues which we inevitably had to leave out in this 
review, including recent progress in our understanding of pomeron loop corrections to small-x evolution 
|124j - |125j and some exclusive processes |129j - P^^ . which can not be measured at RHIC. 

Throughout the paper we will use the following notation for the 4-vectors: for a 4-vector = 
(fo, fi, f2, ws) we define the light cone components by = {vq ± V3)/^/2 and combine the transverse 
components into a two-dimensional vector v = (fi,f2). The metric tensor is chosen in such a way that 

2 Overview of Saturation/Color Glass Condensate Physics 

In this Section we will review the developments and advances of high energy QCD which led to our 
modern understanding of the saturation/Color Glass Condensate physics. 

2.1 Early Developments 
2.1.1 The BFKL Equation 

A milestone in the development of small- a; physics was the derivation by Balitsky, Fadin, Kuraev, and 
Lipatov of what has become known as the BFKL equation j2Z|. The BFKL equation describes the 
behavior of scattering amplitudes and gluon distribution functions at asymptotically high energies. It 
does it by resumming leading logarithms of energy, i.e., powers of the parameter a^lns, where in the 
perturbative QCD regime the coupling constant is small, <^ 1, and, in Regge kinematics, the center 
of mass energy s of the scattering process is much larger than any other momentum scale involved. 

Derivation of the BFKL equation is rather complicated. For derivations in the transverse momentum 
space we refer the interested reader to P7|ll33j . as well as to Sect. l2.3.1l below. Derivation in transverse 
coordinate space can be found in [121 (201 1211 122]; as well as in Sect. 12.3.21 below. Here we are going to 
present a simple physical picture of the BFKL evolution following Mueller in |134j . 

Physical Picture 

Consider an ultrarelativistic gluon scattering on a target at rest. Before scattering on the target. 
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the gluon can emit one or more extra gluons. The emissions are illustrated in Fig. ^ In the infinite 
momentum frame considered here, the gluon moving in the light cone "+" direction has a very large 
typical light cone "+" component of its momentum, which we denote by p'^. The original gluon emits 
another gluon with a much smaller light cone momentum kf p^. A simple calculation shows that 
the emission probability is given by 

= ^ IT ^' 

where is the number of colors (and the SU{Nc) Casimir operator in the adjoint representation) and 
ki is the two-component transverse momentum vector of the gluon #1. Eq. (j^ can be obtained from 
the usual formula for photon Bremsstrahlung by multiplying it with the color Casimir Nc and replacing 
o^EM —>■ tts- Following |134j we assume that all emitted gluons have fixed transverse momentum of the 
order of some scale Q. This allows us to simplify the problem by replacing the transverse momentum 
integral d'^ki/k^ by a constant c. Defining the rapidity variable 

yi = ln|+ (5) 

we rewrite Eq. (jH) as 

a N 

dP, = c^^dy,. (6) 



Using Eq. (jH]) we conclude that given the rapidity interval 

Ayi = (7) 

the original gluon would split into two gluons with probability 1. 

To generalize the picture to many gluon emissions, we have to understand the space-time picture of 
the process. First we note that the typical light cone coherence time of the gluon #1 is 



Ti = X\ 



2 /c]^ 2 



Similarly, for the ith gluon in the gluon cascade, the light cone time is of the order of 

n « ^. (9) 

Thus, if the original gluon emits a cascade of gluons with their longitudinal momenta being progressively 
smaller 

p+ > fc^ > A;^ > . . . > fc]^ (10) 
than their light cone times would also be ordered 

ri > r2 > . . . > r^. (11) 

This cascade is shown in Fig. ^ There the gluons are ordered in time, such that the typical coherence 
time of each emitted gluon is much shorter than the coherence times of all the preexisting (harder) 
gluons. Therefore, for the ith gluon, all the gluons 1, ... ,z — 1 appear frozen in time. The gluon can 
be emitted off any of these preexisting gluons, which is shown by disconnected gluon lines in Fig. ^ 
Indeed, in the transverse direction the gluons are coherent only over short distances of the order of 
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Figure 1: Gluon cascade leading to the BFKL evolution. All gluon emissions are ordered in light cone 
"+" momenta, and, correspondingly, in light cone time. 



Ax± ~ l/Q, such that each gluon is emitted only by some fraction c' of the preexisting gluons. Since 
the colors of all these gluons are random, the ith gluon "sees" an effective color charge 



gi = Vdig (12) 

which results from a random walk in color space of i preexisting gluons. The probability of ith gluon 
emission is ^ 

dPi = cci— — - di/i. (13) 

71 

The rapidity interval required for ith gluon emission is 

TT 



= -77— ^- (14) 



Thus, the rapidity needed to emit gluons is given by 

N N 1 

so that the total number of gluons emitted in rapidity interval Y is given by 



This simple physical picture that we borrowed from [134" gives the right qualitative behavior of the 
BFKL evolution. The exact BFKL evolution also leads to the exponentially increasing number of glu- 
ons, just like we obtained in Eq. (|T6|) . The successive time ordered emissions discussed above and shown 
in Fig. ^resum the leading logarithms of energy, just like the exact BFKL equation: each power of as 
gets enhanced by a power of rapidity Y (which is equivalent to the logarithm of center of mass energy s), 
such that Eq. ()16|) resums powers of UgY . Even the kinematics of successive emissions considered above 
(see Eq. ()10|)) is the same as in the exact BFKL evolution, where extracting the leading logarithmic 
contribution requires ordering of the gluons' longitudinal momenta. 

The BFKL Equation 

Let us now consider a scattering of a bound heavy quark-antiquark state (quarkonium, or simply 
onium) on another quarkonium. The interaction between the onia is shown in Fig. |21 In general, we can 
write down the total onium-onium scattering cross section as a convolution of the onia light cone wave 
functions |135[ I136j with the imaginary part of the forward scattering amplitude of two quark-antiquark 
pairs (the imaginary part is denoted by F) 

atot = 2 Jd^x dz J d^x' dz' %g{x, z) $gg-(x', z') F{x,xl, Y), (17) 
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Figure 2: Scattering of two onia with the interaction mediated by the BFKL evolution. 



where ^qq{x, z) and z') are hght cone wave functions of the onia with transverse sizes x_ and 

x'. In the center of mass frame, each of the onia carries a large light cone momentum, p+ and p'~ 
correspondingly. For high energy eikonal scattering considered here, these two momenta (p"*" and p'^) 
are much larger than any transverse momentum scales in the problem. In each of the onia the quark 
carries the fraction z (or z') of the light cone component of the total momentum of the quark-antiquark 
state. Y = ln(s/4M^) is the rapidity variable, defined as a logarithm of the ratio of the center of mass 
energy of the onium-onium scattering s over some typical transverse momentum scale, like the onium 
mass M used here. 

At the lowest order the interaction in Fig. El is mediated by the two-gluon exchange. The corre- 
sponding imaginary part of the forward amplitude, which we denote hj F = F^^\ is given by 



a1 Cf 



(2 



—il-x 



il-x' 



Al-x' ^ 



where Cp = {N^ — l)/2A'^c is the SU{Nc) Casimir operator in the fundamental representation and / is 
the transverse momentum of each of the gluons. 

The BFKL evolution allows one to calculate quantum corrections to Eq. (fTH|) that bring in powers 
of agY. Apparently, as was shown in ^7\, such corrections preserve the two-gluon exchange structure 
of the interactions and can be summarized by the blob in Fig. |21 Corresponding generalization of the 
amplitude (fTH|) reads 



F(x,x',F) 



(2 - e-'^-^ - e'^-^) (2 - (T'^'^ - e'^-'^) f{l_, Y) 



Jl'-x' 



(19) 



with / and /' the gluons' transverse momenta on both sides of the blob as shown in Fig. |21 At the lowest 
(two-gluon) order the amplitude /(I, i',1^) is given by a delta-function 



f'^'\La = f{ii,Y = Y,) = 5\i-i:) 



(20) 



which, after substitution into Eq. (fT^ readily gives Eq. (fTB|) . {Yq is some initial rapidity, corresponding 
to the center of mass energy where the two-gluon exchange dominates the interaction.) 

The BFKL equation for the amplitude /(I, i',!^) with the initial condition reads ^2f] 



dfUXX) o^sN, 



dY 



d'^k 

(FTTp 



fikLY)- 



(21) 
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Figure 3: Schematic representation of a single rung of the BFKL evolution. 



The equation is illustrated in Fig. El where the dashed vertical line denotes the cut. The equation 
in Fig. El states that the blob can either consist of just a two gluon exchange without any evolution 
(the first term on the right hand side of Fig. EJ, corresponding to the initial condition in Eq. (1201), 
the blob may have small-x evolution corrections included. The evolution corrections can be real (the 
second term on the right in Fig. |21) and virtual (the third and fourth terms on the right of Fig. El). The 
real term contains a gluon in the final state (crossing the cut) and corresponds to the first term on 
the right hand side of Eq. (PT|) . The gluon is emitted off the t-channel gluons via the effective Lipatov 
vertices j57|, denoted by thick dots in Fig. El which represent the sum of all possible emissions. The 
virtual terms in Fig. El contain no gluon in the final state and correspond to the last term on the right 
of Eq. (EH). 




Figure 4: BFKL evolution ClS db ladder diagram. 

The solution of the BFKL equation should contain all iterations of the kernel depicted in Fig. El 
Iterating the BFKL kernel leads to the ladder diagrams, like the one shown in Fig. HI It depicts 
interaction of two onia characterized by typical transverse momentum scales Qi and Q2 interacting via 
a BFKL-evolved amplitude. Fig. |31 shows the structure of the blob in Fig. El in more detail. Again, 
the triple gluon vertices in Fig. 01 are not the usual QCD vertices: they are effective Lipatov vertices 
responsible for the real part of the BFKL kernel (27|. Similarly, the t-channel gluon lines do not 
correspond to the usual QCD gluon propagators: they give the so-called reggeized gluon propagators 

The emissions resummed by a Lipatov vertex are shown in Fig. El There we consider scattering 
of two ultrarelativistic quarks leading to production of a soft gluon with momentum <^ and 
k~ -C p'~ PH ESI 123 1211 EOl Ell- The correct emission amplitude, which is obtained by summing 
diagrams A-E in Fig. El can be written as the first diagram in Fig. El with the effective vertex triple 
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gluon given by [211 EH] 



in the 



p-q 




{k - qf 



+ k ,2q-k 



k- ' k^ 

form with T"" the SU{Nc) color matrix in the adjoint representation. 

C _ „ D 



(22) 




Figure 5: The diagrams contributing to effective Lipatov vertex, which is denoted by a very thick dot. 

To construct an effective reggeized gluon propagator one has to resum all leading logarithmic cor- 
rections to a single t-channel gluon exchange keeping the exchange amplitude in a color octet state. 
The details of this sophisticated resummation procedure can be found in |137j . Here we will only give 
the final result: the leading logarithmic corrections to the t-channel gluon propagator exponentiate, 
modifying the gluon propagator 



q'^ 

where a;(g^) is the gluon Regge trajectory [77| HU] 



(23) 



co{q') 



OsNc 

(2vr)2 



d^kq^ 



k\k-qy 



(24) 



and Ay is the rapidity interval spanned by a given t-channel gluon. 

Eqs. (|^ and give us the rules for vertices and t-channel propagators necessary to construct 
the ladder diagram in Fig. HJ There the leading logarithmic contribution is given by the multi-Regge 
kinematics of the produced s-channel gluons 



p+ > fc]^ > /c^ > . . . > A;^ > 

< /Ej; < < . . . < /C^V < V'~, 



and 



where ki 



k^ ~ ki 



k 



(25) 
(26) 

(27) 



\kA. 



Solution of the BFKL Equation 

To solve Eq. (pT|l we first have to find the eigenfunctions of its integral kernel. The kernel of 
BFKL equation is conformally invariant. It is easy to verify that a complete and orthogonal set of 
eigenfunctions of Eq. (PT|) is formed by the functions (23 1 



(2J 
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where (j) is the angle between vector I and some chosen axes and n is integer. The eigenvalues of the 
eigenfunctions in Eq. (PH|l are 123 1138j 

X{n,jy) (29) 

TT 

where 

X{n, v) = m - \i, (^i±M + _ (^1±M _ (30) 
with ^(z) = d\\].T{z) / dz. Denoting / ~ |/| and ~ |/'| we write the solution of Eq. as 

^ POO 

IXLLy) = / t^z^c„,,(F)ri+2^'^r^-2iv gi„(<^-<A')_ (3i^ 

Substituting Eq. into Eq. (|7T|) and using the eigenvalues from Eq. pU]) we find 

CnAY) = C^Se^>^("'^)^, (32) 
where the coefficient C^l is fixed by the initial conditions (pn|) giving 

^ni = (33) 

Combining Eqs. (IS2I), dSS)) and (jSH) yields 

Eq. (jMj) provides us the solution of Eq. (PT|) with the initial conditions given by Eq. (j^HI)- As one can 
see already from Eq. ()34|1 . the BFKL equation generates amplitudes which grow exponentially with 
rapidity Y. Remembering that F ~ Ins, this translates into a power of energy growth. 

Let us evaluate the amplitude in Eq. ()34|) a little further. Consider the case when / ~ /', i.e., the 
two momentum scales involved in the problem are not very much different from each other. A simple 
analysis of the function x(n, u) allows one to conclude that the dominant contribution to the amplitude 
is given by the n = term in the sum in Eq. (j34j) . Expanding x(n = 0, u) around the saddle point at 
z/ = we get 

X(0,z/) ^ 21n2-7C(3)z/^ (35) 

where ({z) is the Riemann zeta-function. Using Eq. (j35|) in Eq. (j34|) we can perform the //-integration 
obtaining j2ZI 



{ap-l)Y ' 



14C(3)a,F 

where we have defined the intercept of the perturbative pomeron 



(36) 



ap-1 = — ^ ln2 (37) 

TT 

and 

_ _ asN^ 

ots = • (38) 

TT 

The essential feature of Eq. (jHUj) is that it shows that cross sections mediated by the BFKL exchange 
grow as a power of energy 
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This behavior is reminiscent of the Pomeranchuk singularity in the reggeon calculus, and is, there- 
fore, sometimes referred to as the pomeron or the hard pomeron (to distinguish it from the soft non- 
perturbative interaction phenomenon). 



Problems of the BFKL Evolution 



The BFKL equation poses some important questions even in the case of heavy onium-onium scat- 
tering and at the leading order in as in the kernel. 

(i) The power of energy growth of the total cross section (Eq. ()39|1 ) violates Froissart unitarity 
bound, which states that the growth of the total cross sections with energy at asymptotically high 
energies is bounded by 

const ,9 , , , 

a < — In^ s 40 

ml 

with the pion's mass. (For a good pedagogical derivation of the Froissart bound we refer the readers 
to [nS]-) This implies that some new physical effects should modify the BFKL equation at very large s 
making the resulting amplitude unitary. 

(ii) The solution in Eq. (jMj) includes a diffusion term, which is the last term in its exponent. To 
see the potential danger of this term, let us consider a half of the ladder of Fig. HI stretching from one 
of the onia (the top one) to some intermediate gluon in the middle of the ladder carrying transverse 
momentum l± and having rapidity Y/2. Applying Eq. (jHUj) to that half- ladder we see that it includes 
a term 

exp [ ^^MiLl (41) 

14C(3)a. (r/2)J- ^^'^ 

This term is responsible for diffusion of the transverse momenta from the initial perturbative scale 
Qi both to high and low momenta, i.e., into infrared and ultraviolet. It implies that the distribution 
of gluons' transverse momentum in the ladder, while still centered around Qi, may have significant 
fluctuations towards high and low momenta / as shown in Fig. IHl where we plot the typical range of 
transverse momentum I in the ladder of Fig.EJas a function of rapidity of gluons in the ladder. The width 
of the diffusion grows with rapidity Y allowing for larger fluctuations at higher energies. Thus, no matter 
how large the starting scale Qi is, at certain very high energy the momentum of some gluons in the 
middle of the ladder would become of the order of Aqcd leading to the coupling constant as(AQcz)) ~ 1 
and thus invalidating further application of perturbative QCD (a^ <ti 1) and, consequently, of the BFKL 
evolution [22] • The diffusion starts from the scale Qi at one end of the ladder and from the scale Q2 at 
the other end (see Fig. IH)). At high energy the allowed momentum range broadens towards the middle 




Figure 6: Diffusion of the gluon transverse momentum in the BFKL ladder. Non-perturbative region 
is denoted by the scale A. 



13 



and touches the non-perturbative region of low momenta at mid-rapidity. The plot of allowed momenta 
as a function of rapidity shown here in Fig. El is sometimes referred to as Bartels cone (26j. Again it 
hints that the BFKL equation should be modified at higher energies to avoid the problem of running 
into the non-perturbative region. Alternatively, if such modification is not found within perturbation 
theory, we would be forced to admit that high energy asymptotics is dominated by non-perturbative 
physics. Fortunately this is not the will be shown below. 

2.1.2 The GLR-MQ Equation 

As we have seen above, the BFKL evolution leads to exponential growth of total cross sections with 
energy, violating the Froissart bound |25jj. It also leads to exponential growth of the density of partons 
in the onium (or hadron) wave function. To see this let us define an unintegrated gluon distribution of 
an onium by 

(j)ixBj,k^) = / ^i2-e-''--^-e''--^)f{Lk,Y = \nl/xBj)d^xdz<!>gg{x,z). (42) 

The definition (jl^ is illustrated in Fig. [7| To obtain it we have, essentially, truncated the lower onium 
in Fig. leaving two disconnected gluon lines. The diagram in Fig. [7| has to be calculated in A'^ = 
light cone gauge to give the gluon distribution function. The unintegrated gluon distribution function 
in Eq. gives us the number of gluons in the onium wave function having transverse momentum k 
and carrying the fraction Xbj of the onium "-|-" component of the momentum (Bjorken, or Feynman 
x). 




Figure 7: Unintegrated gluon distribution of an onium. 

Using the BFKL solution from Eq. ()36|) in Eq. one can easily see that the gluon distribution 

grows as 

/ 1 \"''"^ 

(f){xB,A^) ~ • (43) 

\XBjJ 

Therefore, the number of gluons rises sharply at small x / high energy, in agreement with the semi- 
quahtative estimate of Eq. (fTHI) . This feature is illustrated qualitatively by the gluon cascade repre- 
sentation of the BFKL evolution presented in the Sect. 12.1.11 There, the gluons are produced in a 
multi-Regge kinematics with comparable transverse momenta (see Eq. (I2Z|))- That means that the 
typical transverse sizes of the gluons, given by ~ 1/fc^ are also of the same order for all the gluons 

rj*- ~ ~ . . . ~ r]^. (44) 

Therefore, the BFKL cascade produces many gluons in the onium or hadron wave function, with roughly 
the same transverse size. As energy increases, more and more gluons are produced in the cascade. The 
gluons overlap in the transverse plane, creating areas of high gluon density. Thus, not only the number 
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Figure 8: BFKL equation as a high density machine. 



of gluons, but their density in the transverse plane increase with energy. This is illustrated in Fig. |H1 
for a wave function of a proton. At some initial value of Bjorken x, corresponding to lower energy, the 
proton's wave function may have one parton with transverse size ~ l/Q in it. As we go to smaller x, 
BFKL evolution would generate many more partons of comparable size, creating a region of high gluon 
density in the wave function |134j . 

However, the gluon density can not rise forever. It is known that in QCD the gluon fields can not 
be stronger than ^ 1/g ioi small coupling g. Therefore, when the gluon field reaches the density 
corresponding to field strength 

17 1 

(45) 



1 
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some new, possibly non-linear effects should become important slowing down the density growth 
This strong field constraint may be intimately connected to the problem of unitarization of cross sections 
at fixed impact parameter, i.e., to the black disk limit (see below). 

To understand how the growth of the gluon distribution can be tamed, Gribov, Levin and Ryskin 
(GLR) P3 I139j considered distribution functions of a "dense" proton or a nucleus. By "dense" proton 
we imply a model of a proton filled with sources of color charge — sea quarks and gluons, which 
were pre-created in the proton's wave function by some non-perturbative mechanism (see, e.g., |14Uj ). 
Gribov, Levin and Ryskin [TJ 1139] argued that for such systems multiple ladder exchanges may become 
important. Since one is interested in gluon distribution, which is a correlation function for two gluonic 
fields, these multiple ladders should come in as the so-called "fan" diagrams. An example of a fan 
diagram is shown in Fig. El There multiple BFKL ladders start from different quarks and gluons 
in the proton or nucleus, shown by straight lines at the bottom of Fig. 121 Due to high density of 
gluon fields, the ladders can not stay independent forever. As the energy increases so does the gluon 
density, eventually leading to recombination of the ladders, as shown in Fig. |H1 Ladder recombination is 
described by effective ladder merger vertices, denoted by blobs in Fig. jUl These vertices are sometimes 
called the triple pomeron vertices, since they connect three different ladders (BFKL pomerons). For 
their calculation we refer the reader to |141j and references therein. 

Gribov, Levin and Ryskin [H I139j suggested that, before the energy gets sufficiently high for all 
nonlinear effects to become important, there could be an intermediate energy region where the physics 
of gluon distributions is dominated by 2 — >■ 1 ladder recombination only. This recombination brought 
in a quadratic correction to the hnear BFKL equation, leading to the GLR evolution equation [H I139j 



d^{x,k^] 
dln(l/x) 



dH 

{k-iy 



k^(f){x,k^) 
f + (A; - ly 



ix.k' 



1 2 



(46) 
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Figure 9: An example of fan diagram resummed by the GLR-MQ equation. 



where, for simphcity, we assumed that the proton or nucleus has a shape of a cylinder oriented along 
the beam axis with the cross sectional area S± = ttB?. As expected, the linear term in Eq. pUj) 
is equivalent to the BFKL equation (j21|) . while the quadratic term, responsible for ladder mergers, 
introduces damping, thus slowing down the growth of the gluon distributions with energy. The growth 
of gluon distributions with energy given by Eq. ()46|) should slow down and saturate at very high energies. 
This phenomenon became known as the saturation of parton distributions. 

The ansatz pUj) of GLR |1, was proven by Mueller and Qiu |2I in the double leading logarithmic 
approximation (DLA) for the (integrated) gluon distribution functions, which are defined as [2] 



dk (h(x, k 



The double logarithmic approximation is a resummation of the powers of the parameter 

as ln(l/x) In(gVA'). 



(47) 



(48) 



The BFKL equation [27 was derived in the leading logarithmic approximation (LLA), corresponding 
to resummation of the parameter a^lnl/x, or, equivalently, agY (see Sect. I2.1.ip . In the limiting case 
of large in the distribution functions (or large k^ in the unintegrated distribution functions) another 
large logarithm becomes important: InQ^/A^, where A is the non-perturbative QCD scale. It becomes 
possible to define a new resummation parameter, Q;sln(l/a;) ln((5^/A^). In principle, the LLA is much 
broader than DLA: it resums powers of aghil/x with any dependence of the obtained terms on Q^, 
not restricting it to the leading logarithmic regime. Leading powers of asln((5^/A^) are, of course, 
resummed by the Dokshitzer, Gribov, Lipatov, Altarelli, Parisi (DGLAP) equation |142j . Indeed, the 
DLA limits of the BFKL and DGLAP equations are identical, since they are resumming the same 
parameter, asln{l/x) ln(Q^/A^). 

Employing the DLA and analyzing diagrams with two merging DGLAP ladders, Mueller and Qiu 
arrived at the following evolution equation [2J (again written here for a cylindrical nucleus) 



92 xG{x,Q^ 



91n(l/x)91n(g2/A2 



TT 



xG(x,Q^)-^ — [xG(x,Q 



2m2 



(49) 



which is known as GLR-MQ equation. Eq. ()49p is in agreement with Eq. (jmj) . and could be be obtained 
from the latter by taking the DLA limit and using Eq. ()47p . Thus Mueller and Qiu 2J proved Eq. ()46|) 
in the double logarithmic limit. 
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Eq. (jl^ allows one to estimate at which the non-linear saturation effects become important. 
To do that we have to equate the linear and quadratic terms on the right hand side of Eq. (j49p . The 
corresponding value of is called the saturation scale and is denoted by Ql- It is determined by 



The non-linear saturation effects are important for all Q ^ Qs, which is known as the saturation region. 

The quadratic damping term in both Eq. ()4(i|l and Eq. was believed to be important only near 
the border of the saturation region, for Q ^ Qs, where the non-linear effects were only starting to 
become important [H 12] • It was expected that higher order non-linear corrections would show up as 
one goes deeper into the saturation region towards Q < Qs (see, e.g., |143j ). In the next Section we 
will talk about the model where all such corrections could be resummed in a particular quasi-classical 
approximation, where the BFKL evolution can be neglected. 

2.2 Quasi- Classical Approximation 

As was suggested by the GLR equation, the effects of high gluonic field strengths come in as multiple 
exchanges. To better understand the effect of multiple rescatterings one can consider a specific model of 
hadronic or nuclear wave functions. In this section we will consider the case of deep inelastic scattering 
on a large dilute nucleus, where multiple rescatterings take place on individual nucleons and, as we will 
demonstrate, can be described by the Glauber-Mueller formula [HZj. In the infinite momentum frame 
this large nucleus can be thought of as a Lorentz-contracted "pancake" in the direction transverse 
to its velocity. Due to a large number of nucleons in the nucleus, the parton color charge density 
fluctuations in this "pancake" are large, described by a hard scale referred to as the saturation scale 
Qs. For very large nuclei, the saturation scale is large Qs S> ^qcd making the strong coupling constant 
small asiQs) ^ 1. At weak coupling the dynamics of field theories becomes classical: therefore, as was 
concluded by McLerran and Venugopalan 0, the gluon field of such a large ultrarelativistic nucleus 
should be given by the solution of the classical Yang-Mills equations of motion |144j with the nucleus 
providing the source current. The resulting small-x nuclear wave function systematically includes all 
multiple rescatterings and exhibits the effects of saturation. 

2.2.1 Glauber-Mueller Rescatterings 

Let us consider deep inelastic scattering (DIS) on a large nucleus. In DIS the incoming electron emits 
a virtual photon, which in turn interacts with the proton or nucleus. In the rest frame of the nucleus, 
the interaction can be thought of as virtual photon splitting into a quark-antiquark pair, which then 
interacts with the nucleus (see Fig. ITUI). Since the light cone lifetime of the qq pair is much longer than 



S±N, 



xG{x,Ql). 



(50) 




Figure 10: Deep inelastic scattering in the rest frame of the target. 
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nucleons in the nucleus 



Figure 11: Deep inelastic scattering in the quasi-classical Glauber-Mueller approximation in df^A^ = 
gauge. 

the size of the target nucleus, the total cross section for the virtual photon-nucleus scattering can be 
written as a convolution of the virtual photon's light cone wave function (the probability for it to split 
into a qq pair) with the forward scattering amphtude of a qq pair interacting with the nucleus 

aZ^'iQ^XBj) = I [<^t{x,z) + <^Lix,z)] d% N{x,b,Y) (51) 

with the help of the light-cone perturbation theory |1H5| lllffij . Here the incoming photon with vir- 
tuality Q splits into a quark-antiquark pair with the transverse separation x and the impact param- 
eter (transverse position of the center of mass of the qq pair) b. Y is the rapidity variable given by 
Y = Inlsx"^) ~ Inl/xBj- The square of the light cone wave function of qq fluctuations of a virtual 
photon is denoted by ^t{x,z) and for transverse and longitudinal photons correspondingly, 

with z being the fraction of the photon's longitudinal momentum carried by the quark. At the lowest 
order in electromagnetic coupling {a em) ^t{2L, z) and $l(x, z) are given by |145| 1130] 

<J>t(x, z) = ^ T.^em {«/ KfixTaf) [z^ + (1 - zf] + m}Ko{xTaff} , (52) 
f 

2N 

^l{x,z) = ^Y.^em ^Q''z\l-zf KlixTaj), (53) 
^ / 

with = Q'^z{l — z) + m'j, xt = |x| and X)/ denoting the sum over all relevant quark flavors with 
quark masses denoted by m^. a^^M ~ ej/47r with Cf the electric charge of a quark with flavor /. 

Our goal is to calculate the forward scattering amplitude of a quark-anti-quark dipole interacting 
with the nucleus, which is denoted by N{x,h,Y) in Eq. (jSH), including all multiple rescatterings of the 
dipole on the nucleons in the nucleus. To do this we need to construct a model of the target nucleus. 
Following Mueller [HZI we assume that the nucleons are dilutely distributed in the nucleus. Let us chose 
to perform calculations in the dfj_A^ = covariant gauge. There we can represent the dipole-nucleus 
interaction as a sequence of successive dipole-nucleon interactions, as shown in Fig. ^2 Since each 
nucleon is a color singlet, the lowest order dipole-nucleon interaction in the forward amplitude from 
Fig. ^2 is a two-gluon exchange. The exchanged gluon lines in Fig. ^2 are disconnected at the top: this 
denotes a summation over all possible connections of these gluon lines either to the quark or to the 
anti-quark lines in the incoming dipole. 

To demonstrate that the ordered gluon exchanges of Fig. are indeed the only possible interactions, 
let us consider a diagram where such ordering is broken, as shown in Fig. ^1 There, for simplicity, we 
consider scattering of a single quark on a pair of nucleons, denoted as #1 and #2, with interchanged 
ordering of quark-nucleon rescatterings. To show that such diagrams are suppressed in the dilute nucleus 
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Figure 12: An example of a diagram which does not contribute to the total dipole-nucleus cross section 
due to incorrect ordering of interactions with nucleons (x^ > xf). 

case, let us consider the scattering in the center of mass frame, where the projectile quark is moving 
in the light cone "+" direction and the nucleus is moving in the light cone "— " direction. In a dilute 
nucleus all the nucleons are spatially separated. For a nucleus moving in the light cone "— " direction 
this translates into a separation between the nucleons along the "+" axis: specifically, let us denote 
xf and X2 the light cone coordinates of the nucleons #1 and #2 in Fig. ^[ choosing X2 > xf for 
certainty. (As we will see, the opposite ordering will be indeed allowed.) We are interested in the 
integration over the /-momentum in the diagram of Fig. El More specifically, we are interested in 
integrating over the "— " component of this momentum. In the high energy scattering case considered 
here, the incoming projectile quark carries large momentum and the nucleons in the nucleus carry 
a large "— " component of momentum, e.g., p'~ carried by a quark in the nucleon #1. In such eikonal 
approximation the quark-gluon vertices do not depend on /, so that all the /-dependence is given by 
the gluon propagators for the / and q — I lines, and by the quark propagator of the k + q — I line. By 
requiring that the outgoing quark in the nucleon #1 is on mass-shell, (p' — /)^ = for massless quarks, 
we obtain = [p' — /)^ ~ —2p'~l~^, such that with the eikonal accuracy /"^ = 0. This means that the 
propagator of the /-line, given in general by — ^//^, becomes and does not depend on /~ anymore. 
Similarly one can show that the propagator for the q — /-line is i/{q — /)^ and is also independent of /~. 
Thus all the /"-dependence of the diagram in Fig. is given by the denominator of the propagator of 
the k + q — I quark line, since the numerator also becomes /"-independent in the eikonal limit. The 
relevant /~ integral can be written as 

r — -ii-{xt-xt) ^54^ 

7-00 {k + q-iy + te ■ 

Since k^ is very large, we rewrite Eq. (jM)) as 

r — -ii'{x+-x+) .ggN 

J-o.2k+{k- + q- -l-)-{k + q-l_y + te ' ^ ' 

Since fc^ > the integrand in Eq. has a pole in the upper half plane. However, as > Xi, the 
exponent in Eq. allows us to close the integration contour in the lower half plane only. Therefore 
expression in Eq. is zero, which means that the diagram in Fig. El can be neglected in the eikonal 
limit. (In fact, Eq. dictates which ordering of the nucleons is allowed: if < Xi the expression 
in Eq. ()55|) becomes non-zero, bringing us back to ordered picture of interactions shown in Fig. ITTl ) 

To resum the diagrams like the one shown in Fig.^^we start by calculating the graphs contributing 
to the scattering of a dipole on a single nucleon, as shown in Fig. El Remembering that we are working 
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Figure 13: Diagrams contributing to dipole-nucleon scattering at the two-gluon order. 

in the d^A^ = gauge and assuming for the moment that all scattering in the nucleon happens on one 
of its valence quarks we write for the sum of the diagrams in Fig. (similarly to Eq. (jl8|) ) 

"'^'^ ^ (56) 



where x is the transverse separation between the quark and the anti-quark in the incoming dipole, as 
shown in Fig. Instead of performing the integration in Eq. (j56p exactly, let us note that gluon are 
correct degrees of freedom responsible for the interactions only if the dipole is much smaller than the 
typical nucleon size, x_l ^ 1/A, with A some non-perturbative infrared cutoff. Expanding the argument 
of the integral in Eq. and integrating over / we obtain (see also jll6j ) 

"'"^"4 = In^. (57) 



As one can show by explicitly calculating the diagram in Fig. [7| without the evolution (or by using 
Eq. (j^mi in Eq. ()42|) . replacing qq dipole by a valence quark in the latter, and using the resulting unin- 
tegrated gluon distribution in Eq. (jlZ|)), at the lowest two-gluon order the gluon distribution function 
of a valence quark is given by 

Using Eq. ^ in Eq. dSH) yields 



x_^XBjG{xBjA/x_L). (59) 

Eq. (|^ can be easily generalized to the case of a nucleon by using the gluon distribution of a nucleon 
xGn in it instead of that of a valence quark 



2Nr 



X^XBjGNixBjA/Xx). (60) 



Eq. ()6()j) gives us the scattering cross section of a quark dipole on a single nucleon. The factor of 2 is 
due to the definition of total cross section in terms of the forward scattering amplitude, as shown in 
Eq. (ini). 

To obtain the multiple rescattering amplitude of Fig. ^2 we have to resum all multiple rescatterings 
of Fig. El (or, equivalently, Eq. (j60|l ). Defining the S'-matrix of a quark-anti- quark pair scattering on 
a nucleus at rest by S{z,x), where z is the longitudinal coordinate of the pair as it propagates through 
nuclear matter with z = at the edge of the nucleus, we can write the following equation 

d p cr'?'?^ 

■^S{z,x) = — S{z,x) (61) 
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with p the density of nucleons in the nucleus (p = y4/[(4/3)7ri?'^] for a spherical nucleus of radius R with 
atomic number A.) Eq. has the following physical meaning: as the dipole propagates through the 
nucleus it may encounter nucleons, with the probability of interaction per unit path length dz given 
by the product of nucleon density p and the interaction probability a'^^^ (2 from Eq. (jUUI) . (In the 
literature devoted to jet quenching in the medium this quantity is referred to as opacity |118| I119j .) 
The initial condition for Eq. ()61|) is given by a freely propagating dipole without interactions, such 
that S{z = 0,a;) = 1. We refer the interested reader to 1,116., .32 | I117j for a more detailed derivation of 
Eq. (jni}. Solving Eq. (jHT|l with S{z = 0,x) = 1 initial condition yields 

f p(T««"^ 1 

S{z,x) = e-x.Y>{ — z\. (62) 

To construct the forward scattering amplitude of a dipole scattering on a nucleus we need to know 
the S'-matrix for the dipole which went through the whole nucleus. To obtain it we need to put z = T{b) 
in Eq. ()62j). where T{b) is the nuclear profile function equal to the length of the nuclear medium at a 



given impact parameter b, such that T{b) =2\J E? for a spherical nucleus of radius R. Identifying 
N{x,b, Y) with the scattering T-matrix, such that = 1 — S" we finally write jHZj 

N{x,b,Y = 0) = l-expl-^-^T{b)\ = 1 - exp j pT(6) xl xb,Gjv(xs,-, 1/4) [ • (63) 

This expression is known as Glauber-Mueller formula EZj , since it presents a generalization of the 
Glauber model of independent multiple rescatterings in the nucleus jHEI to the case of QCD [flTj . 

We put y = in the argument of in Eq. ()63|) to underline that this expression does not include 
any small-x evolution which would bring in the rapidity dependence. Indeed, Eq. ()63p was derived in 
the approximation where the interaction of the dipole with each of the nucleons is limited to a two- 
gluon exchange. At this order the gluon distribution xbjGn is given by Eq. (f^ and is, therefore, 
Xsj-independent, leading to rapidity- independence of the whole expression in Eq. (jHS))- 

Eq. ()63|) allows us to determine the parameter corresponding to the resummation of the diagrams 
like the one shown in Fig. ^2 Using the gluon distribution from Eq. (jSHjl in Eq. (jUHj) . and noting that 
for large nuclei the profile function scales as T(b) ~ A^^^ and the nucleon density scales as p ~ we 
conclude that the resummation parameter of multiple rescatterings is jS] 



a 



A^/\ (64) 



The physical meaning of the parameter A^^^ is rather straightforward: at a given impact parameter 
the dipole interacts with ~ A^^^ nucleons exchanging two gluons with each. Since the two-gluon 
exchange is parametrically of the order we obtain A^^^ as the resummation parameter. 



Defining the quark saturation scale 



QM - i^pT® (65) 



we rewrite Eq. (|?K?jl as 



iVfci.y = 0) = l-exp^-^i«®Ml/^^|. (66) 



I 



The dipole amplitude A^ from Eq. (jUBj) is plotted (schematically) in Fig. as a function of x±. One 
can see that, at small x±, x± <^ ^/Qs, we have N ^ x^ and the amplitude is a rising function of x±. 
However, at large dipole sizes x±>l/Qs, the growth stops and the amplitude levels off (saturates) at 
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Figure 14: The forward amplitude of the dipole-nucleus scattering plotted as a function of the 
transverse separation between the quark and the anti-quark in a dipole x± using Eq. (jU^ . 

= 1. This regime corresponds to the black disk limit for the dipole-nucleus scattering: for large 
dipoles the nucleus appears as a black disk. To understand that = 1 regime corresponds to the black 
disk limit let us note that the total dipole-nucleus scattering cross section is given by 

atl^ = 2 Jd%N{x,b,Y) (67) 

where the integration goes over the cross sectional area of the nucleus. If = 1 at all impact parameters 
b inside the nucleus, Eq. flUTj) gives for a spherical nucleus of radius R 

a!!f = 2nR\ (68) 

which is a well-known formula for the cross section for a particle scattering on a black sphere |146j . 

The transition between the A^ ~ to A^ = 1 behaviors in Fig. 1141 happens at around x_l ~ 1/Qs- 
For dipole sizes a;_L> l/Q^ the amplitude A^ saturates to a constant. This translates into saturation of 
quark distribution functions in the nucleus, as was shown in (as xg + ~ F2 ~ 0-^*1^), and thus 
can be identified with parton saturation, justifying the name of the saturation scale. We will show this 
connection between saturation for the amplitude A^ and for the partonic wave functions using gluons 
as an example in Sect. 12.2.21 below. 

Before we proceed let us finally note that, since T(b) ~ A^^^, the saturation scale in Eq. scales 

as pinij 

~ A^/^ (69) 

with the nuclear atomic number. Eq. ()69j) implies that for a very large nucleus the saturation scale 
would become very large, much larger than Aqcd- If Qs ^ ^qcd, the transition to the black disk limit 
in Fig. happens at the momentum scales (corresponding to inverse dipole sizes) where the physics is 
perturbative and gluons are correct degrees of freedom. Therefore, Eq. (jU^ is of paramount importance, 
since it justifies the approximation we made throughout this Section that dipole-nucleon interactions 
can be described by perturbative gluon exchanges instead of some non-perturbative mechanisms. 
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2.2.2 The McLerran-Venugopalan Model 



Point Charges Approach 

Let us consider a large ultrarelativistic nucleus in the infinite momentum frame. We are interested 
in the small-x tail of the gluon wave function in the nucleus. In the rest frame of the nucleus the small-x 
gluons have coherence length of the order of [TU] 

hoh ~ TT^ (70) 

ZrriN XBj 

where m^v is the mass of a nucleon. If xbj is sufficiently small, the coherence length may become 
very long, much longer than the size of the nucleus. Such small-x gluons would be produced by the 
whole nucleus coherently (only in the longitudinal direction). An example of such interaction is shown 
in Fig. El There the small-x gluon (wavy line) interacts coherently with several Lorentz-contracted 
nucleons. Indeed the nucleons, and the nucleus as a whole, are color-neutral and one may worry that a 
coherent gluon simply would not "see" them. However, the gluon is coherent only in the longitudinal 
direction: in the transverse direction the gluon is localized on the scale x± ^ l/kx with the transverse 
momentum of the gluon. If kr ^ ^qcd, which is a necessary condition for using gluon degrees of 
freedom, the transverse extent of the gluon would be much smaller than the sizes of the nucleons. 
Because of that the gluon would interact only with a part of each nucleon in the transverse direction, 
as shown in Fig. ^1 The color charge in that segment of the nucleon that a gluon is traversing does 
not have to be zero: the gluon may run into, say, a single valence quark there. As a result of such 
interactions, the gluon would "feel" some effective color charge of all the nucleons' segments that it 
would traverse through. In the spirit of Glauber approximation we may assume that all nucleons are 
independent of each other, so that interactions with parts of different nucleons are similar to a random 
walk in color space. If each individual nucleons' segment has a typical color charge g, than, due to the 
random walk nature of the process, the total color charge seen by the gluon at a fixed impact parameter 
would be g y/n, where n ~ A^^^ is the number of nucleons at this impact parameter. 




Figure 15: Small-x gluon sees the whole nucleus coherently in the longitudinal direction and interacts 
with several different nucleons in it. The effective color charge seen by the gluon is a result of a random 
walk in the color space (see text). 

In the infinite momentum frame, due to Lorentz contraction, all the nucleons appear to be squeezed 
into a thin "pancake" of Lorentz-contracted nucleus, as shown on the right of Fig. El One may then 
define an effective color charge density seen by a gluon in the transverse plane of the nucleus IHHEl IH EHj • 
The typical magnitude of these color charge density fluctuations is given by the color charge squared 
divided by the transverse area of the nucleus, {g y/nY / S± = g'^n/S±. The number of color charge 
sources in the whole nucleus is proportional to the number of nucleons in the nucleus, n ~ A. The 
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typical color charge density fluctuations are, therefore, characterized by the momentum scale 

^ A^K (71) 

It is important to notice that the momentum scale in Eq. (j7ip grows with A as A^/^, similar to the 
saturation scale in Eq. (j69p . The important conclusion we can draw from Eq. (j7ip is that for sufficiently 
large nuclei their small-x wave functions are characterized by a hard momentum scale /i, which is much 
larger than Aqcd allowing for a small coupling as description of the process. Field theories with small 
coupling are usually dominated by classical fields, with quantum corrections bringing in extra powers 
of the small coupling constant a^. Therefore the dominant small-x gluon field of a large nucleus is 
classical. It can be found by solving the classical Yang-Mills equations of motion |144j 

V^Ff"" = r (72) 

with the nucleus providing the source current J*^, such that, in the infinite momentum frame 

r = 6''+pix,x-) (73) 

with p{x, x^) the color charge density. This conclusion is originally due to McLerran and Venugopalan 
[2]. Another way of understanding why the classical field dominates is by arguing that large color charge 
density fi implies high occupation numbers for the color charges in the nuclei: high occupation numbers 
lead to classical description of the relevant physics. 

We need to find the gluon field of the nucleus in order to construct its gluon distribution function. 
The latter is most easily related to the gluon field in the light cone gauge of the nucleus. However, 
the classical gluon field of a nucleus is easier to find in the covariant d^j.A'^ = gauge. To do this we 
will, for simplicity, assume that all the relevant color charge in the nucleus is carried by the valence 
quarks. Furthermore, we will specifically chose to consider a nucleus with "mesonic" nucleons made 
out of qq pairs instead of three valence quarks [3]. (This latter assumption would only simplify the 
calculations, with the conclusions being easily generalizable to the case of real nuclei.) Considering 
the nucleus moving ultrarelativistically in the light cone + direction, we label the "valence" quark and 
anti-quark coordinates by Xi,x~ and x!i,x'~. In the recoilless eikonal approximation considered here 
neither one of these coordinates changes due to "emission" of the gluon fields ^3J. In covariant gauge 
the color charge density of such ultrarelativistic nucleus made of valence quarks is given by 

A 

P^ix, x-)=gY: (tt) [Six- - X-) 5{x - - 5{x- - x'-) 5{x - x',)] (74) 

with (t°) the color charge matrix of the quark and the anti-quark in the zth "nucleon" in fundamental 
representation. The classical electromagnetic field of an ultrarelativistic point charge e at a; = 0, x~ = 
in d^A^ = Q gauge can be easily found |147^ Hj 

A'+(x) = -—5{x-) H\x\k) (75) 

271 

with all other field components being zero and A some infrared cutoff. Since this field is localized by the 
delta-function along the light cone, the fields of a ensemble of charges located at different longitudinal 
coordinates x~ do not overlap. Therefore, if we want to generalize Eq. ()75|) to the non-Abelian case 
we need not worry about the non-Abelian effects due to the overlap of the fields from different charges, 
since those do not take place. The non-Abelian generalization can be easily accomplished by replacing 
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the electric charge by its QCD analogue, e —>■ git"'), such that the gluon field of the nucleus in covariant 
gauge is given by ^ 



A 



^ E E^"(^i) k(x--xr)ln(|x-x,|A)-5(x--xr)ln(|x-x:|A)|,A' = 0,A'- = 0. (76) 

a=l i=l 



A 



As one can explicitly verify, Eq. (fTUI) satisfies Eq. (f7^ with the source current given by Eq. (f7H|) with 
the color charge density from Eq. ()74|) . We now have to perform a gauge transformation of this field 
into the light-cone gauge. The field in a new gauge is 



= SA^S-^ - -{d,S)S-\ 
Requiring the new gauge to be the light-cone gauge, = 0, we obtain 

S{x,x^) = Pexp ^—ig J dx'^ A^{x,x'^)^ . 
The field in = light cone gauge is 

A{x,x-)= r dx'-F+^{x,x'~) = r dx'~S{x,x'-)F'+^{x,x'-)S~\x,x' 

J—oo J — oo 



(77) 



(78) 



(79) 



Using the explicit expression for the field in covariant gauge from Eq. ()76|) in Eqs. ()79p and ()78|) we 
obtain ^ 



= ^ E E(^n {s{x,xirs-\x,x-)-^-^^ 

a=l i=l \ l-L Xi\ 



[X — X, 



- S{x,x'r)eS-\x,x'r)^^e{x- -x'r)] (80) 

\x — XiY ) 

with 

A 

S{x, x~) = Y[ exp 
1=1 

Eqs. (lHn|l and (jHT|) give us the classical field of a large ultrarelativistic nucleus moving in the light 
cone + direction calculated in the light cone gauge of the nucleus A'^ = 0. We will also refer to this 
field as a non-Abelian Weizsacker- Williams field, since it is a non-Abelian analogue of the well-known 
Weizsacker- Williams field in electrodynamics |148j . Below we will rederive this result in the approach 
where the nucleus is described by a continuous color charge density p(x) and will finally calculate the 
correlator of two of the fields in Eq. (jHn|l to obtain the classical gluon distribution function of a nucleus. 

Let us first identify which Feynman diagrams correspond to the non-Abelian Weizsacker- Williams 
field of Eq. ()8U|) . Since a detailed analysis of the problem can be found in [5 , here we will only quote the 
answer. An example of the diagram contributing to the classical field of Eq. (jHn|) is shown in Fig. ^1 
Classical fields are usually given by tree diagrams (graphs without loops): this is indeed the case in 
Fig. where the gluons from different sources (valence quarks in nucleons) keep merging with each 
other until they form a gluon field at the point where we "measure" it, which is denoted by a cross in 
Fig. Uni Each nucleon can only emit one (inelastic) or two (elastic) gluons: as discussed in more detail 



• 2 /I _ |\ 

1^ E'-W) In (1^) 
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Figure 16: Classical gluon field of an ultrarelativistic large nucleus: the non-Abelian Weizsacker- 
Williams field. 

in 12] emission of more than two gluons per nucleon allows for diagrams with quantum loops to be of 
the same order as the classical fields. Such graphs has to be discarded because the classical fields does 
not give the dominant contribution at that order. The two-gluon per nucleon limit is the same as the 
one considered in Sect. 12.2.11 giving rise to the resummation parameter for multiple rescatterings from 
Eq. ()64j) . As we will see below classical field from Eq. (jHIH) resums powers of the same parameter a^A^^^. 
The discarded diagrams having more than two gluons per nucleon are suppressed by extra powers of as 
brought in by the extra gluons, which are not enhanced by powers of A. 

Effective Action Approach 

In order to investigate high gluon density effects, McLerran and Venugopalan proposed an effective 
action for QCD at small x and/or for large nuclei (the MV model) 0|. To illustrate the approach, it is 
easier to consider a nucleus (this can be easily generalized to a proton) in a frame where it is traveling 
with the speed of light so that it has a large (light cone momentum). In this frame, the longitudinal 
size of nucleus is Lorentz contracted {2R —>■ 2R/'~f) so that all the valence quarks are localized in a 
longitudinal "pancake" of width 2R/'~f. Furthermore, consider the transverse area of the nucleus as a 
rectangular grid of size r -C 1/Aqcd so that we are not sensitive to details of confinement. Since the 
nucleus is taken to be large, there is a large number of valence quarks (and large x gluons) in this 
rectangular "pancake" . As we discussed above, the small x gluons have large longitudinal wavelengths 
and therefore can not resolve individual valence quarks in the nucleus. Rather, they couple coherently 
to the effective color charge generated by the valence quarks and large x gluons, denoted by p{x~,x). 
This color charge is taken to be independent of the light cone time since the small-x gluons, to which it 
couples, have very short life times (Ax"*" ~ l/xP'^) and see the effective color charge as frozen in time. 

Based on the above physical picture, one can write the following effective action for QCD at small x 
and in the light cone gauge A'^ = OE] (for an alternative but equivalent form of the action see |149j ) 

S = ~J d^xF^^iyx) + j cPxdx~ Trp{x~ ,x)W_oo,oo[A']{x~ ,x) +i J d'^x dx~ F[p{x,x')] (82) 



where 



W_oo,oo[A ]{x ,x) = Pexp 



oo 



ig / dx A^{ 



^3) 



and exp[i{i J d'^xdx F[p{x,x )]}] can be thought as the statistical weight of a given color charge 
configuration. In the original MV model, F was taken to be a Gaussian, given by [HI IH EHl 

J dx-F[p{x,x-)] = ^p\x) (84) 
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while the interaction term was taken to be S{x') p{x)A^ . It turns out that as long as one is interested in 
the evolution of observables with x, it is valid to take p{x, x~) —>■ p{x)6{x~) with logarithmic accuracy. 
However, if one is interested in calculating observables at the same point in x~ , one needs to take the 
structure of the color charge density p in the longitudinal direction x~ into account. 

In order to calculate a physical observable, one first solves the classical equations of motion and 
then averages over the color charge configurations p with the weight function F[p]. 

J[Dp]D[A] exp{-F[p]} 0(A)e^^[^'^l 



< 0{A) >E 



![Dp]D[A] exp{-F[p]}e^^I'''^l 



5) 



As we have just mentioned, in the original MV model, the weight function F is taken to be a 
Gaussian. This is a good approximation as long as the number of color sources is large and if the color 
sources are not correlated, such as in a large nucleus |4i . However, evolution in x will change the weight 
function and it will not be a Gaussian in general. 

The classical equation of motion involving the non-zero component of the current {J~^) is 



d,d+A' + -tg[A,d+A']=gJ' 



and Fij = 0. The solution has the form A* = 6{x~)a^{x) where ai{x) satisfies diai{x) = gp{x)- In the 
original work of McLerran and Venugopalan [3^ , it was argued that the commutator term is zero since 
it involves the commutator of the fields at the same point in x~ . This is not quite correct however, 
and leads to infrared singular solutions. The reason is that this term is very singular due to the presence 
of the so that one needs to understand the structure of the field across this singularity, which in 

turn means that one needs to know the structure of the color source distribution p in x~. 

In pi, 1^! the structure of the color sources in x~ was taken into account. Following 6j the equation 
of motion now becomes 

Di^A\x,y)=gp{y,x) (87) 

where we have introduced the space-time rapidity variable y ~ lnl/a;~. Furthermore, the weight 
function is also modified in order to take the rapidity dependence of the sources into account and is 
given by 



D [p] exp 



dy / (fx 



To find the classical solution, we introduce the path ordered Wilson line, given by 



U{y,x) = Pexp 



dy'A{y\x) 



^9) 



Since F^ = 0, the classical field Ai must be a pure gauge in two dimensions so that it can be written as 



A'{y,x) = -U{y,x)d'U\y,x), 



(90) 



which, substituted into (jHTj) yields the following equation for A: d^A = —g'^Wp U. The classical solution 
can be written as 



A%y,x) 



9 Jy 



dy'Uoo,y'{x) [d'A{y',x)] t/y,oo(^)- 



(91) 



Furthermore, since the transformation between p and A does not depend on p, one can perform the 
color averaging using A. In other words, for any operator O, we have the following relation 



D[p] 0{p) exp 



n 

dy / d'^x 



p'^{y,x) 
p^{y) 



D[A]0{A) exp 



- r ^^ Jd^A{y,x)f 

J g^p^ijj) 



(92) 
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Comparing Eq. (j^T|) with Eq. ([7^1) we can identify matrices U and S and the function A/x~ with the 
gluon field A'~^ in covariant gauge, which, in the case of point charges approach, is given by Eq. (j76p . 
Therefore, the solution from Eq. (PT|) is equivalent to the solution found previously in Eq. (jHUI). 

Classical Gluon Distribution 



Now we can compute the correlator of two gluon fields 



G^^{y,x,y) ^ {A1{y,x)A%y,y)) 

which is related to the unintegrated gluon distribution via 

1 



lWW, 



2 7r2 



d'bd're-'^-^ Tr(A(0) ■ A(r)), 



(93) 



(94) 



where the index WW denotes the quasi-classical Weizsacker- Williams distribution function and A{x) = 
hmy^_^A{y,x). 

The easiest way to compute the correlator (jUHj) is by expanding the path ordered exponentials, 
performing the color averaging and exponentiating the result. The first term of the expansion gives 



1 /■°° 1 
Gijfliy, x; y,y) = — 5ab J dy'g^ fi^{y')didj-^{x, y). 



(95) 



The inverse of the operator is infrared singular and must be regulated. A natural cutoff is the 
QCD confinement scale Aqcd- A more refined treatment shows that the saturation scale Qg provides a 
natural infrared cutoff for the operator Nevertheless, we define 



^(x) = 7(x) = \ii[x'^t^Icd] + 7(0) 



so that the next term in the expansion gives 



%i(?/,^,y) = — — 



Similarly, the n^^ term in the expansion gives 



Ay" 



[didj-f{x - y)] [7(x - y) - i 



G 



ab 



Jab 
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Ay" 



n+l 



didj-i{x- y)Yi{x- y) -7 



We can now sum the series and get the following expression for the gluon distribution function 

4(iV2-l) 1 



Gt:{y,XT) 



gm, xl 



(96) 



(97) 



(9J 



(99) 



where we have defined x{y) = dy' ji^ {y') . Using Eq. (fTIj) to calculate a correlator of two p's in light 
cone gauge we can identify x{y) = pT{b)/4Nc j4j. Alternatively the correlator in Eq. (jUHj) can be 
re-calculated using the gluon field from Eq. (IHUj) In the end one can rewrite Eq. as 



Gl^iy,XT) 



2Cp 1 



as rc xj. 



1 — exp 



X. 



2 Ql ln( l/xTAQCD) j 



(100) 



which looks very similar to the Glauber-Mueller formula ()66|). except that the saturation scale in 
Eq. ()100p is now for gluons, 

Qlib) ^ AiralpTib), (101) 
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Figure 17: Unintegrated gluon distribution in the quasi-classical approximation as a function of gluons' 
transverse momentum (see Eq. (|1UU|) ). 

and is different from the quark saturation scale ()65|) by replacing Casimir operators, Cp —>■ N^. 
Combining Eqs. ()94|) . (j^Hj) and (jlOOj) we obtain the following gluon distribution function 

^^^(XB,,^') = / d%d're-^!^-^-^NG{r,b,Y = 0) (102) 

as 2 7T-^ J r\_ 

expressed in terms of the gluon dipole-nucleus forward scattering amplitude 

Na{rAY = 0) = 1 - exp HI/tt^CCo) ^ ^,03) 

Note that the unintegrated gluon distribution in Eq. ()102j) is proportional to l/a^ when Nq ~ 1: this 
corresponds to gluon fields being Ai ~ 1/g characteristic of all classical solutions. 

It should be noted that is an all twist result for the gluon distribution function, valid in the 
classical regime, which resums multiple scattering of gluons from the target hadron or nucleus, in the 
spirit of the Glauber-Mueller formalism and can be thought of as the initial condition for an evolution 
equation which would take gluon emission into account, such as the JIMWLK equation which we will 
consider below. 

It is instructive to look at the different limits of ()99|) . or, equivalently, Eq. ()102|) . In the limit 
where ^ (perturbative QCD limit), we get G{xt) ~ In(xT) so that in momentum space using 
Eq. (j94j) we have (pixsj,!^) ~ l/k"^ in agreement with pQCD. On the other hand, in the limit where 
xt is large (but is smaller than 1/Aqcd), we get G{xt) ~ l/x|. so that in momentum space we have 
4>{xBj,!}^) ~ In^kxy. The momentum distribution of gluons is shown in Fig. [T7I where the slowing 
down of the infrared divergence of the gluon distribution function is evident. 

It is important to notice that the gluon distribution function (solid line in Fig. IT7|) is now only 
logarithmically infrared divergent, in contrast to pQCD (the first term in expansion of Eq. (|l(J(jp in 
powers of Qg shown by dashed line in Fig. I17|). which would lead to power divergence. In other words, 
the non-linearities of the classical gluon field have (almost) regulated the infrared divergence present 
in leading twist pQCD, making the residual singularity integrable over k. This phenomenon is usually 
referred to as gluon saturation: the increasing gluon distribution slows down its growth in the infrared 
from the power-law scaling to the logarithmic one. It is interesting to note that the saturation scale 
characterizing unitarization of the dipole-nucleus forward scattering amplitude N in Sect. 12.2. II turns 
out to be the same saturation scale (modulo the Casimir factors) as the one governing the saturation 
of gluon distribution functions in Fig. [T7| 
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Figure 18: Phase space distribution of gluons in the transverse momentum space. 



The phase space distribution of gluons is shown in Fig. ^[ where we multiphed the unintegrated 
gluon distribution (l>{xBj,k. ) from Fig. Elby the phase space factor kx- In contrast with leading twist 
pQCD (dashed line), the full classical line in Fig. ^1 has a peak, which is due to the fact that the 
infrared divergence is now regularized. The location of the peak determines the typical momentum of 
the gluons in the small- a; hadron/nucleus wavefunction and is given by the saturation momentum since 
it is the only scale in the problem. The fact that most gluons in the wave function of a hadron or nucleus 
reside in a state with a finite momentum Qg is the reason why this state is also called a condensate. 
It allows us to treat the gluon wave function of a high energy hadron or nucleus in the small coupling 
approximation. 



2.3 Quantum Evolution 

It is known that quantum corrections to the classical results presented in the previous section are 
potentially large. These large corrections arise from emission of gluons, both real and virtual. If the 
available energy is sufficiently high, it leads to a large phase space for gluon emission which gives rise to 
large logs of energy (or of which need to be resumed. In the CGC formalism, the presence of these 
large corrections to the classical approximation was proved in jlSUt I151j . A Wilson renormalization 
group formalism was developed in [71 |H] which allows for resummation of these large logs. The resulting 
equation can be written as an evolution equation for correlator of any number of Wilson lines and is 
known as the JIMWLK equation. This (functional) equation was later shown to be equivalent to a set 
of coupled equations derived by Balitsky [22] • In the large- A'c limit, these equations can be written in a 
closed form, as a single integral equation, independently derived by one of the authors for the correlator 
of two Wilson lines [21], known as the BK equation. 



2.3.1 The JIMWLK Equation 

The starting point is the action given by ()82|). We introduce the following decomposition of the full 
field A^, 

A^ = b^ + 6Af, + (104) 

where 6^ is the solution to the classical equation of motion considered in the previous section, 
denotes the quantum fiuctuations having longitudinal momenta g"*" between P+ and Pn-i, and is 
a soft field having longitudinal momentum where < . The effective action is calculated for 
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field a^, integrating out hard fluctuations 5^4^, assuming that these fluctuations are much smaller than 
the classical field b^. This procedure reproduces the form of the action given in (j82j) with a modified 
functional F', due to inclusion of the hard fluctuations into the color sources. The change of the 
functional F with logl/x leads to an evolution equation for the statistical weight functional, known as 
the JIMWLK equation, which can then be used to derive evolution equations for any n point function 
of the effective theory. 

Expanding the action around the classical solution and keeping first and second order terms in hard 
fluctuations (5^4^, we get 



S 



^G{af - l-6A,[D-\p)rSA, + ga-p' + 0((a-)2) + tF[p] 



4 



where 



with 



p = p + 6pi + 5p2 

-2r'"'a^i{x)SA^{x,x+,x- = 0) 

|/"V(2l) / dy^ e{y+ - x+) - e{x+ - y+) 



(105) 

(106) 
(107) 



dA'^x.y'^.x' = 0) 



and 



5p2{x,x~^,x ) 



= -r'"'[d+5A'l{x, a;+, x-)]6A'r{x, x+, x") 

- I^p'^lai) / dy+6A-%x,y+,x' = 0) J dz+6A-'^{x,y+ ,x- = 0) 

X 0{z+ - y+)0{y+ - x+)tiT''T''T^T^ 

+ e{x+ - z+)e{z+ - y+)iiT^T^T^T'^ 

+ ^(2+ - x+)^(x+ - y+)tTT^T'^T^r 



(108) 



The first term in both 6p1 and 6p2 is coming from expansion of in the action while the rest of the 
terms proportional to p are from the expansion of the Wilson line term. The three terms correspond 
to different time ordering of the fields. Since the longitudinal momentum of a~ is much lower than 
of 6 A, we have only kept the eikonal coupling which gives the leading contribution in this kinematics. 
The higher order terms in a~ do not affect the derivation here and can be reconstructed later using the 
uniqueness of the action. 

At first stage, we integrate over 6Afj_ while keeping both p and 6p fixed. We note that it is sufficient 
to keep only the first and second order terms in correlation functions of 6p, all higher order terms being 
suppressed by powers of as- We therefore define 



asya [x) 

,ab/ 



<5p''{x,x+)5p\y,x+)>sA = asyx^'U,]/) 
with y = ln(l/x). Integrating out the hard fluctuations 6A^ then leads to 

' D[p,p'][Det{x)r'^'exp{-F[p]) 



(109) 



X exp 



D[p,p]exp{-U[p,p]} 



[p'x - Px-OtsV CTx] [x. 



xy J 



p - Py- asyay 



(110) 
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where an integration (summation) over repeated indices is implied. In the next step, we integrate out 
p at fixed p' . This can be done using the steepest decent method since the integrand is a steep function 
of p, peaked around p' . The steepest decent equation is given by 



Px- Px- Oisycrx = Oisy Xxu 



5F_ 

Spu 



Substitution of this into Eq. piO|) gives 



2 dpu 



u 



F + -trln(x) 



(112) 



5pu 2 5pu 



5pv 2 5p^' 



In the above expression all the functionals are taken at p*^ which is the solution of the steepest descent 
equation (|111|) . We also need the correction due to quadratic fiuctuations of p around the steepest 
decent solution po- It is given by 

5^U 



5pJPy 



r]^""^ 5pjpy 
5a„ 5a, 



F + -trln(x) 



(113) 
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Furthermore, there are contributions from the third and the fourth derivatives of U, explicit expressions 
for which are long and can be found in [H]. Putting everything together, we recover the form of the 
action given in (jH^ with the new functional F' given by 



F' = F + 



5^ 



X^ 



SpuSpv 



-F- 



+2 



6F 5x. 
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which can be rewritten as an evolution equation for the functional F in terms of the one and two point 
fiuctuations a and x as 

6^ ^ 6^ Xuv 



dy ~ 2 



Xv 



+2 



Spjpv 
5F 5x. 



Spu 5py 



6F 6F 

r- r- ~ T XuV~? 

OpuOpy opu Op^ 
5au 6F 
dPu dpu 



(115) 



If we define the statistical weight functional Z = e ^ , the above equation takes a very simple form if 
written as an evolution equation for Z known as the JIMWLK equation 

\-^Xuv\ - \Z(J^ 
2 dpuOpv opu J 



d 



dy 



a. 



(116) 



The equation (jll6|) can be used to derive evolution equations for any number of correlators of the 
color charge density p. As an example, we consider the two point function. Multiplying both sides of 
flll6|) with px py and integrating by parts over p gives 

d 

-r- < PuPv >= as < Xuv + PuCTv + PvCTu > ■ (117) 

dy 
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Figure 19: Diagrams contributing to the color charge fluctuation %. 

This equation can be shown to reduce to the BFKL equation in the low density limit which will be 
considered later. 

To complete the derivation of the JIMWLK equation, we need to calculate the one and two point 
fluctuations x ^^id a in terms of the color charge density p. The details are given in 0-112] • Here we 
follow the derivation in JH] which is done in coordinate space and is closer in spirit to the derivation 
of BK equation in the next section. The starting point is x defined in ()109p to which only the color 
charge pi given in ()107|1 contributes. Suppressing the overall coupling constant and rapidity for the 
moment, x is 



Xabix, y) 



id~ 



■6A' 



2% 5A'r 



+i 



5A- 



-P 
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where = 5{x )ai{x) (integration over x and y is implied above). Eq. (jll8|l can be rewritten in 
terms of the propagator of the hard fluctuations 5 A as 

1 , , 



Xah{x, y) 



1 



2J^+' + 2J^+' {x\G' 



id' 



Py 



P. — G-%)2^+* + ip,{x\—G- 



1 



Py 



(119) 



id~ ^ id~ id~ 

The four terms in ()119|) are depicted in Fig. ()19|) and correspond to real corrections to the evolution 
equations. The dashed lines denote the propagator of the hard modes in the background field while 
the dash-dotted lines denote the soft modes for which the effective theory is written. The thick wavy 
lines attached to a circle denote the classical field while the thick wavy lines attached to the solid lines 
represent the color charge density p. 

The hard fluctuation propagators G^" were derived in [T^ 11501 1152j . Here we give the final result 
for X 
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with = 9* — igA^ and V"^^ = d'^^ + igA^ the covariant derivatives constructed with the background 
field where the derivative 9^ is acting on the function to its left. Furthermore, we have used the 
following short hand notation 
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Finally, we note that the above expressions look much simpler when one expresses the background 
field A in terms of the covariant gauge color charge density p. This is allowed since both the functional 
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measure and the weight function are gauge invariant. The classical solution is now simple and reads 
Aa = S'^'^aa where a°- = — T^-p". The non-linear evolution equation (IllbB . written in terms of the 
covariant fields a looks like 
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where we have defined 
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and X; 5" are related to cr via 
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with 



= V\x) = Pexp|i(7 J dz a{z ,a;)|. 
Combining all the expressions above, the final expression for 77"^ is given by 
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For completeness, we also give the final expression for the virtual corrections denoted v 
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Eq. ()123p with the coefficient functionals given by ()128I129|) is known as the JIMWLK equation. To 
conclude this section, we note that it is possible to rewrite the non-linear evolution equation in a more 
compact form, eliminating a in favor of Xi &s was done in Furthermore, it is possible to recast 
the JIMWLK equation as a Langevin equation which lends itself easily to numerical investigation, for 
example, on a lattice ^H]- For a review of the most recent progress in understanding the JIMWLK 
equation and its various forms, we refer the reader to |5Uj . 



2.3.2 The Balitsky-Kovchegov Equation 

Let us approach the problem of resumming quantum evolution corrections from a different side: instead 
of including small-x evolution corrections in the gluon distribution function, as was done in Sect. 12.3.11 
we will consider quantum corrections to the dipole- nucleus cross section from Sect. 12.2.11 In the quasi- 
classical limit, the forward amplitude of the dipole- nucleus scattering is given by Eq. (jMj) [HTj . obtained 
by resumming multiple rescatterings of Fig. ^2 Now let us study how the quantum corrections come 
into this multiple rescatterings picture. 

Similar to the BFKL evolution equation [771 and to JIMWLK equation, we are interested in quantum 
evolution in the leading longitudinal logarithmic approximation resumming the powers of 

1 
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Figure 20: Quantum corrections to dipole-nucleus scattering. 



with Y the rapidity variable. Again we will be working in the rest frame of the nucleus, but this time 
we choose to work in the light cone gauge of the projectile = if the dipole is moving in the light 
cone + direction. This gauge is equivalent to covariant gauge for the multiple rescatterings of Fig. ITTl 
therefore, all our discussion in Sect. 12.2.11 remains valid in this new gauge. 

We need to identify radiative corrections that bring in powers of Og Y. As we have seen in Sect. 
12.2. H instantaneous multiple rescatterings bring in only powers of not enhanced by factors of Y. 
Therefore, additional Coulomb gluon exchanges would not generate any logarithms of xbj bringing in 
only extra powers of a^. We are not interested in such corrections. Other possible corrections in the 
light cone gauge of the projectile dipole are modifications of the dipole wave function. The incoming 
dipole may have some gluons (and other quarks) present in its wave function. For instance, the dipole 
may emit a gluon before interacting with the target, and then the whole system of quark, anti-quark and 
the gluon would rescatter in the nucleus, as shown in the top diagram of Fig. |2I3 The dipole may emit 
two gluons which would then interact with the nucleus, along with the original qq pair, as shown in the 
bottom diagram of Fig. 1201 In principle we can have as may extra gluon emissions, as well as generation 
of extra qq pairs in the dipole's wave function. As we will shortly see, these gluonic fluctuations from 
Fig- EDI actually do bring in the factors of enhanced by powers of rapidity Y, i.e., they do generate 
leading logarithmic corrections. Fluctuations leading to formation of qq pairs actually enter at the 
subleading logarithmic level bringing in powers of Y |153| 11541 1155l| and are not important for the 
leading logarithmic approximation used in this Section. 

The gluons in the dipole wave function have coherence length of the order of (see Eq. (jH}) 



coh 



2k+ 



:i3ii 



if the incoming qq pair is moving in the light cone + direction. If A;+ is large enough, the coherence 
lengths of these gluons would be much larger than the nuclear radius, Icoh ^ -R, so that each gluon 
would coherently rescatter on the nucleons in the nucleus, just like the original dipole in Fig. ^2 This is 
indeed what is also shown in Fig. 1213 Note that gluons are emitted by the incoming dipole only before 
the multiple rescattering interaction (or after the interaction in the forward amplitude). Emissions 
during the interaction are suppressed by the powers of xbj (or, equivalently, by inverse powers of center 
of mass energy of the scattering system) [221 • This can be checked via an explicit calculation in the 
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Figure 21: In the large- A'c limit of the Mueller's dipole model [12112011211 the gluon cascade is represented 
as a color dipole cascade. Quark colors are labeled r, b, g for red, blue and green correspondingly. 

covariant Feynman perturbation theory or by performing the calculation in the light cone perturbation 
theory |135| I136j . In the latter case, the emission of a gluon is allowed and is equally probable at any 
point throughout the coherence length of the parent dipole 1'^^^^^ = 2p^/p'^ with p the momentum of 
the dipole and p^ being very large. The probability of emitting the gluon inside the nucleus is given 
by R/l'^coh ~ Vp^ ~ it is suppressed by powers of the center of mass energy s compared to 

emission outside the nucleus and can be neglected in the eikonal approximation considered here. 

Therefore, our goal is to resum the cascade of long-lived gluons, which the dipole in Fig. 1201 develops 
before interacting with the nucleus, and then convolute this cascade with the interaction amplitudes of 
the gluons with the nucleus. To resum the cascade we will use the dipole model developed by Mueller in 
POIEOIE]- Mueller's dipole model makes use of the 't Hooft's large- iVc limit taking Nc to be very 

large while keeping agNc constant. In the large-Nc limit only the planar diagrams contribute, with the 
gluon line replaced by a double line, corresponding to replacing the gluon by a quark-anti-quark pair in 
the adjoint representation. The diagrams of Fig. |201 translate into the planar large- iVc diagrams shown 
in Fig. 1^ The top diagram there represents emission of a single gluon in the original incoming dipole. 
After replacing the gluon by the double quark line, as shown on the top right of Fig. |^ the original 
dipole splits into two new dipoles, formed by the original quark line combined with the anti-quark line 
in the gluon (blue-anti-blue dipole in Fig. I?T|) and by the original anti-quark line combined with the 
quark line of the gluon (red-anti-red dipole in Fig. . Successive gluon emissions would only split the 
dipoles generated in the previous step into more dipoles, as depicted in the bottom diagram of Fig. |2I1 




A B 



Figure 22: One (real) step of dipole evolution. 

To understand how the dipole model |TOl 1201 12I| works let us first calculate one step of the evolution, 
shown in Fig. |221 There the original dipole consists of a quark-anti-quark pair, with the quark (top 
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line) having transverse coordinate and the anti-quark (bottom hne) having transverse coordinate x^. 
Suppose the initial dipole as a whole carries large light cone momentum p+ and the anti-quark carries 
the momentum kf. It will be convenient to define ratios of these momenta, such that the anti-quark 
would carry the fraction Zi = kf /p'^ of the dipole's light cone momentum and the quark would carry 
the fraction 1 — zi. Now suppose we emit a gluon with transverse coordinate X2 and the light cone 
momentum fraction Z2 = k2 /p'^ as shown in Fig. |221 As can be shown by an explicit calculation, 
such an emission generates a logarithm of energy only if ;22 <C 21, 1 — 21. In this approximation the 
transverse coordinates of the initial quark and anti-quark do not change: the emission is recoilless. 
A straightforward calculation of the diagram in Fig. [221 using the rules of the light cone perturbation 
theory |135j yields [12] 

_ 'igt" I ^ ^ ^ \ ^ (132) 

where we added the diagram where the gluon is emitted by the anti-quark to the graph in Fig. [221 In 
Eq. (jl32j) €3 is the polarization of the emitted gluon, t"" is the color matrix in fundamental representation 
and = x j — Xy Squaring the contribution of Eq. (|132|) . summing over all polarizations and colors of 
the emitted gluon, we obtain the following modification of the initial dipole's wave function 



^^'\xo,X„Zi) = ^^ d^X2 — ^^$^°H^o,2Ii,-2i) (133) 

TT J J Zin Z2 3^20 3^21 



where I and and are squares of the dipole light cone wave functions before and after 

the gluon emission correspondingly. In the spirit of the leading logarithmic approximation we cut off 
the 2;2-integral from above by min{2;i, 1 — zi} to indicate that Z2 -C z\^ \ — z\. however for the gluon 
evolution the exact choice of this cutoff does not matter. (A choice of the cutoff actually does matter 
for formulating the evolution of the non-singlet structure functions, the so-called reggeon evolution, in 
the language of the dipole model |157j .) Zin is some lower cutoff of the Z2 integration which is not 
important for our purposes here. 
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Figure 23: One (virtual) step of dipole evolution. 

In addition to the real emission term of Fig. [221 we need to add virtual corrections, one of which is 
demonstrated in Fig. [2S1 Other virtual diagrams involve the gluon coupling to both the quark and the 
anti-quark lines, as well as only to quark line. The virtual corrections should cancel the real emission 
graph of Fig. [221 in the wave function squared if there is no interaction with the target (for a more 
detailed analysis of real- virtual cancellations see [22]) • Therefore, to calculate the corrections we just 
have to perform the X2-integration in Eq. ()133p inserting the overall minus sign for the virtual term. 
This yields the following modification of the dipole wave function 

*<»(^,^.,,.) = -i^l„f£2l) f"""""';^*'»'fe,^..0. (134) 
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Figure 24: DIS as a dipole cascade interacting with the target (see text). 

where p is some ultraviolet cutoff needed to regulate the X2-integration in Eq. ()133|) : it will cancel out 
for physical quantities. 

Now resummation of the gluon cascade of Fig. 1201 becomes more tractable [23- First of all, in the 
\aige-Nc limit the gluon cascade translates into a dipole cascade of Fig. |^ As we have seen in arriving 
at Eqs. ()133p and ()134|) . in the leading logarithmic approximation the gluon emissions do not change the 
transverse coordinates of the quark and anti-quark lines off which the gluons were emitted. Therefore, 
the color dipoles have the same transverse coordinates throughout the whole process: once they are 
created the transverse coordinates do not change. Now resummation of the dipole cascade reduces to 
the set of diagrams represented in Fig. |^ which is a generalization of Fig. ^2 for the case of quantum 
evolution corrections. (The cascade in Fig. |^ is very similar to the one discussed in Sect. 12.1.11 shown 
in Fig. In DIS the incoming virtual photon splits into a qq pair, which we will refer to as the original 
parent dipole. The parent dipole may emit a soft gluon, like in Fig. |221 splitting itself into a pair of 
dipoles. (In Fig. ED the sum of gluon emission from the quark and from the anti-quark lines is denoted 
by the gluon (double quark) line disconnected from the parent quark and anti-quark lines.) The dipole 
may also emit and re-absorb a gluon, generating virtual corrections of Fig. 1221 The produced dipoles 
may also evolve, generating more and more dipoles, as shown in Fig. (211 

In the end the evolved system of dipoles interacts with the nucleus. In the large- A'^c limit each dipole 
does not interact with other dipoles during the evolution which generates all the dipoles. For a large 
nucleus the dipole-nucleus interaction was calculated in Sect. 12.2.11 vielding the answer in Eq. (jHBj) . 
That result resums powers of a^A^/^^ as shown in Eq. (j64j) . Analyzing the diagrams for the interaction 
of several dipoles with the nucleus we see that interaction of, say, two dipoles with a single nucleon 
is suppressed by a power of A^^^ and is therefore subleading and can be neglected. Interaction of two 
dipoles with two nucleons in the large- A'c limit is dominated by the diagram where each of the dipoles 
interacts with only one of the nucleons (if we require both dipoles to interact). In general one can argue 
that, in the large- A^c limit and at the leading order in A (or, equivalently, resumming powers of a1 A^/^), 
the interaction of any number of dipoles with the nucleus is dominated by independent interactions of 
each of the dipoles with a different set of nucleons in the nucleus through multiple rescatterings shown 
in Fig. ^2 This is depicted in Fig. (211 there each of the dipoles present in the dipole wave function 
by the time it hits the nucleus may interact with different nucleons in the nucleus by exchanging pairs 
of gluons. (Only interactions of some of the dipoles are shown.) Therefore, the dipoles are completely 
non-interacting with each other: they do not exchange gluons in the process of evolution, since those 
corrections would be suppressed by powers of Nc, and they interact with different nucleons in the 
nucleus, which is correct at the leading order in A j24j . 

Summation of the dipole cascade of Fig. now becomes straightforward and is illustrated in Fig. 031 
|24j . There the dipole cascade and its interaction with the target are denoted by a shaded oval. In one 
step of the evolution in energy (or rapidity) a soft gluon is emitted in the dipole. If the gluon is real 
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Figure 25: Diagrammatic representation of the nonlinear evolution equation ()135p . 



(see Fig. I22|) . than the original dipole would be split into two dipoles, as shown in Fig. |25| Either one 
of these dipoles can interact with the nucleus with the other one not interacting, which is shown by the 
first term on the right hand side of Fig. 123 with the factor of 2 accounting for the fact that there are 
two dipoles in the wave function now. Alternatively, both dipoles may interact simultaneously, which 
is shown by the second term on the right hand side of Fig. 123 This term comes in with the minus sign 
due to the minus sign in our definition of the amplitude N in terms of the S'-matrix, S = 1 — N. The 
emitted gluon in one step of evolution may be a virtual correction of Fig. 1221 which is not shown in 
Fig. 123 in that case the original dipole would not split into two, it would remain the same and would 
interact with the target. 

Combining the terms from Fig. 123 with the virtual correction term of Fig. 1221 and using Eqs. ()133|) 
and (I134D we write 
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N{x^2,b +^X2i,y) + N{xi2, b+'^X2o,Y)- N{xo^,b, Y) 



(135) 



Note that since the impact parameter of the original dipole is 6 = (xq + aii)/2, the impact parameters 
of the produced dipoles X02 and X12 are different from b on the right hand side of Eq. (jl35|) . 

Equation ()135|1 is a nonlinear evolution equation, a solution of which gives us the forward amplitude 
of the dipole-nucleus scattering at high energy. The initial condition for Eq. ()135|) at some initial 
rapidity which we denote F = is given by Eq. ()66|) 24 . Therefore, Eq. p35p resums all powers of 
the multiple rescattering parameter A^^^ along with the leading logarithms of energy in the large- iV^ 
limit given by powers of Nc Y. Eq. p35j) was originally derived by Balitsky in the framework of the 
effective theory of high energy interactions in ||2^ and, independently, by one of the authors in [21] 
using the formalism of the dipole model JHl 1201 1^- It is commonly referred to as the BK equation. 
It was also re-derived by Braun in [77| using the expression for the triple pomeron vertex from |141j in 
resummation of fan diagrams in Fig. 13 

Eq. p35j) has a structure similar to the GLR equation [H I139j and appears to correspond to 
summation of the fan diagrams of Fig. |3 If one neglects the impact parameter dependence of the dipole 
amplitude N and neglects the dependence of n on the angle of Xgi, such that N{x[^i,b, Y) ^ N{xqi, Y) 
one can perform a Fourier transform 
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to obtain [211 
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where xif^^^) was defined in Eq. (jHUj) . One can show that Eq. p37|) is equivalent to Eq. ()46|) if one 
identifies 

<P{x,k^) = -^N{k,Y = \nl/xBj), (138) 
which, combined with Eq. ()13(jj) . leads to 

(f)(x,k^) = ^^^^ [ d^xe-'^-^^N(x±,Y = \nl/xBi). (139) 

For from Eq. (jfjfjj) the unintegrated gluon distribution function from Eq. p39|l looks similar to the 
gluon distribution from Eq. ()102|) . However, the important difference is that the Glauber- Mueller dipole 
amplitude Nq which enters Eq. fll02|) is the amplitude for a gluon (adjoint) dipole, while N in Eq. ()139p 
is the amplitude for a quark (fundamental) dipole. Even this formal difference notwithstanding, it is 
not at all clear whether Eq. ()139p . or even Eq. ()102|) . would still be valid for the dipole amplitude 
from Eq. (jl35|) with all the evolution effects included: after all, Eq. (jlU2p was derived only in the quasi- 
classical limit. We will return to this question in Sect. 13. II where we will discuss particle production. 



2.4 Solving the Evolution Equations 

After the JIMWLK and the BK evolution equation were written down there has been a number of 
analytical O EH EHl IH EHj and numerical [ZZl EHl IZOl EDI EH EH] attempts to solve them. While an 
exact analytical solution valid both inside and outside of the saturation region still does not exist, there 
is a number of good approximations for various kinematic regions. Here we are not going to give an 
extensive review of the existing approaches referring the interested reader to the existing review articles 
jUl HSl HO] ■ Instead we will outline features of the solution of the evolution equations which will be 
most important for our discussion of particle production below. 



2.4.1 Linear Evolution 

When the dipole- nucleus interactions are weak and the forward scattering amplitude is small A^ <C 1, 
we can neglect the quadratic term in Eq. (jl35|) and write [121 120] 
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(140) 

This equation is equivalent to the BFKL equation [27], as has been shown in the framework of the 
dipole model in |158j . Defining the the eigenfunctions of the Casimir operators of conformal algebra 

£-(p„,pO = «S_ IJ^) (141) 

with (in general) complex coordinates pi, one can show that they are also eigenfunctions of the dipole 
kernel of Eq. [HHl HHO] 



/ d^X2 [E^^'ixo, X2) + E^'''ix2, xi) - E^^'^ixo, a;i)] = 4 TT xin, u) E"'-(xo, x,) (142) 

-^02 -^12 



with x(n, I/) given by Eq. ([HU]) . Eq. ()142|) demonstrates that the BFKL equation [27] and Eq. ()140p 
have the same eigenfunctions and eigenvalues, and are, therefore, equivalent. 
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Let us write down explicitly the solution of Eq. ()140|) for the case of a large nucleus, where we can 
neglect the 6-dependence in A^(xoi, Y) along with its dependence on the azimuthal angle of Xqi- Using 
again the Fourier transform of Eq. (jlHfjj) we rewrite Eq. ()14()j) as (see Eq. (jTHZj)) 



dN{k,Y) 2asN^ 
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the solution of which is given by 
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where C{k) is some unknown function, which could be determined by using the initial conditions (jUUI) 
in Eq. ()136|) . To determine the high energy asymptotics of the solution ()144j) we write C{k) as a Mellin 
transform 

dv 
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Here Qso is the saturation scale in the quasi-classical approximation given by Eq. (j65p . From here on 
we will use the subscript to distinguish it from the full energy-dependent saturation scale which would 
result from our analysis of the JIMWLK and BK evolution equations. 
Substituting Eq. ()145|) in Eq. ()144|) we obtain 
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where we have defined 
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Substituting the result of Eq. ()14(i|l into Eq. ()13(i|) and integrating over k yields 

/oo _ 
du e2"»>^(°'^)^(xTQso)'+'"^a 
-oo 

where 
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Eq. (jl48|) gives us the solution of the dipole version of the BFKL equation (jl4(jp with the coefficients 
Cu to be fixed from the initial conditions of Eq. ()66|) . Now let us rewrite Eq. ()148|) as 



/oo 
du exp[2asX{0,'^)Y + 2iu In(xTQso)] ( 150) 

-oo 

The integral in Eq. ()150|) can be approximated by the saddle point method in the following two important 



cases. 



XtQso ^ 1 In this case the term 2iu \n{xTQso) in the exponent of Eq. ()150|) is small compared 
to as Y and does not significantly infiuence the position of the saddle point. The saddle point is 
determined mostly by the function x(0, z/), which is shown in Fig. [JHlby thick line as a function of 
—i v: there one can clearly see a saddle point near i/ = 0. We use the expansion of x(0, v) around 
p = Q (see Eq. ((23)) to determine the saddle point 
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Figure 26: The eigenvalue of the BFKL kernel x(0, i^) plotted as a function of —iu (thick line). The 
thin straight line is due to the linear term 2ii> ln{xTQso) in the exponent of Eq. p50|) . The thin curve 
is a sum of the thick line and the thin straight line and it represents the complete expression in the 
exponent of Eq. ()15U|) . 



where the subscript LLA stands for leading logarithmic approximation. Performing the saddle 
point integration we write (cf. Eq. pUj)) 



N{xr,Y) = ^tQsoC^i^J^^^^^ exp 
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with the BFKL pomeron intercept ap — 1 given by Eq. (jHTjl . 



xtQso <^ 1 This case is also known as the double logarithmic approximation (DLA), since here 
transverse logarithms like ln{xTQso) become important, leading to a new resummation parameter 
OsY ln(xT<5so)- Transverse logarithms are of course resummed by the Dokshitzer, Gribov, Li- 
patov, Altarelli, Parisi (DGLAP) evolution equation |142j . The DLA region is where the results 
of DGLAP and BFKL are identical |134j . since there both equations resum powers of the same 
parameter agY Iii^xtQso)- 

As shown in Fig. EHlby the thin curve, when ln{l/xTQso) becomes large, it shifts the position of 
saddle point towards u = —i/2 |il61j . In that region we approximate x(0, i') by 



which leads to the saddle point at 
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Using Eq. p53j] in Eq. p5()j) and performing the integration over v in the saddle point approxi- 
mation yields 

Ar(xT,F) = (xTQ.o)'aB,,^(2«.r)i/Mn-3/^(^^-^^ exp 

(155) 

The coefficient C^*^^^ in Eq. ()155|) is in fact important [62 , since it may modify the prefactor: if 
we require, for instance, that at F = Eq. ()150p reduces to the amplitude given by a two-gluon 
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exchange from Eq. or, equivalently, by the first term in expansion of Eq. (jUBj) . than we would 
get 



= -. (156) 



Using Eq. (fT^ in Eq. ^TT^ yields 



Ni.r.Y) ^ (..Q.o)^^(2Z..y)-3/MnV^(-L-) exp 
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(157) 



The general analytic linear solution found above in Eq. p5()j) allows one to estimate the boundary of 
the saturation region by determining the limit of applicability of the linear equation (jl4(Jj) [T]l I73| 17^. 
Linear evolution is applicable only if <^ 1, and, therefore, breaks down when ~ 1. To estimate 
when that happens, let us rewrite Eq. ()150p as 



/oo 
-oo 



158) 



Following [71] we argue that becomes of order one when the value of the exponent in Eq. ()158|) at 
the saddle point vanishes. This translates into the following two equations 

2a,x'{0,^o)Y + 2ilniQ,o/Qs) = (159) 

(saddle point condition) and 

2asXiO,iyo)Y + {2tiyo + l) HQso/Qs) = (160) 

(vanishing of the power of the exponent at the saddle point). We have substituted ~ l/Qs in Eqs. 
()159p and ()16Up . which corresponds to deGningthe saturation scale Qs by requiring that N{l/Qs, Y) ^ 1. 
One can see that this definition is consistent with the quasi-classical approach of Sect. 12.21 if we put 
y = in Eqs. ()159|) and ()160|1 . which corresponds to removing the small-x evolution, we obtain 
Qs = Qso as the solution, meaning that the full saturation scale Qs maps onto the quasi-classical 
saturation scale Qso in the small rapidity limit. 
Solving Eqs. and (fMUl yields [Zi 

Qs{Y) = g,oexp|2a,;^^^r| (161) 

with 

= (162) 

giving 

z/o ^ -i 0.1276. (163) 

Using this value of uq in Eq. p6H) we get 

QsiY) ^ g.oe^-^^^=^. (164) 
A more detailed determination of the saturation scale specifying the prefactor in Eq. p6ip can be found 

in nans]. 

Eq. ()161|) . or, equivalently, Eq. ()164|) . give us the dependence of the saturation scale Qs on rapidity or 
energy. We see that as rapidity Y increases, so does the saturation scale. Eq. ()161|) generalizes Eq. 
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Figure 27: Our understanding of high energy QCD interactions plotted in the plane of rapidity Y = 
Inl/xBj and InQ^ (see text). 

for Qso by including the effects of small- a; evolution, which bring in rapidity dependence. Combining 
the two equations one writes 



Eq. ()165p demonstrates that for a large nucleus and at high energies the saturation scale is proportional 
to 



with c some constant, c = 4.88 in Eq. ()164j) . Therefore, Qg is an increasing function of both the atomic 
number A and of center of mass energy s. It only gets larger as we go to higher energies, making 
the coupling constant as{Qs) smaller. Thus the saturation approach is better justified as we go to 
progressively higher energies. 

Our understanding of high energy scattering at this point in the article is summarized in Fig. 1771 
There we depict the plane of two variables essential to the high energy scattering: the typical transverse 
momentum scale and the rapidity Y = Inl/xsj- Right away the region of Q < Aqcd should be 
excluded from our analysis here, since there the coupling is large and not much can be understood 
using the perturbative methods we describe here. At Q ^ ^qcd we may distinguish two important 
regions. One region is given hj Q > Qs{Y), where the amplitude is small and the linear BFKL 
evolution equation p40|) 27 applies, along with the DGLAP evolution |142j . as shown in Fig. |2ZI In 
the region with Q < Qs{Y) the amplitude A^ becomes of order 1 and the non-linear effects of the BK 
and JIMWLK equations [Zj-[II], [2S1I21] become important: this region is called the saturation region 
and is given by the shaded area in Fig. 123 Saturation region is also called Color Glass Condensate. It 
is important to point out that all this non-linear dynamics takes place at Qs^Q ^ ^qcd, i-e., in the 
perturbative region where the strong coupling constant is small and our calculations are justified. 

2.4.2 Geometric Scaling 

Let us now analyze the behavior of the solution of Eq. ()135|) deep inside the saturation region, where 
the non-linear effects are very important. Deep inside the saturation region, when the dipole size xt 




(165) 



(166) 
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becomes large, xt ^ (but still xt <^ I/Aqcd), the quasi-classical Glauber-Mueller amplitude 

from Eq. approaches 1. Analyzing Eq. (jl35|) we easily see that = 1 is also a stationary solution 
of that equation. Therefore, we conclude that 

Ar(x,6,y) = 1, XT->llQs{Y), (167) 

corresponding to the black disk limit of Eq. ()68j) |146j . Now let us determine the asymptotic approach 
to the black disk limit of Eq. p67|) following [721 ESI- To do that, let us write 

N{x,b,Y) = I- S{x,b,Y) (168) 

where S is the S'-matrix of the dipole-nucleus collision, which is small when the system approaches the 
black disk limit (A^ approaches 1). Substituting Eq. p68p in Eq. ()135|) and keeping only terms linear 
in S yields [HI IS Hi 

dS{xQi,h,Y) _ a^Cp f ,2 ^< 



dY TT 



2 



2 

f rf2^2^^^(xoi,6,y), (169) 
J X2Q X21 



where the integral over dipole sizes goes over Xo2,a;i2 > l/Qs{Y). Hence we have to replace the 
ultraviolet (UV) cutoff from Eq. (fil^ with l/Qs{Y) obtaining 

= InMsiY)] S{x,„b,Y). (170) 

Defining the scaling variable 



TT 



In 



4mY) (171) 



with (cf. Eq. ^EB) 
we rewrite Eq. ()170p as 



_ d In [xl,Q',{Y)] 
^^^=dY = dY 

f ^ -> 



Solution of Eq. ()173|) can be straightforwardly written as [721 113 

S{^) = Soe-^'/' = ^oe-"^'Ni«'(^)]/= (I74) 
with 5*0 < 1 a constant. Corresponding dipole amphtude A^ is given by 

N{^) = 1 - Soe-^'/' = 1 - 5oe-'°'["i«'(^)]/^ xt > l/Qs{Y). (175) 

Eq. ()175p is knows as Levin- Tuchin formula |72j . 

Note that the S'-matrix and the amplitude A^ of the dipole-nucleus scattering given by Eqs. (|174j) 
and (|175|) are functions of a single variable ^, or, more precisely, of the combination x±Qs{Y). This 
phenomenon is known as geometric scaling. The original argument for geometric scaling has been given 
by McLerran and Venugopalan in fS^, where it was suggested that small- a; nuclear or hadronic wave 
functions are described by a momentum scale Qs being the only scale in the problem and, therefore, all 
transverse coordinate (or momentum) dependent physical observables should depend on the combination 
x±Qs{Y) (or kx/Qs)- Geometric scaling has been demonstrated in an analysis of the HERA DIS data by 
Stasto, Golec-Biernat and Kwiecinski in [71], presenting one of the strongest arguments for observation 
of saturation phenomena at HERA. These results are shown here in Fig.|2Hlfrom [7T|, where the authors 
of [TT] plot combined HERA data on the total DIS 7*p cross section a^ot for xsj < 0.01 as a function 
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Figure 28: HERA data on the total DIS 7*p cross section plotted in jJIj as a function of the scaling 
variable r = jQ^^XBi^- 
oi the scaling variable r = Q"^ /Ql{xBj)- One can see that, amazingly enough, all the data falls on the 
same curve, indicating that a'tot is a function of a single variable Q'^/QKxBjY This gives us the best 
to date experimental proof of geometric scaling. (For a similar analysis of DIS data on nuclear targets 
see and the first reference of [HT].) 

The fact that geometric scaling is a property of the solution of the BK equation has been later 
demonstrated in I73j . In the next Section we will also show that such scaling is also valid outside of 
the Q < QsiX) saturation region. 

2.4.3 Extended Geometric Scaling 

In j73] lancu, Itakura and McLerran noticed that the geometric scaling of [72], which was originally 
derived as a solution of the BK equation inside the saturation region, also applies in an area outside 
of that region. To see this let us first remember that outside of saturation region {xt < l/Qs{Y)) the 
nonlinear BK equation reduces to the linear BFKL evolution (jl4()|l . Therefore, the dipole- nucleus 
amplitude is still given by Eq. ()158|) . which we rewrite here as 

/CO 
rfz/e^(^)a (176) 
-oo 

with 

P(z/) = 2a,xiO,u)Y + i2tu + l)\n{xTQso). (177) 
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In order to perform the integration in Eq. ()176|) around the saddle point from Eq. ()162p we expand 
P(z/) around z/q 

P(z/) = P(z/o) + P'ii^o) + \ P"{i^o) {u-uof + .... (178) 

Following [73] we first assume that the integral in Eq. ()176|) is dominated by its value at the saddle 
point z/q. Then noticing that 

P{uo) = (2ZZ/0 + 1) \n[xTQs{Y)] = Q + ^^o) ^ (179) 

we obtain 

N{XT,Y) ^ aoe^^'^") = a„ MsiY)]'^''"' = aoe(^^''"")^. (180) 



Eq. ()180|) obviously exhibits geometric scaling, since the amplitude in it is a function of ^ only. 
Therefore geometric scaling also works outside of the saturation region for xt < l/QsiY) {Q > QsiY)) 
[73j . This effect is called extended geometric scaling. 

To determine the applicability region of Eq. p81|) we use the formula ()178p . truncate the series in 
it at the quadratic order, substitute the result in Eq. p76|) and integrate over u obtaining 



iV(.T. Y) = C„, exp [PM - , (181) 

Since 

P'{uo) = 2z In [xtQs{Y)] = (182) 

and 

P"{uo) = 2a,x"{0,uo)Y (183) 

the second (diffusion) term in the exponent of Eq. p81|) along with the prefactor in Eq. ()181|) violate 
geometric scaling. These scaling violation corrections are small and can be neglected as long as 

M « Pi.). (184) 

Using Eqs. (fT79|l . (fT82|l and (fT83|l in Eq. (fT8^ yields 

1 



Vsl^ ) 

Replacing xt by 1/Q this translates into 

Q < Q^(y)e~^»>^"(0'-«)(2^-o+i)y _ Q^(y) e^o-^"^^. (186) 



However, the condition in Eq. ()18(ij) . while being necessary for extended geometric scaling to apply, is 
not sufficient. What we learn from Eq. ()18H1 is that geometric scaling is violated when the regions of 
ly further away from z/q start to contribute in the integral of Eq. ()176p . This is also demonstrated in 
the DLA approximation: clearly Eq. p57|) can not be written as a function of a single variable ^ and 
thus violates geometric scaling. We may, therefore, define the border of the extended geometric scaling 
region by the transition point between the LLA saddle point of Eq. (jl5ip and the DLA saddle point of 
Eq. (jl54j) |i73j . The point of closest approach of the two saddle points in Eqs. (jl51|) and (jl54|) is at 

In — ^ ^ 3.28 asY (187) 
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Figure 29: The summary of our knowledge of high energy QCD interactions plotted again in the plane 
of rapidity Y = Inl/x^j and InQ^. 

such that the border of extended geometric scaling region is defined by the scale 

kgeom = Qsoe^'^^°°^ (188) 
which, using Eq. (|164|) can be rewritten as [TSl 102] 

fcgeom = QsiY) (^P) ■ (189) 

Therefore, extended geometric scaling is valid up to 

Q < fcgeom, (190) 

which is a more restrictive condition than Eq. (jl86|) . The exact value of fcgeom niay still be slightly 
different from the one given by Eq. ()189|) [TSj [731 UM- what is important is that at large Y this scale 
is much larger than the saturation scale, Qs(Y) <^ kgeom, allowing for a parametrically wide region of 
extended geometric scaling. 

The summary of our current knowledge of high energy interactions is shown in Fig. 1221 One can see 
there that, on top of the saturation region Q < Qs{Y) of Fig. (271 one now has the extended geometric 
scaling region Qs{Y) < Q < kgeom, where the linear BFKL evolution still applies with its solution 
having the property of geometric scaling due to the presence of saturation region. 

Indeed the property of extended geometric scaling was derived here for the case of DIS where a 
small projectile {qq dipole) scatters on a nucleus. Saturation effects were only present in the nuclear 
wave function. At extremely high energies gluon saturation may take place even in the projectile's 
wave function: to take this into account pomeron loop resummation needs to be performed jl24j - [T^H] . 
While such resummation is still an open problem, it has been argued in |163| I164j that such pomeron 
loop effects along with energy conservation constraints may potentially lead to a violation of geometric 
scaling at very high energies. 
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The interested reader, who would hke to learn more about how to obtain the above analytical 
solutions in a more complete way keeping all the prefactors and on how to include the running coupling 
effects into the problem is referred to [73 EU [7^] for further reading. An elegant solution of the 
BK equation, reducing it in the geometric scaling region to Kolmogorov-Petrovsky-Piskunov (KPP) 
equation, which has traveling wave solutions, was found in 

2.4.4 Numerical Solutions 

As we have mentioned before, both the BK and the JIMWLK evolution equations have been studied 
numerically in [ZZl UHl IZHl IHOl IMl QH|- Here we are not going to give a comprehensive review of these 
solutions, but will just show the results of one of them which summarizes all the important qualitative 
features of the nonlinear evolution and shows how it resolves the problems of the BFKL equation stated 
in the end of Sect. 12.1.11 

As we have already seen in Sect. 12.4.21 the solution of BK evolution equation approaches = 1 at 
very high energy. Using Eq. (j67p we see that this behavior corresponds to the black disk limit for total 
cross sections given by Eq. (jUHj) . The total cross section in Eq. (jUHj) is constant and, therefore, does 
not violate the Froissart bound of Eq. ()40|) . Therefore the nonlinear evolution appears to resolve the 
unitarity problem of the BFKL equation, posed as question (i) in Sect. EHIU ES EZl El 113 IHIZl EHl 
17^ ISUl IHH EH] • (More precisely, unitarity is restored at a given impact parameter b: integration over b 
in Eq. (fHTjl may still give a cross section growing with energy due to the diffusion of the boundary of the 
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Figure 30: The contour plot of the numerical solutions of the BFKL and BK evolution equations in 
momentum space [kx N{k, Y)) as a function of transverse momentum k and rapidity Y = In l/x^j from 
|7H|. 
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black disk which increases its radius R and the area vri?^ and, hence, the total cross section atot = 2 tt i?^ 
(see |1651 1166] and the second reference in [75j).) 

To answer the question (ii) in Sect. 12.1.11 regarding the diffusion into the infrared shown by the 
Bartels cone of Fig. IHl we will present one result from the numerical solution of the BK equation done 
in [ZHj- In Fig. EOl the authors of [ZHj plot the lines of constant value for the numerical solution of 
the BFKL and BK equations in momentum space. Namely they make a contour plot of kx N{k,Y) 
from Eq. p36|) scaled down by the maximum value reached by that function in the phase space region 
considered |j7Sj as a function of transverse momentum k and rapidity Y = Inl/xsj. One can see that 
the solution of the BFKL equation (dashed lines in Fig. IHUj) spreads out as the rapidity increases [xej 
decreases). This is the diffusion discussed in Sect. 12.1.11 [26]. which is dangerous because it generated 
non-perturbative low-fc-p gluons, for which our small coupling treatment would not apply. Fig. IHUl shows 
that the nonlinear BK evolution (shown by solid lines in Fig. IHUj) avoids this problem. The effect of non- 
linear term in Eq. ()135|) is to drive the constant value lines of the solution towards higher momenta, 
which is consistent with increase of the saturation scale in Eq. fll61j) . and to eliminate the diffusion 
spread of the solution: as one can see from Fig. 1201 the width of the fcj-distribution of the BK solution 
is roughly independent of rapidity. This solves the IR diffusion problem posed as question (ii) in Sect. 
12.1.11 Similar results were obtained by other numerical simulations of the BK [771 IZ3 HOI IM] and the 
JIMWLK [IH] evolution equations. 



3 Particle Production in pA Collisions 

In this Section we apply the formalism of saturation/Color Glass Condensate physics to the problem 
of particle production in proton-nucleus (pA) collisions, which is directly relevant to deuteron-gold 
collision experiments at RHIC advertised in our paper's title. The Section is structured in the following 
way: first, in Sect. 13.11 we calculate gluon production in the quasi-classical approximation (McLerran- 
Venugopalan model) of Sect. 12.21 We will then continue in Sect. 13.21 bv including the effects of BK and 
JIMWLK quantum evolution from Sect. 12.31 into the obtained expressions for gluon production cross 
section. We then calculate valence quark production in Sect. 13.31 in the approach where inclusion of 
quantum evolution is straightforward. We move on to electromagnetic probes in Sect. 13.41 bv deriving 
prompt photon and dilepton production cross sections. In Sect. 13.51 we study two-particle correlations 
by calculating two-gluon, gluon-valence quark and qq pair production cross sections. 

3.1 Gluon Production in the Classical Approximation 
3.1.1 Gluon Production Cross Section 

Here we are interested in calculating the inclusive single gluon production cross section in pA collisions. 
We assume that, in the center of mass frame, the gluons light cone momentum components are much 
smaller than that of the proton and the nucleus. If in the center of mass frame the proton has a large 
p~ component of its light cone momentum, we will be interested in produced gluons with k^ ^ p^. 
Conversely, if a nucleon in the nucleus has a large p'~^ component of its momentum, then k'^ <^ p'^. 

Just like for the gluon field of a single nucleus in Sect. \'2.'2.'2[ the problem of gluon production 
in pA collisions in the quasi-classical approximation can be formulated as the problem of finding the 
classical gluon field from the Yang- Mills equations [Eq. (f72|) ] with the source current now given by both 
the nucleus and the proton |2H1 1201 1201 1211 122] • The outgoing gluon line of the resulting gluon field 
can be truncated and the gluon can be put on the mass shell, giving the production amplitude needed 
for constructing production cross section. The classical gluon field of the proton and the nucleus, or, 
equivalently, the resulting production cross section, would resum powers of the parameter A^^^ (see 
Eq. just like the classical gluon distribution from Eq. fll02|) or the dipole-nucleus forward scattering 
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Figure 31: Diagrams contributing to gluon production in pA collision in = (proton) light cone 
gauge in the quasi-classical approximation. 



amplitude in Eq. (jUnj) . This program has been carried out in jSH] and the classical gluon field produced 
in pA has been found there. However, here we use a different strategy: instead of calculating the 
classical gluon field, we calculate the production cross section directly by summing diagrams, following 
El (see also |H2l Ell II321 ) • 

We begin by choosing the gauge: it turns out to be much simpler to work in the light cone gauge of 
the projectile (proton). Let us agree that, in the center of mass frame, the proton is moving in the light 
cone — direction: that means we will be working in A~ = light cone gauge. Again we will use the 
light cone perturbation theory (LCPT) |135| 1136] . The diagrams contributing to gluon production in 
A~ = light cone gauge in the LCPT framework are shown in Fig. |^ where one should also add (but 
we do not show) the complex conjugate of the graph in Fig. 131b. which is obtained from the graph in 
Fig. 131b by mirror-refiecting it with respect to a vertical axis. 

The physical picture of the gluon production is the following: the incoming proton may already 
have a gluon in its light cone wave function before the collision with the nucleus and the system of the 
proton and gluon would multiply rescatter on the nucleons in the nucleus. Alternatively the proton 
can emit the gluon after the multiple rescatterings in the nucleus. Just like in derivation of the BK 
equation in Sect. 12.3.21 the diagrams where the gluon is emitted during the proton's passing through 
the nucleus are suppressed by powers of its large light cone momentum p~ , i.e., by powers of the center 
of mass energy of the system (see the discussion following Eq. (I131|) ). and can be neglected in the 
spirit of the eikonal approximation. The graph in Fig. I^Tt \ corresponds to the square of the amplitude 
given by the case when the gluon is present in proton's wave function before the collision. As one can 
see in Fig. I^TIA- in the gluon-proton system only the gluon interacts with the nucleons in the nucleus: 
interactions of the proton cancel by real- virtual cancellations • Moving a gluon exchanged between 
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the nucleus and the proton across the cut does not change the momentum of the produced gluon in 
Fig. I^TR. but does change the sign of the whole term, causing the cancellation. That is why we have to 
consider only the interactions with the gluon in Fig. I^TR. The diagram in Fig. 131b gives the interference 
term between the amplitude from Fig. IHTR with the gluon present in the proton's wave function before 
the interaction and the amplitude in which the gluon is emitted by the proton after the collision. Of 
course a diagram complex conjugate to Fig. 131b should also be included. Real-virtual cancellations 
of Fig. I3TK do not happen in Fig. l3Tt3. Moving an exchanged (Coulomb) gluon across the cut would 
force us to move it across the gluon emission vertex for the produced gluon on the right hand side, thus 
changing the momentum of the produced gluon. Thus the interactions of the nucleons with both the 
proton and a gluon have to be included in Fig. 131b. On the right hand side of the diagram in Fig. 131b 
only the interactions with the proton are possible. It can also be shown that the square of the diagram 
with late gluon emission (after the proton passes through the system, as shown on the right hand side 
of Fig. 131b) does not have any interactions in it and can be neglected. (The gluon exchanges between 
the proton and the nucleus cancel by real- virtual cancellations 32 .) 

Note that in the quasi-classical approximation depicted in Fig.|^the interaction is modeled by single 
and double gluon exchanges. Similar to Fig. ^2 and Sect. 12.21 we have to impose the limit of no more 
than two gluons per nucleon [5]. If a particular nucleon exchanges a gluon with the incoming proton 
and/or gluon in the amplitude then it has to exchange a gluon in the complex conjugate amplitude to 
"remain" color neutral. Alternatively the nucleon can exchange two gluons in the amplitude (complex 
conjugate amplitude), but then it would not be able to interact in the complex conjugate amplitude 
(amplitude). This is done in the spirit of the quasi-classical approximation resumming all powers of 

A^^^, as was discussed in Sect. 12.21 

The calculation is easier to perform in transverse coordinate space. We assume that the outgoing 
gluon has a transverse coordinate x to the left of the cut in Fig. ED and a transverse coordinate y to the 
right of the cut. To calculate the production cross section we will first calculate the amplitude M{x) 
and its conjugate M*{y) in the transverse coordinate space, Fourier-transform them into transverse 
momentum space and take their product, which would be the square of the amplitude in the momentum 
space M{k), giving the production cross section 

|M(^)P = ——r / d^x^ye-'^-^^-y) M{x)M*{y). (191) 



(Pk(Phdy 2(27r)3' 2{2Tif . 

To obtain the answer for the gluon production cross section in pA in the quasi-classical approximation 
we have to convolute the wave function of the proton having a soft gluon in it with the interaction 
amplitude of a gluon scattering in the nucleus. Below we will model the proton by a single valence 
quark: generalization back to the real proton case will be manifest. The soft gluon wave function of a 
single quark is given by the first term in the parenthesis of Eq. fjl32j) corresponding to the diagram in 
Fig. 1221 where we should disregards the anti-quark (lower) line. The interaction amplitude between the 
gluon and the nucleus squared can be calculated directly following [321 1116] . However, it is easier to 
use crossing symmetry |l(i7j and "reflect" the gluon having transverse coordinate y from the complex 
conjugate amplitude with respect to the cut placing it into the amplitude. One would then see that 
the interaction amplitude squared would reduce to the forward amplitude of a gluon dipole formed by 
gluons at x_ and y scattering on the nucleus. In the quasi-classical approximation such amplitude is 
given by Eq. ()103p . Combining the first term in Eq. ()132|) squared with the amplitude from Eq. ()lU3p 
yields the following contribution to production cross section coming from the diagram in Fig. I^Tt \ 

4^ = - I d%d'xd^-\,^^e-'^-'^-y) (e-(^-# MiMA)/4 _ 1) (192) 
d^k dy Tc J [ZTT)'^ n ary^ ^ ' 

with the gluon saturation scale Qso given by Eq. ()10H) . where now the subscript denotes the satu- 
ration scale of the quasi-classical McLerran-Venugopalan model without quantum evolution/rapidity 
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dependence. In Eq. fll92j) we have also assumed that the proton's transverse coordinate is 0. 

The diagram in Fig. 131b is calculated in a similar way. The incoming proton wave function there 
is the same as in Fig. I^Tl\. To resum the interactions with the nucleus we can either follow direct 
calculations of P^lll6j or reflect the quark line from the complex conjugate amplitude with respect to 
the cut into the amplitude, making it an anti-quark in the amplitude. The reflected quark line would 
then be on top of the quark line in the amplitude, since they both have the same transverse coordinate 
0. However the quark and the anti-quark do not cancel each other: in fact, since they emit a gluon, the 
pair of quarks should be in the color octet state. We again have an interaction of a color-octet dipole, 
formed by the gluon at x and the quark-ant i- quark pair at 0. Adding the diagram complex conjugate 
to Fig. 131b we obtain the following contribution to the cross section 

= - [ d^d'xd'y^^^ 
d^k dy n J (2vr)^ vr x^y 

Combining the terms in Eqs. ()192j) and ()193p yields 

da^^ 1 f ,2, ,2 ,2 1 asCp^-y 



- I d%d'^x d^y ^ ^'^^ - ' - e-'^-^^~y) fl - e"^" QloM^M^)/^ + i _ e'^ Q?oi°(Vly|A)/4\ _ 

(193) 



d^k dy TT I ^ (27r)^ n x^y'^ 



(1 



Q2^1n(l/|x|A)/4 _ g-f Q2„ln(l/|y|A)/4 _^ ^-{x-yf Q^^, ln{l/|x-y|A)/4^ _ j^^g^^j 



Eq. p94|) gives us the single gluon inclusive production cross section for a scattering of a quark on a 
nucleus in the quasi-classical approximation, resumming multiple rescatterings of Fig. |^ Before we 
proceed to derive the properties of Eq. p94|) let us first show which diagrams it corresponds to in a 
different gauge, the light cone gauge of the nucleus, and demonstrate an interesting duality between 
multiple rescatterings of Fig. |^ and the non-Abelian Weizsacker- Williams wave function of the nucleus 
on the light cone from Fig. first observed in [^2] ■ 



Initial Versus Final State Interactions 



Let us now consider pA scattering in = gauge, which, in our convention, is the light cone gauge 
of the nucleus. We also boost the system into a frame where the proton is at rest and the incoming 
nucleus is ultrarelativistic. The analysis of pA process in the light cone gauge of the nucleus is sometimes 
referred to as Ap scattering [5^ . 

The light cone gauge diagrams contributing to the gluon production cross section in proton-nucleus 
collisions are depicted in Fig. |221 Again we are going to perform the calculation in the framework of the 
light cone perturbation theory |135t 113b'] . Similar to A~ = gauge case considered above the incoming 
nucleus can emit a gluon in its wave function either before or after the collision with the proton. The 
one gluon light cone wave function of an ultrarelativistic nucleus is given by ^^'^{x) ■ e, with the 
non-Abelian Weizsacker- Williams field of the nucleus given by Eq. (IHU)) of Sect. 12.2.21 with suppressed 
x~ dependence (taking x^ oo) and with e the gluon polarization vector in A'^ = light cone gauge. 
The correspondence between the light-cone wave function and the classical gluon field was shown above 
to be true at the lowest order, which could be seen by comparing x~ —>■ oo limit of the lowest order 
term in Eq. ()80p with the light cone wave function in Eq. ()132j) . Diagrammatically the light cone wave 
function corresponds to the same set of diagrams as was depicted in Fig.^J The fields of the nucleons in 
the nucleus "gauge rotate" the Weizsacker- Williams field of one of the nucleons [Sj. Therefore, here we 
have strong reasons to assume that this light cone wave function can be obtained from the non-Abelian 
Weizsacker- Williams field of Eq. (jHUj) by taking x~ ^ oo in it and multiplying the resulting field by 
e. The interaction of this non-Abelian Weizsacker- Williams (WW) wave function with the proton can 
only be by the means of single or double gluon exchanges, since having more gluons in interaction of 
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the WW gluon field and nucleons from tlie wave function of Fig. with the same proton would bring 
in extra powers of not enhanced by A^^^, which would be subleading in terms of resummation of the 
parameter a 





Figure 32: Diagrams contributing to the same process of gluon production in pA collision as in Fig. |^ 
but now working in A'^ = (nucleus) light cone gauge in the quasi-classical approximation. 

Before colliding with the proton the nucleus can develop the Weizsacker- Williams one gluon light 
cone wave function which then interacts with the proton by means of one or two gluon exchanges, 
according to the rules of the quasi-classical approximation. The square of the graph corresponding 
to this scenario is shown in Fig. IH^. As in Fig. |^ the direct interactions of the proton with the 
nucleons in the nucleus cancel through the real-virtual cancellation leaving only the interactions with 
the gluon line. One may notice that the final state interactions, where another WW gluon merges with 
the outgoing gluon in Fig. 1221 are left out in the diagram of Fig. |22K, but as we will show below we do 
not need them to reproduce the contribution of the graph in Fig. HWi. which implies that they cancel. 

The second possible scenario corresponds to the case when there is no gluon in the nuclear wave 
function by the time the collision happens and the gluon is emitted by the nucleus after the interaction 
with the proton. Then the nuclear wave function without an emitted gluon corresponds to the fields of 
the nucleons "rotating" the current (quark line) of one of the nucleons in the nucleus. This is shown 
on the right hand side of Fig. 132b. The nucleon then interacts with the proton by exchanging one or 
two gluons with it. After that the nucleus can emit a gluon to be produced in the final state. Another 
possibility which is not shown in Fig. 132b but which contributes to the gluon production corresponds to 
the case when the Weizsacker- Williams gluon is present in the nuclear wave function by the time of the 
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collision, similar to Fig. IH^. but after the interaction with the proton the gluon merges into the quark 
line of one of the nucleons, which later re-emits the gluon. We could not find an a priori argument 
prohibiting an emission of the whole Weizsacker- Williams field after the interaction. However, as we 
will see below one needs to emit only one gluon to be able to reproduce Eq. ()194|) . The square of the 
diagram on the right hand side of Fig. 132b is zero since the interactions cancel due to real-virtual 
cancellation The only contribution we get from it is the interference term depicted in Fig. 132b. 
There on the left hand side we have the same diagram as in Fig. l32l \ except that now interactions of 
the proton with the "last" nucleon in the nucleus do not cancel. One can show that the diagram of 
Fig. 132b provides us with the contribution equal to that of the graph in Fig. 131b. 

Let us now calculate the diagrams in Fig. 1^21 The contribution of Fig. |22K can be obtained by 
convoluting the correlation function of the fields on both sides of the cut with the gluon-proton interac- 
tions amplitude, which can be obtained from Eq. by replacing Cp ^ Nc'vn the gluon distribution 
function from Eq. (j3H|) . The result yields 

Employing the correlation function of two WW fields from Eq. ()1U0|) in Eq. ()195|) and defining new 
variables z_ = x — y and b = ywe obtain 

We will need the following mathematical formula 



= - / O -^/" T (197) 
z^A^ Tc J y [z + yY 

where the y integration is cut off by 1/A in the infrared near y = Q_ and y = —z_. Inserting Eq. ()197p 
into Eq. ()196|) and comparing the result to Eq. p92j) one can see that 



(Pkdy (Pkdy ^ ' 

Thus we have shown that the contribution of the diagrams in Fig. is equal to the contribution of 
the diagrams in Fig. IMK. 

The calculation of the graphs depicted in Fig. 132b is a little more complicated. We refer the 
interested reader to ^2], where the estimate is done in detail. After performing the calculation one 
demonstrates that the contribution to the inclusive gluon production cross section from the diagram 
in Fig. 132b is equal to the contribution of the diagrams in Fig. 131b. We can now conclude that the 
diagrams in Fig. |221in the = gauge reproduce the production cross section from Eq. (|194p . 

Comparing Fig. ED with Fig. |221 we see that, indeed, depending on the choice of gauge, different 
sets of diagrams become dominant in the production cross section. Diagrammatic representation of the 
process is gauge dependent. In the light cone gauge of the projectile, A~ = 0, the gluon production in 
proton-nucleus interactions is dominated by multiple rescatterings, as we see from Fig. |^ These mul- 
tiple rescattering are (almost) instantaneous, happening only during the passage of the proton through 
nuclear matter. One can refer to them as final state interactions. On the other hand, the same process 
of gluon production looks different in the light cone gauge of the nucleus, A'^ = 0, as shown in Fig. |221 
The multiple rescatterings of Fig.l31lbecome incorporated into the WW wave function of the nucleus in 
Fig. 1221 Therefore, final state interactions from one gauge become initial state interactions in another 
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gauge. This duality (or dichotomy) was first observed in 
kT-factorization 



It is interesting to show that one can recast Eq. ()194|) in a fcr-factorized form. Following |62] we 
first use the amplitude from Eq. p03p to rewrite Eq. ()194|) as 

= - f d% £x £y -i-^ "^^^ f e-^^<^-y) \NG{x,b,0) + NG{y,b,0)-NGix-y,b,0) , 
d'^k ay rr J l^^j ^ a;^?/"^ l — — J 

(199) 

where we have for simplicity assumed that Nq is a slowly varying function of the impact parameter b. 
(In principle gluon dipoles x, y and x — y have slightly different impact parameters, the spread in which 
we neglected in Eq. (I199p by putting them all equal to b.) Integrating over y in the first term in the 
brackets of Eq. ()199|) . over x in the second term in the brackets of Eq. ()199p . and over, say, x in the 
third term in the brackets of Eq. p99|) (keeping x — y fixed) yields 



daP^ 1 asCp 



d^k dy tt 



d%d'^ze-'^-^ 



z-k , 1 
2 1 ——TT — m ■ 



NG{z,b,0), (200) 



where z_ was chosen to replace x, y and x — y in each of the terms in the brackets of Eq. (jl99|) 
correspondingly. We also use zt = \z\. In arriving at Eq. (j200|) we have used Eq. (jl97|) along with 

"^'^ e-^-^ = -—4- (201) 



(27r)2 z2 27rfc2- 

Using the fact that Ng{z = 0, 6, 0) = we write Eq. ^M)^ as 

Eqs. fl59p and (j3H|) allow us to derive the forward scattering amplitude of a gluon dipole on a single 
quark (or proton), which we will denote ug, such that 



/ d%'nG{z,^,y = 0) = vra^^Mn^. (203) 
J Zt a 



ztA 

Eq. ()203|) corresponds to the two gluon exchange interaction between the dipole and the proton. With 
the help of Eq. fM^ we rewrite Eq. as ^ 



daP^ Cp 1 



d^k dy Us n {2n)^ k^ 



J d^Bd'^bd'^zVlnG{z,b-B,0)e-'^-^VlNG{z,b,0). (204) 



Now 5 is the impact parameter of the proton with respect to the center of the nucleus and b is the 
impact parameter of the gluon with respect to the center of the nucleus as shown in Fig. EHl 
Defining the unintegrated gluon distributions for the nucleus 

Mx.k'') = — ^ / d% d'z e-^^-^ Vl NgU b,y = \n 1/x) (205) 
as {2Tcy J 

and for the proton 

0p(x,^2) = [ d%d^ze-'!^-^VlnGiz,b,y = \nl/x), (206) 
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Figure 33: Gluon production in pA collisions as seen in the transverse plane. To make the picture easier 
to read the gluon is placed far away from the proton which is a highly improbable configuration. 



we transform Eq. (|2U4p into 



dgP^ 
(Pk dy 



2a, 1 



Cf k 




iM^i^y = 0)4>A{k-q,y = o). 



(207) 



Eq. ()207j) is the well-known /Ey-factorization expression for gluon production cross section |168[ 
I139[ll6l't l9()j. It has been proven in |168j for a gluon production from a single BFKL ladder of Fig.lH A 
graphical representation of /cT-factorization implied by Eq. ()207p is shown in Fig. |^ fc^-factorization, 
as shown in Eq. ()207p . is the statement of separation of the inclusive cross section into a convolution of 
unintegrated gluon distribution functions of the target and of the projectile with the square of Lipatov 
vertex for gluon production from Eq. (|22p . This is what is shown in Fig. |^ for pA collisions at hand. It 
is very surprising to obtain this factorization formula after resumming multiple rescatterings of Fig. 1^ 
which seem to explicitly violate factorization of Fig. |^ It may appear that the representation of gluon 
production in the light cone gauge of the nucleus, shown in Fig. |221 niay help cast the process in a 
factorized form of Fig. EU Unfortunately that is not the case: the reason is that the unintegrated 
distribution function used in Eq. (j207|) . given by Eq. (j205|) . is different from the unintegrated gluon 
distribution in Eq. (jl02p . which we obtained by calculating the correlator of two WW fields. The 
"factorization" of Fig. appears to require the nucleus to be described by the gluon distribution of 
Eq. ()102j) (see Fig. 135b). which is the correct correlator of two WW fields which provides us with the 
number of gluon quanta in the wave function. Unfortunately Eq. ()204j) demands that one uses the 
gluon distribution from Eq. ()205|) in the factorization of Eq. (120 7j) instead of Eq. ()102|1 . Therefore, the 
y4+ = light cone gauge representation of the interaction in Fig. 021 does not appear to lead to Eq. ()207j) 



Figure 34: /cT^-factorization in pA collisions, as implied by Eq. ()207p . Unfortunately in no known gauge 
do the contributing diagrams look as shown. 
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Figure 35: A. The definition of unintegrated gluon distribution relating it to tlie gluon dipole cross 
section from Eq. ()205|1 . Tlie exclianged gluon lines can connect to either gluon in the dipole. B. The 
non-Abelian Weizsacker- Williams gluon distribution obtained by calculating the correlator of two WW 
fields, given by Eq. (fTTl^ . 

with the gluon distribution from Eq. ()2()5j) . In general, no gauge is known to the authors in which the 
dominant diagrams responsible for gluon production could be cast in the form shown in Fig. 

Finally, we would like to point out the difference between the two unintegrated gluon distribution 
functions in Eqs. ()205|) and ()102|) . We show the definitions of both distributions in Fig. ESI As was 
mentioned above, the WW gluon distribution from Eq. ()l()2j) is obtained by constructing a correlator 
of two WW gluon fields of Eq. (jHOI)- It is shown in Fig. 135b. It is this distribution which counts the 
number of gluon quanta and, in that sense, is the true unintegrated gluon distribution function. (The 
distribution from Eq. ()102|) should not be confused with the one we defined in Eq. fjl39|) to connect the 
BK and GLR equations: the former is defined in terms of a gluon dipole, while the latter employs the 
quark dipole.) The gluon distribution from Eq. (j2(J5|) was required to recast the gluon production cross 
section (|194|) in the fc^-factorization form of Eq. (|2U7|) . However, its diagrammatic representation exists 
and is shown in Fig. There we consider the scattering amplitude of a gluon dipole on a nucleus. 
To define the gluon distribution from Eq. ()205|) one has to require that the interaction happens only via 
a two-gluon exchange, with some complicated gluon distribution in the nuclear wave function. It is not 
clear whether diagrammatic representation of the dipole-nucleus interaction from Fig. ESK is achievable 
in any known gauge. The bottom line is that it is rather puzzling that the /cT-factorization of Eq. (j2(J7|l 
arises for the single gluon production cross section in pA, and it is also not clear why the unintegrated 
gluon distribution used in that formula should be given by Eq. ()205|) . which the calculations lead to, 
instead of the physical gluon distribution of Eq. ()102|) . 

Eq. ()205|) . along with Eqs. ()206p . allows one to generalize the gluon production cross section ()194|) . 
which was originally derived for quark-nucleus scattering, to the case of proton-nucleus scattering. This 
could be done by using the gluon dipole-proton amplitude uq instead of the dipole-quark one from 
Eq. (|2(J3|) . However, the resulting two gluon distribution functions are different only by an overall 
prefactor, which is not important for our discussion below. 

3.1.2 Nuclear Modification Factor: Low-p^ Suppression 

To analyze the properties of the particle production it is convenient to construct nuclear modiGcation 
factor, which is defined experimentally as 



cPk dy 



(208) 



Ncoll 



d'^k dy 
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where Ncoii is the number of proton-nucleon colhsions in the pA scattering process, dN^^/(Pkdy is the 
particle multiphcity in pA scattering and dN^^/cPkcly is the particle multiplicity in pp. BF^ can be 
expressed in terms of particle production cross section as (for subtleties of this redefinition see |lti9j ) 

R'\k.y) = -f^. (209) 

d?k dy 

In this and the following Sections we will study given by the quasi-classical gluon production 
cross section derived above in Sect. 13.1.11 BF^ for quark production, that will be derived later, can be 
constructed by analogy. To construct gluon BP^ using Eq. ()2()9|1 we will need the gluon production cross 
section ()194j) , which we will use in the form given by Eq. ()2()()|1 , since it is more convenient for deriving 
various asymptotic regimes. Eq. ()2U9p also requires us to know the gluon production in pp collisions, 
which could be obtained by expanding Nq from Eq. ()103|) to the lowest non-trivial order, corresponding 
to rescattering on a single nucleon in the nucleus, putting A = 1 in it and using the resulting expression 
in Eq. ()2()()|1 . The result yields 



with A some infrared cutoff. In Eq. ()210|) we have included the factor of A on the left hand side to 
make it ready to use in Eq. ()209p . 

In this Section we are interested in the behavior of R^^{k, y) at low values of the gluon transverse 
momentum kx, namely for k^'^Qso- Analyzing Eq. (j2(J(J|) we note that low kx^Qso corresponds to 
dominance of large dipole sizes 2;^ > 1/ Qsq in the integral. For > 1/ Qsq the cross section is likely to be 
black, and we can approximate Nq in Eq. ()200j) by NG{z_,b, 0) ~ 1. Integrating the resulting expression 
over z we obtain 

Before we proceed, let us point a very important property of the classical gluon production cross 
section (|194j) . Eq. (|21U|) in fact reflects a high-fcy asymptotics of the gluon production cross section: at 
high-fcr it scales as ~ 1/k^, in agreement with the lowest order perturbative calculations [1312312211201 
IMJSni- This result is badly infrared divergent, and demonstrates a problem of lowest order perturbative 
calculation of gluon production, shown in Fig. |Sj if one wants to obtain the integral production cross 
section da/dy from it, one would have to introduce an infrared cutoff A for fcy-integration. The result 
would depend on the cutoff as ~ l/A^, making the resulting da/dy very sensitive to the value of the 
cutoff and demonstrating that the lowest order cross section is dominated by non-perturbative infrared 
effects, which we do not know how to account for theoretically. On the other hand the full classical cross 
section of Eq. ()194j) is (almost) free of such problems. Eq. ()211|) describes the low-k^ asymptotics of the 
cross section in Eq. ()194|) . As one can see it leads to the cross section scaling as ~ 1/k^. Thus, similar 
to the classical gluon distribution of Sect. 12.2.21 fsee Fig. ITTj) . the multiple rescattering effects regulate 
the infrared singularity, changing it from ~ l//c|, to ~ 1/k^, and making the resulting cross section less 
infrared divergent. The residual infrared divergence ~ 1/A;|. is due to the fact that we have assumed that 
saturation behavior with multiple rescatterings takes place only in the nuclear wave function: there is 
no saturation in the proton wave function in the approach used here. The case of multiple rescatterings 
in both the target and the projectile is usually studied in the context of nucleus-nucleus (AA) collisions. 
There the residual ~ l/fc|, divergence is regulated completely EH E] . 

Using Eqs. (ETT|l and (ETn|) in Eq. (1^ yields 
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where, for simplicity, we assumed again that the nuclei are cylindrical in the beam {z) direction, such 
that the impact parameter integration over b in Eqs. (j211|) and (j21(J|) would give a trivial overall factor 
of the nuclear transverse area S±, which cancels in Rp^. From Eq. (j212|) we see that [5 ^ IMl IS^ lU^ 1^ 

ESlEaiHHllHni 

RP^{k,y)<l, kT<Qso- (213) 
Therefore we conclude that classical gluon production in pA leads to suppression of Iow-Zct gluons. 

3.1.3 Nuclear Modification Factor: Cronin Effect 

Our goal now is to study the nuclear modification factor R^^ for high-fej- particles with > Qso- 
Starting with gluon production again we note that, in the quasi-classical case one can prove the following 
sum rule [221 



d'ke 



daP^ 



A / d^ke 



daPP 



(214) 



d'^k dy J d^k dy^ 

which could be obtained from Eq. ()200|) using the fact that, for the quasi-classical amplitude ()103|) and 
the two-gluon exchange amplitude from Eq. ()203p the following relation is true 

hm {[v^iVG(z, 6,0)1 -^1/=^ \vlnG{z,hM] = 0. (215) 

Eq. (j215|) and the sum rule of Eq. (j214p are valid in the quasi-classical case of McLerran-Venugopalan 
model only. (When quantum small-x evolution is included, Eq. ()215p ceases to be valid.) 

The sum rule fl214|) insures that if the quasi-classical gluon production cross section in pA collisions 
is, in some region of kx-, smaller than A times the gluon production cross section in pp than there should 
be some other region of k^ in which their roles are reversed. For Rp^ defined in Eq. (j209|) that means 
that if, in some region of /c^, it is less than 1 there must be some other region of kx in which it is 
greater than 1. Since in the previous section we have proven that Rp^ < 1 for kr < Qso, the sum rule 
fl214|) demands that RP"^ > 1 for some region of kr in the kx > Qso interval. Enhancement of particle 
production in pA collisions leading to Rp^ > 1 was originally observed by Cronin et al in ^3] and is 
usually referred to as Cronin effect. Eqs. (j214p and (|213p prove that the cross section (jl94p leads to 
Cronin enhancement of produced gluons in pA collisions. 
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Figure 36: The ratio RP"^ for gluons plotted as a function of k^/QsQ in the quasi-classical McLerran- 
Venugopalan model as found in [221 IH2l EH EHl E2] ■ The cutoff is A = 0.2 Q,. 
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To analyze the behavior of RP^ for kx > Qso we follow the approach originally introduced in | 
and find the first correction to the high-Zc^ asymptotics of Eq. (jl94j) given by Eq. (j21(jp . Expanding Nq 
from Eq. (jl(J3p to the second non-trivial order, substituting the result in Eq. (|2UUp and performing the 
integration over z_ we obtain 
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(216) 



Substituting Eq. ()216|) along with Eq. ()210p into Eq. ()209p and keeping only the leading logarithmic 
(ln(A;|n/A2)) terms in the parentheses of Eq. ()216|) yields for a cylindrical nucleus 



RP^ikj 



1 + - 



2 k' 



ln^ + ... 



k. 



00. 



(217) 



Eq. (|217|) indicates that approaches 1 from above at high k^, which is typical of Cronin enhancement 
|84j . We therefore conclude that in the framework of the quasi-classical approximation employed in 
arriving at Eq. ()194j) the ratio is less than 1 at small kx^Qso and displays Cronin enhancement 
at high kr > Qso- This conclusion has been reached in|HainillH3inaEaiHnilHZllHHllH3- 

To summarize let us note that, if one neglects the logarithms in the exponents of Eq. ()194|) as a 
slowly varying functions of transverse separations, writing for instance x^\n{l/xTA.) x^, it would 
become possible to perform the x and y integrations in Eq. ()194|1 exactly, obtaining [^21 lElj 
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(218) 



where Ei(a;) is the exponential integral. Expanding Eq. (j218|) would give the gluon production cross sec- 
tion in pp, which can be obtained from Eq. ()210j) by dropping ln(fc|./A^) in it. The nuclear modification 
factor resulting from Eq. ()218p is 
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The ratio R^^^kx) is plotted in Fig. EH for A = 0.2 Q^q. It clearly exhibits an enhancement at high-fc^^ 
|HlinilH3iniinaEnilHIllHHlEni typical of Cronln effect [H. The height and position of Cronin 
maximum are increasing functions of Qg, and, therefore, of collision centrality 



3.2 Gluon Production Including Quantum Evolution 

Here we are going to show how to include the effects of quantum BK evolution equation from Sect. 12.31 
into the gluon production cross section derived in Sect. 13.11 We will then analyze the impact of small-x 
evolution on the nuclear modification factor R^"^, demonstrating that it eliminates Cronin enhancement 
leading to suppression of gluon production at all transverse momenta kx- 



3.2.1 Gluon Production Cross Section 

Starting with the gluon production cross section, let us first solve the problem of gluon production in 
deep inelastic scattering (DIS) instead of pA collisions. The advantage of DIS is that we know explicitly 
how to include the non-linear small-x evolution into the dipole wave function, as shown in Fig. |^ and 
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Figure 37: Gluon production in DIS in the quasi-classical approximation. The produced gluon may 
be emitted either off the quark or off the antiquark: lines both in the amplitude and in the complex 
conjugate amplitude. Only one connection is shown. 

given by Eq. ()135|l . The answer we will obtain this way would be easy to generalize to the case of pA 
collisions. 

The problem of including the effects of BK evolution into the quasi-classical gluon production cross 
section from Eq. ()194|) could be thought of as producing one of the gluons in the non-linear cascade 
of Fig. 1^ However, the cascade we needed for calculation of the total DIS cross section is somewhat 
different from what we need now for the inclusive cross section. The difference is in the fact that the 
cascade of Fig. |^ develops from early times, with light cone time r = x~ = —oo, till the interaction 
with the nucleus at r = 0.^ Since Fig. |21] represents the forward scattering amplitude, the evolution 
of the cascade after the interaction, during the times between r = and r = +oo, is restricted to 
one simple case when the system returns back to the initial state and nothing is produced. This is 
the definition of the forward amplitude and this is all we need to calculate the total cross section. To 
calculate the inclusive cross section we need to square the production amplitude, as was shown in, for 

^Remember, all the multiple rescatterings with the nucleus in Fig. 1111 happen very fast, hence the whole multiple 
rescattering process can be viewed as instantaneous when compared to the time it takes to develop the cascade of 
gluons/dipoles in Fig.|^ therefore we denote the time at which all the multiple rescatterings take place as t = 0. 
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instance, Fig. |^ for the quasi-classical case. There, already in the amplitude, we had to include the 
diagrams where the incoming proton would "cascade" , generating a gluon both before the interaction 
in the time interval r = — oo ... and after the interaction, at r = . . . + oo. Therefore, to construct a 
gluon cascade making sure that we always have the gluon we are trying to produce in the final state, we 
have to include gluon emissions bringing in logarithms of energy both at r = — oo . . . and r = . . . + oo 
in the amplitude and in its complex conjugate. 

To do that let us first generalize the quasi-classical expression from Eq. (|194|) to the case of DIS. 
The relevant diagrams are shown in Fig. |S3 They are similar to the case of quark-nucleus scattering 
pictured in Fig. |^ The major difference is that now we have an incoming qq dipole instead of just a 
single quark. The produced gluon can, therefore, be emitted off of either quark and anti-quark lines on 
both sides of the cut. In Fig. EH we show only one particular way of emitting the gluon. Diagrams A 
and B in Fig. |S7| correspond to diagrams A and B in Fig. |^ Another difference now is the diagram in 
Fig. EZP, which used to cancel through real- virtual cancellations in quark-nucleus scattering case and 
was not even shown in Fig. EH is now non-zero and has to be included. This is due to the fact that 
moving a t-channel exchanged gluon across the cut in Fig. EZP would change the color factor of the 
diagram, thus preventing the cancellation. Adding all the diagrams in Fig. EH and summing over all 
possible emissions of the gluon off the quark and anti-quark lines, we write for the gluon production 
cross section in the dipole-nucleus scattering |132| Elj 



X 



-Ng (^Xi-Xj, 
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(220) 



where, just like in Sect. I2.3.2[ the quark and the anti-quark have transverse coordinates Xq and x_i 
correspondingly. The rest of the notation of Eq. ()220|) is explained in Fig. EH Ng in Eq. ()220|) is taken 
from Eq. ()103p . and now we put the correct impact parameters for all the dipoles. 

To include the effects of small- a; evolution we should add to the graphs in Fig. EZI the diagrams 
with more gluon emissions before and after the interaction. Let us denote the rapidity of the target 
nucleus as and the rapidity of the incoming qq dipole as Y . The produced gluon would have rapidity 
y. Following [33] we will divide all possible extra gluon emission into two categories: the gluons can 
have rapidities larger (harder) or smaller (softer) than y. 

First we analyze the emissions of harder gluons. Some of the relevant diagrams are shown in Fig. EH 
Since we are trying to apply quantum evolution from Mueller's dipole model jl9t 1201 1^ we are again 
working in the large- iVc limit. Similar to Fig. |2Slthe gluons are denoted by double lines which do not 
connect to any particular quark or anti-quark lines in the dipole in which a particular gluon is emitted: 
this denotes summation over both connections on both sides of the cut. In Fig.EHthe gluon 3 is the one 
that we measure in the final state, the gluon 2 is a harder gluon emitted before gluon 3. Solid vertical 
lines in Fig. IHHl denote the r = +oo final state, while the dashed vertical lines denote the multiple 
rescattering interactions with the target from Fig. II II taking place at r = on both sides of the cut. 

In Sect. 12.3.21 we have demonstrated that successive emissions of progressively softer gluons before 
the interaction (for r = —oo ... 0) give the leading logarithmic contribution, in the sense of giving 
a factor of Y for each factor of as resumming the parameter asY , as is accomplished by Eq. ()135|) . 
Interestingly enough, this ordering has to be reversed for emissions after the interaction (for r = 
0... + oo): to obtain the leading logarithmic contribution (a^l^), softer gluons have to be emitted 
before harder ones at r = . . . + oo |34|- More details on how this happens could be found in j34j . 
Indeed, if emission of a harder gluon happened after the emission of the softer one, the harder gluon 
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Figure 38: Emission of a harder gluon in the dipole evolution, as considered in Gluons are denoted 
by double lines in the large Nc limit. The produced gluon is marked by a cross. Dashed lines comprise 
multiple rescatterings with the target from Fig. which are instantaneous on the time scale of the 
dipole evolution. 



can not be emitted off the softer one: it is emitted by hard quark lines of the original dipole [21] . Since 
the diagrams E-I in Fig. EHl violate such inverse ordering, they can be discarded. Only the diagrams 
A-D in Fig. EH contribute. Those are real contributions to the dipole evolution leading to production 
of the dipole in which gluon 3 was emitted. Note that all of these surviving diagrams have early 
time (r = — oo ... 0) emissions only: they correspond to standard dipole evolution of Fig. |^ A similar 
analysis shows that out of the rest of the diagrams with possible emissions of gluon 2, the ones that bring 
in leading logarithmic contribution are the ones which give virtual corrections to the dipole evolution at 
light cone times r = — oo . . . [Mj. Emissions of other harder gluons do not change the conclusions: the 
surviving diagrams are the ones which contribute to dipole evolution leading to production of the dipole 
in which gluon 3 was emitted. Such evolution is described by the quantity ni{xjQi, Xjq,ii, B_ — b,Y — y), 
originally introduced in fOllinij which has the meaning of the number of dipoles of size Xj^/y at rapidity 
y and impact parameter b generated by evolution from the original dipole having rapidity Y and 
impact parameter 5 . This quantity obeys the dipole model analogue of the BFKL equation j2Z] , which 
is just the linear part of Eq. (fnH|l [Wl 1^ 



dni{xj^^,Xj^,^,,b,y) asNc 



dy 



2 7r2 



j (fx2 



X, 



01 



^2 2 

•^20 -^21 



"-1(^02, ^o'i',^+ ^X2i,y) + ni{xu,Xj^'i',b+ ^X2o,y) 



- ni{xo^,Xo,y,b, y) 



(221) 
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with the initial condition 

"-i(^oi,^o'i',^,y = 0) = 6{xqi- XQ,-^,)6{b). (222) 

The inclusion of the effects of harder gluons in Eq. ()22()|1 is then accomplished by replacing the cross 
section in Eq. ^TM by |M1 

d^kdyd^B ^'^'^^J'^'^'^''"" M^u-.'v,B-b,Y-y) ^^^^(^oaO- (223) 

Now we have to resum emissions of softer gluons having rapidity less than the rapidity of the 
produced gluon y. Such emissions are demonstrated in Fig. EHl where now the gluon 2 is softer than the 




Figure 39: Diagrams including one softer gluon produced after the emission of the measured gluon |34j . 

gluon 3. Fig. ISnionly has emissions where gluon 2 interacts with gluon 3 and with either the quark or 
the anti-quark lines in the dipole. As was shown in such diagrams cancel pairwise via real- virtual 
cancellations: 

A + D = 0, B + C = 0, F + G = 0, E + H = 0. (224) 

The diagrams where gluon 2 connects only to the quark and the anti-quark lines cancel analogously 
pij . The only remaining diagrams are where the gluon 2 is emitted and reabsorbed by gluon 3 only. 
Using the crossing symmetry [167|| we reflect the line of gluon 3 into the complex conjugate amplitude 
(denoted by 3') with respect to the cut, obtaining a gluon dipole in the amplitude, just like we did in 
Sect. l3.1.1l to obtain Eq. ()194j) . Then it becomes manifest that emission of gluon 2 can be thought of as 
one iteration of dipole evolution in the gluon dipole 33'. Successive emissions of even softer gluon would 
not modify this conclusion. Therefore, to include the effects of softer gluon emissions in Eq. ()220j) we 
have to replace [SI] 

NG{x,b,0) ^ NG{x,b,y) (225) 

on its right hand side, where now NG{x,b,y) is the amplitude of a gluon dipole interacting with the 
nucleus including the small-x evolution of Eq. p35|l . In the large- A^^c limit Nq can be expressed in terms 
of the quark dipole amplitude from Eq. ()135|1 as 

Ncix^ky) = 2N{x,b,y)-N\x,b,y), (226) 

since a gluon dipole in the large- Ac limit can be thought of as a pair of quark dipoles, with either one 
of them or both quark dipoles interacting. 
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Combining the prescriptions for including quantum evolution for hard ()223p and soft 
we obtain the following expression for inclusive gluon production cross section in DIS ||34i 



(Pk dy(PB 



feoi) 



SxQiy ni(xoi, Xq/i/, 5 -h,Y -y) 



(2vr) 



emissions 



i,j=0' 



_iy+j -1 - 



-y ,-yj 2 /y 2 

±.1 \±.2 d^jl 



Ng\Zi~ Xj, -Ui + Xj),y] + 



+ Ng (^Z2 - Xi, ^(^2 + Xi), y^ - Ng (^zi - Z2, ^Ui + ^2), - Ng (^Xi - Xj, ^{Xi + Xj),y^ 



(227) 

Eq. ()227|) provides us with the single inclusive gluon production cross section for scattering of the 
quark-ant i- quark dipole 01 on the target nucleus, which includes all multiple rescatterings (powers of 
a^A^^^) and small-x evolution corrections (powers of asY and 




Figure 40: Interpretation of the inclusive gluon production in pA collisions from Eq. 
fan diagrams. 



in terms of 



If one associates Eq. p35|) with summation of fan diagrams from Fig. El than one can represent 
Eq. ()227p in the way shown in Fig. 1321 There the interaction between the quark-ant i- quark dipole and 
the target is mediated by exchange of fan diagrams, where we trigger on one of the produced gluons, 
which is denoted by the cross. The part of the diagram between the produced gluon and the projectile 
dipole is described in Eq. (j227|) by the quantity ^i(xoi, Xq'i/, S — 6, F — y), which satisfies linear BFKL 
evolution equation ()22H) . In terms of ladder and fan diagrams, that corresponds to having just a single 
ladder exchange between the dipole and the produced gluon, as shown in Fig. HHl In principle one could 
draw gluon production diagrams where the ladder between the dipole and produced gluon would split 
into several ladders interacting with the target. However, Eq. (j227|) implies that such graphs cancel. 
The same cancellation is expected if one applies AGK cutting rules to QCD diagrams |170t I17H 1^ . 
It is intriguing to see AGK rules work in perturbative QCD in the framework of the dipole model. 
They have also been found to work previously for single diffractive cross sections in DIS calculated in 
the dipole model approach |131j . However, the AGK rules appear to break down in QCD for inclusive 
two-particle production cross sections |1()3| |l04j . 

In the limit of a very large target nucleus the momentum transfer to the nucleus is cut off by inverse 
nuclear radius and is very small. We can therefore take the scattering amplitudes Ng in Eq. ()227p 
at t = 0, which in coordinate space is equivalent to neglecting (factoring out) the impact parameter 
dependence since the shifts in impact parameter in Eq. ()227|) are small compared to the nuclear radius 
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jl66U172] . We, therefore, put all the impact parameters in iV^'s of Eq. 
the steps outlined at the end of Sect. 13. 1.11 we arrive at [SI] 



to be equal to b. Repeating 



(Pk dy 



2as 1 



d^q (Pqqil, y -y) <pA{k - q, y) 



(228) 



with (pAik — q,y) given by Eq. ()205|) and the unintegrated gluon distribution in the dipole given by 

[3^rin3] 



Mq, y) 



1a,CF f d%d''xo're-'!l-^o'i' -^n^^Xoi,Xo,,,,b,y). 



TT 



V2 

^O'l' 



(229) 



Interestingly enough, as we see from Eq. ()228|) . the small-x evolution does not lead to break-down of 
the /cT^-factorization formula (|2(J7|) ! In fact, Eq. (j228p . with the gluon distribution functions given by 
Eq. (|229p and by Eq. (j2(J5j) . only with twice the quark dipole amplitude 2 N instead of Nq, was postu- 
lated as the answer for inclusive gluon production cross section by Braun in [^0] using /cT-factorization 
approach. Still, Eq. ()228|) . obtained in [31], differs from the result of [HO] by the difference in Nq and 
2 (see Eq. fl226p ). which in the standard pomeron language appears to corresponds to emission of the 
produced gluon from the triple pomeron vertex [01] . This difference may become numerically significant 
for low transverse momenta of the produced gluon, as was shown in [HI]. 

Now generalization to the case oi pA becomes manifest: we have to replace the unintegrated gluon 
distribution function of a dipole, 4>qqi'lyy)j by the one for the proton, 4>p{q,y), given by Eq. ()206|) with 
no satisfying the linear part of Eq. ()135|1 



dnG{xjoi,b,y) _ a, Nc 
dy ~ 2 7r2 



J d^X2 



X, 



01 



20 -^21 



nc (^02 ,b+^X2i,y)+nG{x^2,t+^X2o,y)-nG{Xf^i,b,y) 



(230) 



with Eq. ()2()3j) as the initial condition, 
becomes in the case of pA 



To summarize, the inclusive cross section from Eq. ()228j) 



dgP^ 
d^k dy 



2 a. 1 



Cf k^ 



d^q (j)p{q, Y -y) (j)A{k - q, y) 



(231) 



with the unintegrated gluon distributions given by Eqs. (|2(J5|) and (j2(J6|) with the dipole cross sections 
evolved with Eqs. ()135p and ()230p correspondingly. 

Indeed Eq. ()228|) was derived in the large- A^^^ limit. It appear hard (though maybe not impossible) 
to generalize the large- A^^c dipole model to the case of projectile proton: therefore, we can not prove the 
suggested generalization to the proton case. However, this generalization is probably a good approx- 
imation of the exact answer for the following reasons. We showed that /cT-factorization works in the 
quasi-classical case to all orders in Nc (see Eq. ()207|) ). We have also shown that it survives evolution 
corrections (Eq. ()228p ). Therefore, factorization formula ()23H) . which expresses the production cross 
section in terms of gluon dipole amplitudes scattering on the projectile and on the target, is proba- 
bly valid in the general case of including small-x evolution beyond the large-A^^^ limit. The results of 
numerical solution of the full JIMWLK evolution equation for the dipole amplitude were shown to be 
amazingly close to the solution of BK equation [TH] . The linear BFKL evolution analogue in the dipole 
model [I9J gives the same value for the pomeron intercept as the exact BFKL equation [2Z]. Therefore, 
it is probably safe to use Eq. (j231|) for the gluon production in pA, even though strictly speaking it 
should work only for the case when the proton is approximated by a di-quark. Both pA and DIS pro- 
cesses are considered here as scatterings of an unsaturated projectile (proton or qq pair) on a saturated 
target (nucleus), and, therefore, gluon production cross sections for them should be similar. Therefore, 
we will use Eq. (12311) for the analysis outlined below of the gluon production in pA. 
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3.2.2 Nuclear Modification Factor: Double Logarithmic Region 



Let us now analyze what happens to nuclear modification factor (j!2(J9|) when quantum evolution is 
included. To do that, we use Eqs. (j2(J5|) and (j2(J6|) to rewrite Eq. (j231|) as 



daP^ 



(Pk dy as vr (27r)'^ A;^ 



d'B d% d'z Vi ncU, B-b,Y -y) e''^'^ Vi Ng{z, b, y) . (232) 



For simplicity, let us neglect the impact parameter dependence of No throughout the nucleus and of no 
over the transverse area of the proton. Suppressing the impact parameter dependence in the arguments 
of Ng and hq in Eq. (j232j) we rewrite it as 



Sp Sa 



d^k dy as vr {2'kY k 



d'z Vi naU, Y-y) e~''-^ Vl NgU, y) 



(233) 



with 5*^ and 5*^ the cross sectional areas of the nucleus and the proton correspondingly. In this Section, 
and in the one below, we will be interested in gluon production with large kT^Qsiv), such that nonlinear 
effects could be neglected in Ng- Outside of the saturation region we write, similar to Eq. 



NGizT,y) 



du e2"»^(°'^)^(;zrQ. 



\l+2 iu qA 



V 1 



(234) 



where is determined by the initial condition of Eq. (jl()3|l . Similarly, for ug one writes the general 
6-independent solution of Eq. ()23()j) as 



nG{zT,y) 



where, to satisfy the initial condition given by Eq. ()2()3j) . 

1 



cr, 



and 



Substituting Eqs. 
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Ana — . 



(235) 

(236) 
(237) 



and 



into Eq. 



and integrating over z_ yields 



daP^ _ CpSpSA 



d^k dy as (Stt)' 
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X 



Q 



k 



T 



l+2iv 



l+2iu' 



,2«, x{0,v) y+2as Xi^y) (y-y) 



(238) 



with as given by Eq. ()147|) . 

In this Section we are interested in the double logarithmic region (DLA) where kT>kgcom ^ Qs{y)- 
Evaluating the u and z/' integrals in Eq. ()238|) in the stationary phase approximation around the saddle 
point from Eq. (jl54|) with x(0, z/) (and x(0, z^')) given by Eq. (jl53p we obtain 



daP^ 
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2 as (Y-y) In ^ 



(239) 
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In arriving at Eq. ()239|) we made use of Eq. ()236|) . and have also used the fact that near the saddle 
point of Eq. (jl54|) is approximately given by Eq. (j236|) as well |62j . 

To calculate i?^" we also need the gluon production cross section in pp collisions in the same kinematic 
DLA region, which can be obtained from Eq. ()239j) by replacing Sa Sp and Qso Ap in it jFS] 



d?k dy 



DLA as (27r)4 2 a, 1/3/4 (y _ ^)3/4 



exp 



h 



(240) 



Substituting Eqs. 
yields 



and (IMH) in Eq. ^M^ and remembering that Sa = A^/^Sp and Q^^ = A^/^Al 



f In In ^ 
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The behavior of RP^^kx, y) from Eq. 




kT'»Qs 
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at high y is dominated by the exponent in it 



exp 




2J2asy 
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< 1. 



(241) 



(242) 



Since Qso S> Ap, the exponent in Eq. ()242j) becomes less than 1, driving RP^(kT,y) below one as well 
|()2j . We therefore conclude that in the DLA region with /cT^^geom, the nuclear modification factor 
becomes smaller than 1 due to quantum evolution. This indicates suppression of gluon production due 
to quantum evolution. A suppression of Rp^ by small-x evolution effects was originally suggested in 
|61j . though the authors of [61 analyzed particle (gluon) production only in the geometric scaling region 
fcT<fcgeom, which is what we are going to do now. 



3.2.3 Nuclear Modification Factor: Extended Geometric Scaling Region 

In the extended geometric scaling region, QsiX) < A^T^^geom, multiple rescatterings are still unimpor- 
tant and Eq. ()238|) from the previous Section still applies. To obtain the leading behavior of Eq. ()238p 
in the geometric scaling region we will evaluate the i/-integral in it by simply replacing u with z/q from 
Eq. (|163|) in the integrand, as it was done in obtaining Eq. (|18(J|) . We obtain 



daP^ Cp Sp Sa ^a ^ ^ ^^^^y ^2.u, ( ^) 



d^k dy as (27r)3 ' \ kr 

X r du' C^, {1 + 2tu')' 2'^^' /M'""'" g2^.x(oy)(y-.)_ (243) 



T[l-t{uo + iy')] Kkr 
The gluon production cross section in pp collisions is given by 

d'^k dy as[2'Ky^ J-oo T[l — z[h' + u' 



A ^ l+2iiy ^ A V l+2ii/' 

tlp\ (tlP\ ^2asX{0,u)y+2asX{0y)iy~y)_ (^24:4:) 

kj' J xk^ J 
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To evaluate the i^'-integral in Eq. ()243|) and both v and v' integrals in Eq. ()244|) we need to know 
whether we are inside or outside of the extended geometric scaling region for the proton. Recall that, 
while Eq. (j244|) and Eq. (j231|) are valid only outside the saturation region of the proton, such that 
the evolution between the produce gluon and the proton can only be linear, that evolution may still 
be dominated either by the double logarithmic saddle point v^di^a from Eq. p54|) or by the leading 
logarithmic saddle point vq from Eq. p63|) . The dominant saddle point is determined by scale of 
the geometric scaling in the proton, fc|eom) which can be obtained from Eq. ()189j) by using proton 
saturation scale in it, which, in turn, is obtained from Eq. ()l(i5j) by putting A = \. In the problem 
we are considering here, both the atomic number A and the rapidity intervals y and Y — y are large. 
Therefore, in principle, /c|eom can be above the nuclear saturation scale Qs{y)- Since we are interested 
in the extended geometric scaling region of the nucleus, QsiX) < ^T^^geom, we have to consider two 
cases: (a) fc|eom ^ < ^geom and (b) Qs{Y) < fcT<fc^eom- 

(a) ^kT< kgeom If of the produced gluon is outside of the proton extended geometric scaling 

region, then z/'-integral in Eq. ()243|1 and u and z/' integrals in Eq. ()244j) have to be evaluated around 
the saddle point z^Ij^^^ from Eq. (jl54|) . which in momentum space for the case of the proton is 
given by 

- - vrSQ ■ 

Evaluating both Eq. and Eq. using Eq. (j2SSl), substituting the results into Eq. piH]) 

and dropping the slowly varying prefactors, we obtain the nuclear modification factor ^2j| 



1 S, < fcgcom Q 



exp 



sQ 



I kr In' ^ 



A 14C(3)a.2/ 



(246) 



To estimate the value of R^^ in Eq. (j246p in the region k^^^^ ^kT< kgeom we substitute kx = fcgeom 
into Eq. ()246|) with kgeom from Eq. ()189|) . The result yields an asymptotically small value 

RP^ikT = kgeom,y) ~ e-'-^'^^y « 1, (247) 

where we used A = 197 for gold nucleus. For other values of A and for other values of fcy in 
the region fcgeom ~ ^ ^geom one still gets exponential suppression for RP^{kT,y)- Therefore we 
conclude that RP^{kT,y) < 1 in the region fc|eom~^T^fcgeom- 

(b) Qs{Y) < /cT^^geom of the produced gluon is inside the geometric scaling regions of both the 
nucleus and the proton, we evaluate the i/'-integral in Eq. ()243p and u and z/' integrals in Eq. ()244j) 
by putting u = u' = uq in them with uq given by Eq. p63p . Noting that almost all the prefactors 
in Eqs. fl243|) and ()244|) would then be identical, we write for the nuclear modification factor 

-^^A-\[^\ . (248) 



R'\kT.y) 

As QIq = A^/^Al Eq. dUHl) leads to 



R'^{kT,y) 



V 



A-H^. (249) 



Since z/q is just a number, C^^jCI^ does not depend on rapidity and varies weakly with A. The 
nuclear modification factor in Eq. ()249p is, therefore, driven by the power of A in it. Using the 
numerical value of z/q from Eq. ()163j) we obtain [H2j 

R^^{kT,y) ^ A-H"^ = ^ ^_ ^250) 
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This is the suppression of particle production which was originally predicted in j61j . 



We conclude that small- a; evolution leads to suppression of in the whole extended geometric 
scahng region Qs{y) <kT< /Cgeom jEH E21 EH ES] ■ 

3.2.4 Nuclear Modification Factor: Saturation Region 

Above we have shown that for kx > Qs{y) quantum evolution introduces suppression of W^. Here 
we will study the saturation region, kx^Qsijl)- The important question now is to understand what 
happens to Cronin maximum from Fig. IHElas the effects of small-x evolution become important. Since 
the Cronin peak is located at kr = Qs{y), we need to find out what happens with EF^ikx = Qsiv), v) 
as y increases. 

Inside the geometric scaling region {zt > l/A^goom) the dipole amplitude Nq depends only on the 
parameter ^ from Eq. p7H) . Therefore one can write 

NG{zT,y) = / du[zTQs{y)]'~''"'C^ (251) 



with some unknown coefficient. Substituting Eq. ()251|) along with Eq. ()235|) into Eq. ()233p and 
integrating over z yields 



d^kdy 



CfS.Sa r dudu'Cl,Ct{l + 2iuf{l + 2iu'f2''^''^-''''^ r[z(z/ + z/')] 



X (— ^ g2a.x(oy)(y-y)_ (252) 

\Qs{.y)) 

(In arriving at Eq. ()252|) we have assumed that extended geometric scaling region zt > l/^geom dom- 
inates in the integral of Eq. ()233|) : this is explicitly shown in the Appendix of [HI]-) To evaluate the 
z/'-integral in Eq. (j252|) we again need to know whether kx of the produced gluon is above or below the 
scale of extended geometric scaling in the proton k^^^^. This may be an issue when kP^^^ < Qsiy)- As- 
suming for simplicity that A;^eom > Qs{y) we evaluate the i/'-integration in Eq. ()252p by putting u' = uq 
with i/Q given by Eq. p63|) . (The case when the evolution between the proton and the produced gluon 
is dominated by DLA saddle point is examined in [^21 with the same conclusion as we will reach below.) 
Eq. ()252j) becomes 



daP^ 



d^k dy 



= ^^^7^ Cl,CA{^ + 2iv,f2'-^ ,2^.x(o,.o)(r-.)^ (253) 

s(y) as(2vr)^ \Qs[y) J 



where we have defined 

Ca - / dvCt{l + 2^.f2-- I'^'^^'.^^'f (254) 

It is important to note that all the information on the solution of the non-linear evolution equation 
given by Nq is now absorbed into a constant from Eq. ()254j) . which does not depend on rapidity and on 
atomic number of the nucleus, and is thus not important for behavior of RP"^ at high energy /rapidity. 
Similarly, evaluation the v and v' integrals in Eq. ()244j) around vq yields 

d^kdy ~ (27r)3 ^^-oJ (i + 2^z/oj 2 _ 2,^^) ^ • ^^^^) 
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Putting = Qs{y) in Eq. ()255p . and using it together with Eq. ()253|) in Eq. ()209|) yield for the nuclear 
modification factor 



R'^{kT = Qs{y),y) = ^ (l + 2zz.o)-^2-^-o (^^) " e-^-'^^^-'^o)^. (256) 

Remembering that Qsiv) is given by Eq. ()16ip with Qso = A^^^Ap we can recast Eq. ()256p into the 
following form 

R^^kT = Qs{y),y) = ^ (i + 2..o)-^ 2-^-0 r^-^ (257) 

The prefactor of Eq. ()257|) is just a number. The dynamical information is carried only by the power 
of A in it, which, for the value of uq from Eq. (jl63j) gives 

RP^{kT = Qs{y),y) - = A-o-^^^ < 1 for A > 1. (258) 

Eq. (j258p allows us to conclude that small-x evolution tends to reduce the Cronin peak turning enhance- 
ment into suppression. Since the suppression in Eq. ()258p is (parametrically) of the same order as given 
by Eq. ()250p for the extended geometric scaling region, we conclude that both at kx = Qsiy) and 
in the extended geometric scaling region is smaller than 1 at all values of kx in that region, approaching 
the same (energy-independent) limit at very high energies/rapidities. 

Deep inside the saturation region, for kx -C Qs{y), one can use the above techniques to show that 
RP^ there is also suppressed at least at the level of Eq. ()258|l and likely even more than that. We 
refer the interested reader to |62j for details on this derivation. Here we can just argue that saturation 
effects with or without small-x evolution, reduce the ~ infrared-singular scaling of the inclusive 
production cross section of Eq. ()210|) to the more infrared- regular ~ scaling of Eq. 1)2111) . Hence 
the suppression of R^"^ at kx <^ Qs{y), shown in the quasi-classical limit by Eq. ()212|) . will still be valid 
in the case of quantum evolution, though it will be slightly modified by the anomalous dimension uq. 

With the help of Eq. ()258)1 we conclude that at high rapidities/energies the Cronin maximum 
decreases with energy and centrality, with RP^{Qs{y),y) becoming less than 1. Eventually, at very 
high energy, the Cronin peak flattens out and saturates to an energy independent lower limit given by 
Eq. (1258)) ■ which is parametrically suppressed by powers of A. The height of RP^{Qs{y),y) is also a 
decreasing function of collision centrality/atomic number A. 



3.2.5 Nuclear Modification Factor: Overall Picture 

In the three above Sections we have demonstrated that small-x evolution leads to suppression of the 
nuclear modification factor R^^ for the transverse momenta kx of the produced gluon in the double log- 
arithmic region /cr^^geom (see Eq. ()242|) ). in the extended geometric scaling region Qs{y) < kx^kgeom 
(see Eq. ()25Up and Eq. ()247|) ) and inside the saturation region kx^Qsiy) (see Eq. (I258|) ). To summarize 
the effects of quantum evolution, we illustrate how RP^{kT,y) varies with rapidity/energy in Fig. |^ 
which is constructed using a toy model of lU^, that incorporates the correct qualitative behavior of 
RP^kr^y). 

The top line in Fig. |^ is borrowed from Fig. ESI and represents the Cronin enhancement of particle 
production in the quasi-classical approximation. As energy (or rapidity) increases the suppression begins 
to set in as shown by dash-dotted and dashed lines in Fig. |^ The asymptotic energy-independent limit 
of maximum suppression to be reached at very high energies/rapidities given by Eqs. ()250p and ()258p 
is depicted by the lower solid line in Fig. |^ We conclude that as the energy/rapidity increases the 
upper solid line in Fig. would decrease eventually turning into the lower solid line. This result was 
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Figure 41: The ratio RP"^ plotted as a function of kx/Qsiv) for (i) the quasi-classical approximation of 
Fig-EHl which is valid for moderate energies/rapidities (upper solid line); (ii) the toy model for very high 
energies/rapidities from [52] (lower solid line); (iii) an interpolation between the two to intermediate 
energies (dash-dotted and dashed lines). The cutoff is A = 0.3 Qs- The figure is from jHS]- 

obtained by numerical simulations in jHSj and by an analytical analysis similar to the one outline above 
in (see also [64, 92]). 

Indeed Fig. ^2 is not yet ready to be compared to data. In d + Au collisions at RHIC the forward 
rapidity high-p^- particle production receives a substantial contribution from valence quark production, 
along with the gluon production which led to Fig. ^2 On general grounds, since the gluon suppression 
in Fig. ^2 is driven by small- a; evolution, we expect valence quark production cross section to also be 
suppressed at high-p^- To investigate this further and to be able to make quantitative predictions we 
will calculate forward production of valence quarks in the next Section. 

3.3 Valence Quark Production 

We now consider production of valence (high x) quarks in proton-nucleus collisions, due to scattering 
of a valence quark on the nucleus jSZj. In this section, we treat the target nucleus as a Color Glass 
Condensate while the projectile proton is taken to be a dilute system of quarks and gluons in the spirit 
of standard pQCD. This is the appropriate approach if one is probing the large x components of the 
projectile proton wave function which is the case in the very forward proton-nucleus collisions. 

The leading order diagram is shown in Fig. (|42j) . Using the LSZ formalism, the scattering amplitude 




Figure 42: Scattering of a high x quark from a dense target. 



{(l{(l)out\q{p)in) = {out\bout{q)bln{p)\in) 



(259) 
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which can be written as 

{out\bout{q)bL{p)\in)=^ f rf^xdV"*^^'^"'''^^M(g)Mv -m]{out\Tij{y)i:{x)\in)[-i0^ -m]u{p) (260) 



where m is the quark mass and Z2 is the fermion wave function renormahzation factor while u{p) is 
the quark spinor with momentum p. The quark propagator Sp in the background of the classical color 
field is defined by 

{out\Tilj{y)ijj{x)\in) = —i{out\in)SF{y,x). (261) 

The amplitude then becomes (setting Z2 = 1 and {out\in) = 1 since we are working to leading order in 
as and our background field is time independent) 

{(l{(l)out\q{p)in)=if d'^xd^ye-'^P-''-'^-^^u{q)[i -rnlSpiy, x)[-i 0^ -m]u{p). 



In momentum space, the fermion propagator Sp can be written as 

Spiq,p) = {2T^f5\q-p)S%{p)-igSl{q) J ^^i{k)Sp{q + k,p) (262) 

where 4 = ^^lii is the classical background color field, and is the free fermion propagator. It is 
useful to define the interaction part of the fermion propagator from ()262j) as 



5^(g,p) = {27,Y5\q-p)S%{p) + Sl{q)T{q,p)S%{p). (263) 

Substituting ()263|) into the amplitude leads to 

{(l{(l)out\q{p)in) =u{q)T{q,p)u{p). (264) 

This simple relation relates the amplitude for scattering of a quark or anti-quark from the Color Glass 
Condensate and the quark propagator in the background color field of the nucleus. 

The explicit form of the quark propagator in the background of a classical color field was calculated 
in |152| l88]. The interaction part, as defined in ()263|1 is given by 



^(?,P) = (2vr)(5(p -q )l / d^z 



V{z) 



e^d-py^ (265) 



where \^ is a matrix in the fundamental representation which includes the multiple scattering of the 
quark on the nucleus, given by p27p 



V{z) = Pexp 



+00 



(266) 



tg dz-A+{z-,z)e 
J —00 

and t"- are in the fundamental representation. Using ()265|) in the scattering amplitude ()264p gives 

{q{q)out\q{p)^n) = {27f)6{p- - q-)u{q)ru{p) ( d'z[V{z) - l]e^(£-E)-^ (267) 



The presence of the delta function in the amplitude is due to the target being (light-cone) time inde- 
pendent which leads to conservation of the "minus" component of momenta and can be factored out 
(for a rigorous treatment of this using wave packets, we refer the reader to P^ ) 

{qiq)out\qiphn) = {2Tx)5{p- -q-)M{p,q) (268) 
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leading to 



da = l ^(27r)5(2gV-g|)^^(g+)^(27r)5(p--g-)|M(p,g)r. (269) 

One can set the transverse momentum of the incoming quark pt to zero without any loss of generality. 
We then get the differential cross section for production of a quark with transverse momentum qt per 
unit area 

dfji^^i^ 2 

^2^=(^^.^oZe&^) (270) 

where the dipole cross section (which is, via the optical theorem, the imaginary part of the forward 
scattering amplitude) is defined as 



F 

^dipole 



{b,q) = J d^re'^-^^TT,(v\b + r/2)V{b-r/2) - l) ^ = - J d^r e'^'^ N{r,b,Y) (271) 



and we have neglected terms which do not contribute at finite quark transverse momentum. It is trivial 
to extend this calculation to the case of incoming on shell gluons 139, scattering on the dense target. 
The only difference is the presence of adjoint matrix U rather than V which have different color factors 
in the exponent, appropriate for the given representation. It is given by 

^^gA^gX _ 2 ^ 

where the adjoint dipole cross section is defined by 



<yipoieik,q) = / d\e^^-^-jJ—^i:i,{u\b + T/2)U{b-T/2) - l)^ = - | d\ e^T^- Nc^L^kY). (273) 



Our result in ()270[ I272|l can be understood as the generalization of the leading order quark-gluon 
and gluon-gluon hard scattering cross section convoluted with the target gluon distribution function 
in leading twist pQCD to include multiple scattering of the incoming quark on the target. To include 
the effects of quantum evolution in the target at high energies, one should resum large logs of energy 
or 1/x. This can be accomplished by solving the JIMWLK (or BK) equations for the (fundamental or 
adjoint) dipole cross section and using the result in ()270l 1272]) . 

The quantum evolution of the quark-nucleus scattering cross section could have very interesting 
observable signatures in the forward rapidity region. At forward rapidities the nuclear saturation scale 
is large due to small-x evolution. A typical valence quark in the incoming proton would have a transverse 
momentum of order Qs after scattering from the nucleus. This momentum is large and the quark will 
tend to fragment into hadrons independent of the other valence quarks. This would reduce the number 
of baryons produced in the forward rapidity region of pA collisions as compared to pp collisions, while 
increasing the number of produced mesons, so that the meson to baryon ratio could increase in pA 
compared to pp |173j . 

We note that the qualitative behavior of the nuclear modification factor R^a for quark production 
should be very similar to that of gluons. This is due to the fact that the behavior of RpA is determined 
by the properties of the of the underlying (fundamental vs. adjoint) dipole cross section. There may 
however be quantitative differences due to the different color factors for the fundamental and adjoint 
dipoles which would also affect the rate of quantum evolution for the two dipoles. For example, one can 
show that in the classical approximation to quark production, the nuclear modification factor RpA for 
quarks ^U] exhibits an enhancement similar to that of gluons even though the magnitude of enhancement 
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for gluon RpA is larger because of the larger value of the gluon saturation scale due to the color factors. 
As one goes to smaller values of x where quantum evolution in Inl/a; becomes important, the nuclear 
modification factor RpA becomes smaller than unity in analogy to RpA for gluons even though the rate 
of decrease of RpA for quarks is smaller than the rate for the decrease of RpA for gluons. 

To relate this to hadron production in proton-nucleus collisions, one will need to convolute these 
hard cross sections with the quark and gluon distribution functions in a proton and quark-hadron and 
gluon-hadron fragmentation functions. At this level, the distribution and fragmentation functions are 
delta functions of x and z respectively, given by the parton model. One can include DGLAP evolution 
of the projectile parton distributions and hadron fragmentation by considering radiation of a gluon or a 
quark anti-quark pair and then absorbing the resultant collinear divergences into the bare distribution 
and fragmentation functions which lead to DGLAP evolution of the parton distribution functions in 
the proton and the hadron fragmentation functions. The detailed results are given in |174j . 



3.4 Electromagnetic Probes 

In this Section we consider the production of photons and dileptons in proton-nucleus collisions 
inHinSlinHI- photon and dilepton production has the advantage that there is no hadronization involved 
in the final state, unlike quark and gluon production. Therefore, photons and dileptons are a cleaner 
probe of the Color Glass Condensate. Nevertheless, they have the disadvantage of low production rate 
due to the smallness of the electromagnetic coupling vs. the strong coupling involved in quark and 
gluon production. 

3.4.1 Photon Production 

We consider the process q{p)A 'y{k) q{q) X where A stands for a hadron or nucleus which is treated 
as a Color Glass Condensate and the produced quark and photon are real. The incoming quark has 
momentum p while q and k denote the momenta of the outgoing on shell quark and photon respectively. 
To get the photon production cross section, we will integrate over the momentum of the produced quark 
at the end. The amplitude is given by 

{qiqh{k)out\q{p)m) = (Ooutloout (A:)6out(g)&L(p)|0in) (274) 
which, using the LSZ formalism and the definition of the quark propagator, can be written as 

{q{q)l{kUM{p)in)=ej d^xd'^ydhe'^'-^+^-'-P-y^uiq) 

{i 0^ -m)Sp{z,x)iSp{x,y){i 0y -m)u{p) (275) 

where 5*^ is the quark propagator in the background of the classical color field given in ()262|) and e^{k) 
is the polarization vector of the produced photon. Using the decomposition of the quark propagator as 
a free part and an interacting part as in ()263p . the amplitude can be written as 



{q{qh{k)out\q{p)m) = -eu{q) 



(27r)V(A; + g- 
+ T{q,p-k) Si {p-k)i + i Si {q + k)T{q + k, p) 

u{p) (276) 



+ / ^,Tiq,l)Slil)^Slik + l)Tik + l,p) 

where r is defined in ()265p and Sp is the free quark propagator. The various terms in ()276p are shown 
in Fig. 1331 where the thick wavy lines attached to the full circle signify multiple scattering of the quark 
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Figure 43: Real photon production in quark-nucleus scattering. 

from the nucleus and photon is shown by a dashed line. The first term corresponds to the case when the 
incoming quark does not scatter on the target and radiates a photon. This is kinematically impossible 
so that this term is identically zero (as can be seen by noticing that the delta function has no support). 
The second term corresponds to the case when the incoming quark scatters from the nucleus and then 
radiates a photon while the third term corresponds to the case when the incoming quark radiates a 
photon and then scatters from the target nucleus. The last term corresponds to the incoming quark 
scattering from the target, radiating a photon and then scattering from the target again. In the high 
energy (eikonal) limit considered here, the last term is also zero (or more precisely, suppressed by 
powers of center of mass energy) and is therefore dropped (for a proof of this, we refer the reader to 
|93j). Extracting the overall delta function, the amplitude is 



M.{q,k:p) = —ieu{q) 



u{p) 



{p — kY — rn? {q + ky — m? 
d^x e*(£+^-£)-^ (y{x) - l) . (277) 

To get the cross section, we need to square the amplitude, include the phase space and flux factors and 
average (sum) over the initial (final) state degrees of freedom. We also set the transverse momentum 
of the incoming quark and the quark mass to zero for simplicity. We get the differential cross section 
for production of a photon with transverse momentum kx and a quark with transverse momentum qt 

da"^-'^^'' _ [l + (l-z)^](g + fc)V|,,,efe£ + fe) . . 

d?h(Pk(Pqdz {2Tifkl z [q + k-k/zf ^ ' 

where z = k" /p~ and [1 + (1 — z^j/z is the standard leading order photon splitting function and the 
dipole cross section is defined in ()271|) . To get the cross section for a hadron + photon cross section, one 
would need to convolute the above with the quark-hadron fragmentation function which can then be 
used to investigate the correlation between the produced hadron and photon. It is also remarkable that 
the cross section is written in terms of the same degrees of freedom as in the valence quark production 
cross section. 

To get the single inclusive photon production cross section, we integrate over the transverse mo- 
mentum of the produced quark qx- A quick inspection of Eq. ()278|) shows that the integration over qx 
is divergent. This happens when the transverse momenta of the photon and the final state quark are 
parallel, g = (1 — z)k/z. This is nothing but the standard collinear divergence in pQCD calculation of 
the so called fragmentation photons (as opposed to direct photons), which occurs when the photon and 
produced quark are parallel. This contribution is formally of order of as compared to direct photon 
contribution. However, the large collinear log can be identified as being of order l/us and therefore it 
is customary in pQCD to include fragmentation photons at the same order of calculation as in direct 
photons [HZj. This collinear divergence is not affected by the multiple scattering of the quark on the 
nucleus. 

To do the qx integration j9Bj, we write the dipole cross section back in the coordinate space and 
shift the quark transverse momentum q ^ q + k{l — z)/ z which leads to the following integral 



r ^ e^l-- ^ vr In 5/A|^^ (279) 

Qt' 
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where s is the center of mass energy of the quark-nucleus system. This logarithm is typically written as 
sum of two pieces; In Q'^/Aqi^^ and In s/Q^ where Q is the factorization scale, usually taken to be the 
transverse momentum of the produced photon. With leading log accuracy, only the In Q'^/Aq^j^ piece 
is kept inZ] and absorbed into the quark-photon fragmentation function D^/q{z,Q^) which satisfies 
the DGLAP evolution equation, evolving with Q^. Identifying the quark-photon splitting function 
P^/^ = ^"'"^ z"^'* ' write the photon production cross section as 

f^jj^^-^T^ 1 1 

MkTz = J^2^D,/,iz,kr)al^,Mk/z) (280) 

where the Leading Order (LO) photon fragmentation function is defined as 

D^/,iz,Q^) = P^/g InQVM (281) 

It is straightforward to include quantum evolution of the target wave function at high energies in the 
photon production cross section here. Again, the dipole cross section is the universal object present in all 
single inclusive particle production cross sections in proton-nucleus collisions. By solving the JIMWLK 
equations for the fundamental dipole cross section and using the solution in the above cross section, 
one incorporates small x evolution of the target wave function with Leading Log (of 1/x) accuracy. 

We note that the process considered here is the dominant process at high energy and forward rapidity. 
At fixed energy, as one goes towards mid rapidity, one probes smaller x in the projectile proton where 
gluons become as (or more) prominent as valence quarks. One then can have a gluon split into a quark 
anti-quark pair, one of which would then radiate the produced photon. This process is parametrically 
suppressed by a^. However, this suppression may be partially compensated by the dominance of the 
gluon distribution function in the proton in mid rapidity and may be numerically important. Also, as 
one goes to very high pt at fixed center of mass energy, the valence degrees of freedom in the target 
wave function will become important. This is beyond the realm of application of CGC at its present 
form and would require extending CGC to high x region. 



3.4.2 Dilepton Production 

We now consider production of dileptons in proton- nucleus collisions [HH EHl ISHl ESI HUUj • This follows 
simply from our derivation of the photon cross section. We just need to allow the photon to be off shell 
and include the effect of the massive photon splitting into a dilepton pair. We again start with the 
production amplitude for this process 

q{p) +A-^ q{q) + + /"(^a) + X (282) 

where /ci, ^2 are the momenta of the two leptons. The amplitude can be written as 

(0out|&out(g)&L(p)cout(/i;2)rfout(/i;i)|0in) = j d^xd^d^zi Z2 

(Oout|T^(x)?(2/)^(zi)^(z2)|Oin)(^ +m){i Id,^ \m)v{kx)u{v) (283) 

where and \Ef are the quark and lepton fields while u and w [v) denote quark and lepton (anti lepton) 
spinors. The scattering of the quark on the target is identical to the photon production case considered 
earlier. The only difference with the photon production case is that we need to replace the photon 
polarization vectors by the virtual photon propagators which amounts to the following substitution in 
the squared amplitude 

^ lek^ — le 
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with k = ki + k2 the virtual photon momentum and the leptonic tensor LP"{ki,k2) is the imaginary 
part of the one loop leptonic contribution to the photon self energy. Furthermore, we are interested in 
the dilepton pair rather than the individual leptons. Therefore, the leptonic tensor can be written as 



Extracting the overall delta function again, the amplitude becomes 



(285) 



M{p\qk) 



-iequ{q) 



(p — ky — m? 



(g + kY — Tn? 



u{p) 



(286) 



The two terms of the amplitude are shown in Fig. where again, as in the real photon production 
case, the virtual photon can be radiated either before or after the multiple scattering from the target. 
Note that this amplitude looks exactly as the amplitude for photon production in ()277j) . The difference 





Figure 44: Dilepton production in quark-nucleus scattering. 

is that, here the photon 4-momentum k is off shell {k"^ ^ 0). Squaring the amplitude, including the 
flux factors and averaging (summing) over the initial (final) state degrees of freedom, the cross section 
is given by 



i+i- X 



d%d^kd\nM^ dz 
'l + (l~zY 



37r 



(2vr) 



i'^dipole\'± 

Z Jjrp 



(q + k) 



-z(l-z) M' 



+ M2(l - z)][{k -zl_y + M2(l - z)] 
1 1 



[kl + M^{l-z)] [{k-zl_y + M^{l-z)] 



(287) 



where is the dilepton invariant mass squared, i = q + k and we have ignored quark masses. As 
before, z is the fraction of the incoming quark light cone energy carried away by the photon. 

Since the dilepton production cross section is proportional to the electromagnetic coupling constant 
almi the production rate is quite small, specially at higher transverse momenta. One way of getting 
around this problem is to consider transverse momentum integrated dilepton production cross section. 
One can then investigate the behavior of this cross section with varying dilepton mass. To do this, 
we write the dipole cross section in transverse coordinate space, and integrate over the transverse 
momentum of the dilepton pair kx in ()287|1 . We get 



' d^dM^dz 



l-z 



37r2 



[1 + (1-^)1 



-Mrr] + 2{1 ~ z)K^[ 



'-Mm 



(288) 
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To relate this to proton-nucleus scattering, we need to convolute ()288p with the quark (and anti- 
quark) distributions q{x,M'^) in a deuteron (proton). As shown in ^98. .99(177^] . this can be written in 
terms of the proton structure function F2 



daP 



A^l+l- X 



d%dM^dxF 



Gvr^ Xg + x„ 



1 r 1 
dz I dr^ — 



Fi{xq/z)aip^i^{xg,b, tt) 



[1 + (1 - zf] Kll^^^Mrr] + 2(1 - z) K',[^ 



'-Mn 



z 



(289) 



where 



Xn 



X a 




+ Xf 



Xf 



(290) 



with XF = ^[ey - e-y] and 

PP ^Y.^[qf{x,M^) + qf{x,M^)] 
f 

is the proton structure function. 

Inclusion of small x evolution of the target wave function is again straightforward. As in the case 
of photons and elastic quark and gluon production, the dilepton production cross section depends on 
the dipole cross section which satisfies the JIMWLK or BK evolution equations. Using the solution 
of JIMWLK or BK for the dipole cross section cr^^pQi^ in (j287j) and (j288j) would automatically incor- 
porate small X evolution of the target nucleus wave function. Finally, it is worth noting that there 
is a deep relation between dilepton production in proton-nucleus collisions and in DIS of virtual pho- 
tons on a hadron/nucleus due to the approximate crossing symmetry of the amplitudes in the Color 
Glass Condensate [101 1175j (the crossing symmetry of the amplitude is approximate due to the eikonal 
approximations made). 



3.5 Hadronic Two-Particle Correlations 

In this Section we will discuss two-particle inclusive production in pA collisions. We will review produc- 
tion of two gluons, of a valence quark and a gluon at forward rapidity and of a qq pair at mid-rapidity. 



3.5.1 Two-Gluon Production 

Here we are interested in calculating an inclusive production cross section for two gluons with transverse 
momenta k^, ^"^^ rapidities y2 in pA collisions or in deep inelastic scattering. A diagram contribut- 
ing to two-gluon production in DIS is shown in Fig. |^ Again the target nucleus has rapidity and 
the projectile (proton) has rapidity Y. The gluons are produced at central rapidities ^ yi,y2 ^ Y- 
For simplicity we assume that ?/2 ^ Hi- generalization of the results we will present below to yi ~ y2 
appears to be conceptually straightforward but extremely cumbersome technically. 

Without going through details of the calculation, which was worked out in |103j . we will just quote 
the answer for two-gluon production cross section at mid-rapidity for gg-nucleus scattering including 
small-x evolution. After a tedious calculation one obtains 

(^(j9<?A qqGiG2X 

d^ki dyi d'^k2 dy2 



80 



nucleus riQ ^ j Q ir 

nucleons 



Figure 45: Two-gluon production in DIS on a nucleus including multiple rescatterings. 

X (fxi (fxi (fx2 d^X2 s{xi,xi,ki, Ui) s{x2, 212, 2/2) + ^i(2i0) 3io5 Y'j 2ii5 7/2) c^^a^i c^^a;]; 



a. 



{2n) 



(fx2(fx2'e -2-2 



"T" ~ "T" 1 ■ ( "T^ ~ I M{^2,X2',xi,y2;ki,yi) S'(x2, X2/, 2/2) 
,3;2i \a;2'i ^2'i/ 



X21 X^j 



2'1 



^2'1 _ ^2'1 \ ^21 
^•^2'1 *^2'1' 



^21 



X M(Xi,X2/,a>2/2;^l,2/l)'S'(Xi,X2/,?/2) + ^ ■ ^M(Xi,Xi,Xi,?/2;^l,l/l) + (1^1) 



(291) 



X21 X^ii 

Eq. ()29H1 gives us the two-gluon inclusive production cross section for a scattering of a qq dipole on a 
nucleus. (Generalization to "pA is discussed in |l()3j .) The transverse coordinates of the quark and anti- 
quark in the dipole are Xq and Xq, with the impact parameter B_ = (xq + Xo)/2. ni(xo, Xq, Xi, xj, 2/2) 
is given by Eq. 1)2211) : to convert it to the notations of Eq. ()22H) we write 



ni{xo,x^,Y;xi,xi,y2) ^ rii(xoo, x^, 6 = -(xq + Xq - Xi - xj), F - ?/2)- (292) 

^2(210? 2i05 ^! !^ii^2Li, Z/i, 2i25 2i25 2/2) is the number of pairs of dipoles formed by quark and anti-quark at x^ 
and xj with rapidity yi and by quark and anti-quark at X2 and X2 with rapidity y2 in the original dipole 
Xqq having rapidity Y. n2 obeys the following equation [201 



ri2(xo,Xo, y ;Xi,Xj;,|/i,X2,X2,|/2) = 7^ dye ^ p '"^ ^> d 

In JmSiX{yi,y2} Jp 



2 ^00 

P X30Xgg 



X [ni(xo,X3,2/;xi,Xj;,?/i)ni(x3,Xo,?/;x2,X2,i/2) + ni(xo, X3, X2, X2, 2/2) ni(x3, Xq, Xi, xj, ?/i) + 

+ ^^2(^0, ^3,1/; ^1,^1,1/1,^2, ^25 1/2) + n2(x3,Xo,y;Xi,Xi,?/i,X2,X2,?/2)] 

which is linear and can be solved after one finds rti from Eq. ()221)1 . 

The quantity s in Eq. ()291|) is defined in terms of gluon dipole amplitudes as 



(293) 



■s(xi, 212, ^, y) 



a. 



d^Zi d^Z2 e 



i+j 



Z_Y X.2 —"l —j 



I /y 12 I .-y /y I 2 

\Al X_i\ \±.2 X j I 



X 



Ng {zi,Xj, y) + Ng {z2,Xi, y) - Ng (^i, ^2, v) - Ng (x^, x^, ; 



(294) 
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Figure 46: Diagrams describing one step of evolution for M. 



and is an integral component of single inclusive production cross section from Eq. p99|) . S in Eq. ()29H) 
is the S'-matrix of a quark dipole scattering on the target nucleus, related to the quark dipole amplitude 

N by 

S{xo,x^,Y) = l-N{xo,x^,y), (295) 

where A^, in turn, can be found from the nonlinear evolution equation (jl35p . Note that the gluon dipole 
amplitude Nq is defined by Eq. (j226|) . 

Finally, the quantity M{x2,X2',xi,y2', ki,yi) in Eq- f)29H) has the physical meaning of production 
amplitude of a gluon with transverse momentum and rapidity yi in an off-forward scattering of dipole 
21 into dipole 2'1 on the target nucleus. The evolution equation for M is |IT)H] 

M{x2,X2',x^,Y;ki,yi) = e v " J d{x2,X2f,x^,ki,yi) + 



a 



H — - I d X4, dye 

2tT J Jyi 



Y 
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X 



2~ 



40' 



X' 
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^2 
J, 42 X 



2i42 2i42' 



M{x^,X2,,XjQ,y;k^,yi) [1 - A^(x2, X4, y)]- 

, 3^42' 2:42 / 
/ 2i42 



- ^ ) M{x2,x^,XQ,y;ki,yi) [1 - N{x2',X4,y)] \. 



X2 
, 42 



X2 
42' 



X2 
^ 42' 



^40 



(296) 



It includes the (nonlinear) evolution between the incoming dipole and the produced gluon, which is 
shown in Fig. using the same notation as in Fig. EHl The initial condition for the evolution in 
Eq. ()296|) given by d{x2,X2i,x^,ki,yi) has to be calculated separately by solving a linear evolution 
equation (see Eqs. (27) and (28) in |l()Hj ). 

Eq. ()29H) is represented graphically in Fig. |17| in terms of traditional Feynman diagrams. The first 
term in Eq. ()291|) corresponds to splitting of the original linear evolution in two, which is described 
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Figure 47: Feynman diagrams corresponding to double gluon production cross section given by 
Eq. (j291|) . Emitted gluons are denoted by crosses. 

by 712- Then each of the two ladders independently produces a gluon with all the possible splittings 
happening afterwards. This is illustrated in Fig. IT7l\. The second term in Eq. (j291|) corresponds to 
nonlinear evolution successively producing both gluons, after which all possible splittings are allowed, 
as shown in Figs. 147b and ITflC. where we have divided the nonlinear evolution into the linear (Fig. 
147b) and nonlinear (Fig. |37p) parts. The linear evolution in this second term in Eq. ()291|) is given 
by rii and by M from the linear part of Eq. (j296|) . This linear evolution leads to production of both 
gluons #2 and #1 and is illustrated in Fig. 147b. The initial conditions for Eq. (|296|) are nonlinear. 
They include ladder splittings and are pictured by the fan diagram in the lower part of Fig. 147b. One 
should note, however, that Eq. ()296|) . while being linear in M, has extra factors of 1 — on its right 
hand side. That means that evolution of M includes ladder splittings between the gluons #2 and #1, 
one of which is shown in Fig. WTC. There the evolution leading to creation of gluon #2 is still linear 
since it is still given by rii in the second term on the right hand side of Eq. (j29H) . However, since the 
evolution in the rapidity interval between the emitted gluons (evolution of M) is nonlinear, splittings 
are allowed between the gluons #2 and #1, as depicted in Fig. HTp. 

The diagrams A and B in Fig. |37| are the same as would have been expected from AGK cutting 
rules |17Uj (see also [S^ for similar correspondence between the dipole model results and AGK rules 
expectations). However, the diagram C in Fig. |37l while being included in Eq. ()291|) . is prohibited 
by AGK cutting rules. Therefore, we seem to observe a direct violation of the AGK rules in QCD. 
Since AGK rules have never been proven for QCD, one should not be too surprised that they do not 
work here. It is interesting to note that AGK violation sets in at the level of the 2-gluon production: 
single gluon inclusive production cross section calculated in [SI] adheres to AGK rules and so does the 
diffractive DIS cross section calculated in |131j . 

The violation of AGK cutting rules in Eq. (j29ip is due to non-linear terms in Eq. (|296|) . which 
are in turn due to late time (after the interaction) gluon emissions at light cone times r > 0. These 
terms were not important for the calculation of the total cross section in the dipole model 
there they were found to cancel (221 ■ Thus if one would try to construct an analogy between the fan 
diagrams [T] and dipole calculations P3] based on the correspondence of total cross sections, one would 
omit such terms. Since the fan diagrams seem to adhere to AGK rules, this omission would lead to the 
erroneous conclusion that AGK rules should work for the two-gluon production cross section. However, 
as we have seen above, these late time emissions are important for single and double inclusive 
gluon production, violating the AGK rules for the latter. What appears to fail here is the one-to-one 
correspondence between the fan diagrams and dipole calculations. 

Another difference between the result ()291|) and the direct application of AGK rules to calculating 
inclusive cross section done in [90j is that nonlinear splittings may start exactly at the point in rapidity 
when the softer of the produced gluons is emitted in the diagram B or exactly at the point of emission 
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of both gluons #1 and #2 in the diagram A. Similar discrepancy was already observed when comparing 
the single gluon inclusive cross section calculated in [34, to the results of ^90j. 

It is interesting to note that Eq. ()29H) can not be cast in the form adhering to fcy-factorization, 
as it was done for Eq. ()227p transforming it into Eq. ()228p |lU3j . This indicated a breakdown of kx- 
factorization for two-particle production. Below we will observe the same failure of fc^-- factorization for 
the case of quark production ll(J4j . 

Eq. ()29H1 gives a complete result for two-gluon production cross section in saturation/Color Glass 
Condensate formalism. It can be used to describe two-particle correlations in dAu collisions at RHIC 
and in pA collisions at LHC. Unfortunately Eq. 1)2911) is rather complicated. In order to make Eq. ()29H) 
easier to implement, it is highly desirable to find some way of simplifying it. At this time we are not 
aware of any simplification of Eq. ()29H) in the general case. Nevertheless, in certain kinematic regimes 
Eq. ()29H) may be simplified. For instance, if the center of mass energy of the collision is not too high 
or if the transverse momenta of the produced gluons are sufficiently large 1^2! ^Qs), the nonlinear 
saturation effects, such as ladder splittings, could be neglected. This implies that the diagrams in Figs, 
im A and C are small with the linear part of the diagram B dominating the cross section. This is the 
well-known leading twist result |176j (see also |102p . 

In jlUlj , using /cT^-factorization formula, Kharzeev, Levin and McLerran suggested that back-to-back 
particle correlations would be suppressed in p{d)A collisions, when the rapidity interval between the 
jets is large, e.g., when one of the jets (1/2) is going in the proton fragmentation region and the other 
one is at mid-rapidity at RHIC (the so-called Mueller-Navelet jets |177j ). This prediction will be 
discussed in more detail in Sect. 14.2.21 While preliminary data from STAR collaboration appears to 
confirm this prediction |178j . it would be important and interesting to generalize the approach of |l()lj 
beyond the /cT-factorization formula, which is valid only at high kj- ^ Qg. Eq. 1)2911) provides us with 
an opportunity to carry out such a generalization. 

3.5.2 Gluon— Valence Quark Production 

Here we consider production of a high x (valence) quark and a gluon from scattering of a quark in the 
proton wave function on the target nucleus, treated as a Color Glass Condensate. This is expected to be 
the dominant process in the projectile proton fragmentation region where one probes the high x degrees 
of the freedom in the proton and the small x degrees of freedom in the target nucleus. This process is 
the as correction to the scattering of the high x quark on the target considered in the previous sections. 

Production of a valence quark and a gluon was considered in [103j. Here we go through the main 
part of the derivation and refer the reader to jlU3j for more detail. The starting point is the amplitude 
for scattering of a incoming massless on-shell quark on the target, producing an on-shell quark and 
gluon shown in 



where Sp, G^p are the quark and gluon propagators in the classical field background and Df'^ is the 
"inverse" propagator (which amputates the external gluon line), defined by 



q{p)A^ q{q)g{k)X 



(297) 



which is given by 




(298) 




(299) 
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and G^p is the free gluon propagator and the quark and gluon hnes with a thick circle represent the 



propagators in the background field as illustrated in Fig. 021 The quark and gluon propagators in the 
classical background field are known |152| I179j . Here we write them as an interacting part and a free 
part as before 



SF{q,p) ^ {27r)'S\p-q)SUp) + SUq)rf{q,p)SUp) 
G'^''{q,p) ^ {27rY6\p - q) G'^^{p) + Gf{q) T,{q,p) G'^'^ip) 

where the free propagators are 



sUp) 



i— and 
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(300) 



(301) 



and is the light cone gauge vector so that t] ■ A = A~ = defines the gauge we choose to use. One 
of the advantages of using this gauge is that the interaction part of the gluon propagator in Eq. (|300p . 
denoted here by Tg{q,p), is diagonal in Lorentz indices (proportional to Qi^^^) which is why we can write 
it as in Eq. ()300|) . The interacting parts of the quark and gluon propagator Tf and Tg are defined as 



y{(l,P) = {2n)6{p -q )Y 
'Tgiq^p) = 2p^ {2tx)6{p' - q' 



1] 



(302) 
(303) 



where V is the multiple scattering matrix in the fundamental representation defined in ()127|1 and U is 
the analogous matrix in the adjoint representation. The scattering amplitude then can be written as 



M{q, A, k- p) = ef{k) [Mf + + + M^] 
where e'^^\k) is the polarization vector of the produced gluon and 



(304) 



Ml- 
My. 



-g u{q)ie SUq + k) Tf{q + k, p) u{p) 
-9 Tfiq, p-k)S%{p-k)^ u{p) 
-g u{q) 7.t' r'^ik, p - q) u{p) G',^{p - q) e,{k) 



(305) 
(306) 
(307) 

(308) 

The different terms in Eqs. (|305|) - (j308p are depicted in Fig.ESlwhere the solid line denotes the incoming 
quark, the thick wavy line attached to a full circle denotes multiple scattering from the target and the 
dashed line represents the radiated gluon. 



M, = -g dH uiq)Tfiq, p - 1)SUp - 07 l)uip)Gril)e,ik). 



q(p) 



g(k) . 



q(q) 






Ml 



M- 



Figure 48: Production of a quark and a gluon in quark-nucleus scattering. 

Ml corresponds to radiation of a gluon before the multiple scattering before the quark and the 
target while M2 corresponds to the incoming quark scattering on the target and then radiating a gluon. 
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Ms depicts the case when the radiated gluon multiply scatters from the target, while in M4 both the 
radiated gluon and the final state quark multiply scatter from the target. To get the production cross 
section, we need to square the amplitude. This is straightforward but tedious. First, one can do the l~ 
integration in M4 using the delta function 6{l^ — k~) and the integration via contour methods using 
the p — I pole which sets /~ = k~ and = —l\/2q^. 

Contribution of Mi and M2 terms is identical (up to color matrices) to the photon production case 
except for the extra contribution coming from the pieces of the gluon propagator proportional to the 
gauge vector rjp which is not present in the case of the inclusive photon production. This warrants some 
clarification. In single photon production, photons are emitted from the incoming quark and can not 
be emitted from the target (CGC) which, to leading order in as is comprised of gluons only. Radiation 
of a photon from the target is suppressed for two reasons; first it is a higher order in as correction since 
one would need a gluon in the target wave function to split into a quark anti-quark pair which would 
then radiate a photon and second, it is suppressed in the soft photon limit (it is proportional to z where 
z is the fraction of the parent quark the radiated photon carries away). Therefore, in the case of photon 
production, only the g^u part of the propagator contributes. In case of gluon production, one needs 
to consider radiation of gluons from the target as well as from the projectile quark. One can avoid 
these contributions by working in the light cone gauge as shown by |j32j. This however necessitates 
inclusion of the contribution of the extra terms, proportional to the gauge vector, in the propagator. 
The contribution of Mi and M2 is therefore 



(Mf + Ml^) 



1] ■ k 



(M^ + M^) = 16p-p- 
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(309) 



where Tr denotes trace of color matrices. A tedious but straightforward calculation of the rest of the 
diagrams leads to 
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[U'^ik- L) - 5'"'i27rY5\k- l)]Trt\V{q + L) - {27ry5\q + l_)]t'' 



\mIm, 




Trt'[V\q + l) - (27r)V(£ + i)r[V(g + i) - (27r) V(g + i)] 



\MlM2 



Trf[V\q + V) - {2'Kf5\q + V)][V{q + k) - {2'Kf5\q + k)]t' 




[U''\k-V)-5''\2'Kf5\k-V)]TrtH''[V\q + l) - {2'Kf5\q + Vi\ 
[V{q+T)-{2'nf6\q + l_)] 



(310) 



where for any interference term, there is an equal conjugate contribution. To get the cross section, one 
needs to average (sum) over the initial (final) state degrees of freedom and include the phase space and 
flux factors. The invariant cross section for production of a quark with transverse momentum qx and 
energy fraction z and a gluon with transverse momentum kx (and energy fraction 1 — z) is given by 



where |Mp is defined in ()304|) and given in ()309I310|) . This cross section looks quite more complicated 
than single inclusive quark or gluon production. Nevertheless, the degrees of freedom are again Wilson 
lines in fundamental or adjoint representation as before. Here, one has correlators of more than two 
Wilson lines (which were not present in single inclusive production of quarks and gluons). To do a 
quantitative study of two hadron production in proton-nucleus cross section, one needs to convolute this 
cross section with the distribution function of the quark in a proton and the two hadron fragmentation 
function. 

One can use the two particle cross section derived here (after the convolution with distribution and 
fragmentation functions) to investigate two hadron correlations in proton-nucleus collisions at RHIC 
and LHC The expectation is that due to classical multiple scattering (no quantum evolution), the 
correlation function will become wider. Inclusion of quantum corrections can be achieved by solving 
the JIMWLK equations for the higher point functions of the Wilson lines and using the solution in the 
above cross section. This should be valid as long as the produced quark and gluon have the same (or 
similar) rapidities. One expects that inclusion of quantum evolution would lead to reduced correlation 
between the produced hadrons. The effects of quantum evolution on the two hadron correlation should 
be most dramatic when the produced hadrons are widely separated in rapidity, as considered in the 
previous section. 

3.5.3 Quark— Anti-quark Pair Production 

We finally consider production of a quark anti-quark pair in proton-nucleus collisions. This has been 
investigated in detail in (TUS CHHl IMl IM 1103 • Here, we follow the derivation of the [10311001 since 
it is most closely related to the Color Glass Condensate formalism. For technical details and more 
references, we refer the reader to |lU8j . The starting point is the classical background field of a single 
nucleus in terms of the nuclear color charge density given by 




(311) 




(312) 
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in the covariant gauge and its Fourier transform in momentum space 

A'^^iq) = 27196'^- 5iq+)^ . (313) 

Qt 

Note that ()313|) is hnear in the color charge density even though it does describe the classical field 
created by a dense nucleus. This is a property of the covariant gauge chosen here and is not true in a 
general gauge. 

Since proton is treated as a dilute object while the nucleus taken to be a Color Glass Condensate, 
one needs to keep only the linear term in the color charge density of the proton while in the nucleus, 
the color charge density is kept to all orders. The general solution of the classical equations of motion 
for a proton-nucleus collision was found in jlU8t and is given by 

A'^iq) = A>^^^iq) + 5^-A-^giq), (314) 

where A'^^^ is given by 
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and the singular term is given by 
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where U and U (note the factor of 1/2 in the definition of U) are matrices in the adjoint representation 
and 
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and T"- are the generators of the adjoint representation of SU{N) and ki is the momentum transfer 
from the proton and = g — fci is the momentum transfer from the nucleus. The coefficients C^{q,ki) 
and C^j.^g{q) are defined by 



C^{q,k,T)^ 



h.2 1.2 _ 2 

,C-iq,k,r)^ '^ 



q +te 



-2k{ 



(318) 



2 

and C^.^^g{q) and C^{q) are related through C!^,^^g{q) = C^{q) To get the amplitude, one again 

uses the quark propagator in the background field of the nucleus which includes multiple scattering 
from the target nucleus and includes multiple scatterings from the quark, anti-quark and gluon lines. 
The final amplitude is given by 

d'^ki d^k pp ,a{k^) f^2^^2y^ik-x^i(2\^k-k^yy 
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where V is in the fundamental representation and p and q are the momenta of the quark and anti- 
quark and m denotes the quark mass. Note that the Wilson line U (with the 1/2 in the exponent) has 
dropped out of the expression for the amplitude due to cancellation between the regular and singular 
contributions. To get the cross section, we need to square this amplitude, sum (average) over the final 
(initial) state degrees of freedom and color average over the sources which is tedious but straightforward. 
Here we just give the final result |lU9j 



da 
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d kid k2 — ~ ~ 
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where 
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and Tgq, T'^, Tg are related to Eqs. (|321|) via 



Tqq — Tqq{kl^k) j T q — Tqq{k_l^kt ) ) ^3 — Tg{k_i] 



and yj's are the Fourier transforms of two, three and four point functions of Wilson lines jl09j . 
example. 
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There are a few remarks in order here. First, the quark anti-quark production cross section is again 
expressed in terms of the Wilson lines just like the other cross sections in Color Glass Condensate 
approach to high energy proton-nucleus collisions. Furthermore, the standard /cT^-factorization does 
not hold in general, at least in its most common form, even though one does recover the /cT^-factorized 
form of the cross section in the leading twist region. For a quantitative study of the breakdown of /ct 
factorization in this context, we refer the reader to |ll()j . We also note that one can integrate over the 
momenta of the final state anti-quark (or quark) in order to get the cross section for single inclusive 
quark production. Finally, to include quantum corrections of the form Inl/x in the cross section 
above one needs to solve the JIMWLK equation for the n point function of Wilson line and use them 
in the cross section in ()320|) . This may then be used to quantitatively investigate production of heavy 
(or light) quark anti-quark pairs in proton-nucleus collisions at RHIC and LHC. 



4 Results from dAu Collisions at RHIC 

Now we are going to review the experimental data obtained from d + Au experiments at RHIC. We 
will show that at mid-rapidity d + Au collisions provide a control experiment giving us evidence of jet 
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Figure 49: Nuclear modification factor for charged hadrons and neutral pions produced in d + Au and 
Au + Au collisions at mid-rapidity, as reported by (clockwise from left) PHENIX ^51j, STAR and 
BRAHMS mH collaborations. 

quenching due to strong final state interactions indicating formation of quark-gluon plasma. We will then 
proceed by reviewing the results from forward rapidity in d + Au collisions: they confirm the qualitative 
expectations of saturation/ Color Glass Condensate physics, summarized in Fig. ^2 Therefore we argue 
that RHIC experiments have seen the first experimental confirmation of Color Glass Condensate in 
nuclear collisions. 

4.1 Mid-Rapidity 
4.1.1 Cronin Effect in dAu 

The experimental program at RHIC involved d+Au collisions at the center of mass energy y/s = 200 GeV 
per nucleon. While the above analysis of particle production was, strictly speaking, designed for pA 
and DIS, we believe that its results can be applied to ci + Au collisions as well, d + Au scattering 
also involves a scattering of a small projectile on a large nucleus, similar to DIS and pA. The only 
assumption we have made in our above analysis was that the saturation scale of the projectile is much 
smaller than transverse momentum of produced particles, which allowed us to neglect saturation effects 
in the projectile's wave function. The fact that deuteron has the atomic number of A = 2 increases its 
saturation scale compared to that of a proton by a small factor of A^^^ ^ 1.26. Therefore, the region of 
transverse momentum above the deuteron saturation scale, where our analysis is applicable, would still 
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be very wide in d + Au collisions. It is, therefore, a good approximation to apply our above conclusions 
to d + Au collisions. 

Nuclear modification factor for charged hadrons produced in ci + Au collisions at mid-rapidity was 
measured by all four experimental collaborations at RHIC 13^ llllj . The data from BRAHMS 

[TTT] . PHENIX ini and STAR |53j collaborations is shown in Fig. Hi and the data from PHOBOS 
collaboration [32] is shown in Fig. EDI One can see that i?''^" at mid-rapidity at RHIC exhibits Cronin 
enhancement jHl] , in agreement with theoretical expectations of multiple rescattering models, shown in 
Fig. El ^ Ea Ea ES .88» ^ El El • 




Figure 50: Nuclear modification factor for charged hadrons produced in d + Au collisions reported by 
PHOBOS collaboration in pseudo-rapidity range of 0.2 < rj < 1.4 for four bins of centrality. 

The PHOBOS data, presented in Fig. EOl allows one to study i?'^^" as a function of centrality. 
Analyzing Eq. ()219p . along with its asymptotic high-p^ limit of Eq. ()217p . we see that high-p^ R^^'^ 
in the quasi-classical multiple rescattering approach is an increasing function of A, or, equivalently, 
centrality. These theoretical expectations are confirmed by the data in Fig. EOl where i?'^^" is clearly 
an increasing function of the collision centrality. 

An intriguing property of the PHOBOS data from Fig. EHl |12] is that the height of the Cronin 
maximum in it is much lower than in the data from BRAHMS, PHENIX and STAR collaborations 
shown in Fig. Eni[llH lKT | IHH]. The origin of this disagreement appear to be in the fact that PHOBOS was 
collecting data in the pseudo-rapidity interval of 0.2 < rj < 1.4, whereas the other three collaborations 
presented the data at r/ = 0. This decrease of Cronin maximum with increasing rapidity can already be 
interpreted as a precursor of the suppression predicted by saturation physics shown above in Fig. 1^ 
However, we will postpone this discussion till Sect. I4.2.1l where the data at more forward rapidities will 
be presented. 

Another interesting feature of the d + Au data at mid-rapidity is that other observables, such as 
integrated charged hadron multiplicity as a function of rapidity j3U ES], appear to be in agreement 
with saturation/CGC models. In Fig. |^ we show a prediction of saturation model from [031 
total charged particle multiplicity dNch/drj, compared to the data from [33] • The model in jHS] uses kj- 
factorization formula for gluon production from Eq. ()231|) integrated over all fcy with the unintegrated 
gluon distribution on the nucleus "frozen" at a constant at transverse momenta below Qs and falling 
off as ~ l//c|. for kr > Qs (for details see [Ml)- As one can see the agreement between the prediction 
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Figure 51: Predictions of saturation model from |£5J (solid and dashed lines) for the total charged 
particle multiphcity dNch/drj in d + Au collisions compared to the experimental data from [S3] as a 
function of pseudo-rapidity. 

from [HH] and the data is very good. 

4.1.2 dAu as a Control Experiment for Quark-Gluon Plasma Production 

dAu collisions can be used as a control experiment for formation of quark-gluon plasma (QGP) in 
AuAu collisions at RHIC. As was originally suggested in |115| llltil 11171 I^THl I119j . creation of a dense 
hot medium, such as QGP, in high energy heavy ion collisions should lead to a depletion of high-p-r 
particles produced in the collisions. This phenomenon is known as jet quenching |118j . The underlying 
physics is rather straightforward: high-p^ particles produced in the medium would loose energy in 
interactions with the medium. As a result there would be less high-p^ particles produced in the heavy 
ion collisions than one would naively expect by scaling the number of high-p^ particles produced in a 
pp collisions by the number of elementary nucleon-nucleon collisions Ncou- 

One usually measures this suppression by analyzing the nuclear modification factor for AA collisions, 
which is defined as 



Jet quenching due to energy loss of particles in the medium would result in R'^^ being suppressed. The 
calculations of such suppression were carried out in |116t 11171 1118j . where both multiple rescatterings 
in the medium and medium- induced particle emissions were taken into account. 

The data on i?^^ collected at RHIC was reported in |112[ I113| I114[ llllj . Here those data are shown 
in Fig. EHl along with the nuclear modification factor for d + Au collisions. As on can see from Fig. 
jet quenching was observed in Au + Au collisions at RHIC with < 1 for all measured pt- 

For jet quenching to become a signal of creation of quark-gluon plasma one has to prove that the 
observed suppression is really due to energy loss of jets in the produced matter, which would be a final 
state effect. In principle, jet quenching can be accounted for by the small-x evolution effects in the 
nuclear wave functions, similar to those described for pA collisions in Sect. 13.21 leading to suppression 
shown in Fig. |^ Such suppression would be the initial state effect. In fact, it has even been suggested 
in |61J that the suppression observed in jll21 11131 ITTT] could be due to initial state effects induced 
by small-x evolution. The initial state suppression should be observable already in d + Au collisions, 
while the final state energy loss is typical only for AA collisions and should not manifest itself in d + Au. 
In that sense d + Au scattering is a control experiment for QGP formation 
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The data presented in and shown here in Figs. |3H1 and EOl clearly indicate absence of 

suppression in W^^"^. The Cronin enhancement observed in these data is in agreement with the classical 
gluon field dynamics for pA collisions of Sect. 13.11 the same initial state classical fields would lead to 
Cronin-like enhancement of R^^ in AA collisions |18Uj . One can therefore conclude that jet quenching 
observed for AA collisions in |112[ I113t I114[ llllj is due to strong final state interactions leading to 
energy loss of partons in the medium. 




Figure 52: Two-particle correlation functions for pp, d + Au and central Au + Au collisions measured 
by STAR collaboration at mid-rapidity |181j . 



The case for production of hot and dense medium is strengthened by observation of two-particle 
correlations. In pp collisions a high-p-r jet is usually accompanied by another high- or intermediate-pT jet 
balancing its transverse momentum. This is known as hack-to-hack jet correlations. In the presence of 
dense medium absorbing many of the produced jets, such back-to-back correlations should be reduced, 
and possibly wiped out completely. To see whether this takes place, and to eliminate possible initial 
state nuclear wave functions effects, one has to study two-particle correlation functions in Au + Au, 
d -\- Au and pp collisions. 

The data on such correlations was reported by the STAR collaboration in jl81j and is shown here 
in Fig. Fig. |S21 depicts the correlation function (with constant background subtracted) between a 
trigger particle (hadron) with 4 GeV < pxitrig) < 6 GeV and and associated particle with 2 GeV < 
Pt < Pritrig) for pp, d + Au and central Au + Au collisions as a function of the azimuthal angle 
A(f) between the two particles. First of all one should note that back-to-back correlations in Au + Au 
collisions (with A0 = vr) are completely wiped out in Fig. |S21 Secondly, the back-to-back correlations in 
pp and d+Au are, approximately, the same. Therefore one concludes that disappearance of back-to-back 
correlations in Au + Au collisions is due to final state effects. 

Final state interactions have to be quite strong to eliminate more than half of the jets (see Fig. 
and completely wipe out their back-to-back correlations (see Fig.l52|l. Such strong interactions are likely 
to lead to thermalization of the produced dense medium: if partons interact strongly with high-p^ jet 
particles, they would interact strongly with each other reaching thermal equilibrium. It is therefore 
quite plausible that d -\- Au scattering experiments indicate the formation of QGP in Au -\- Au collisions. 
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4.2 Forward Rapidity 



4.2.1 Suppression at all p^: Evidence for the Color Glass Condensate 

The results on nuclear modification factor R'^^'^ at forward rapidity were first presented by BRAHMS 
collaboration in [JE] . The final version of the data published in [2Zj is shown in here in Fig. EHl There 
the nuclear modification factor is plotted as a function of transverse momentum for charged 

hadrons at rapidities = 0, 1 and for negatively charged hadrons at rapidities rj = 2.2, 3.2. The data 
of Fig. ESI demonstrates suppression of hadron production as one goes from central towards forward 
rapidity. We can clearly see that Cronin maximum gradually disappears and suppression sets in as 
rapidity becomes larger: this clearly conGrms the qualitative prediction of Fig. ^2 based on Color Glass 
Condensate/saturation physics [HH ESJ E3] • 
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Figure 53: Nuclear modification factor for charged hadrons at rapidities 77 = 0, 1 and for negatively 
charged hadrons at rapidities 77 = 2.2,3.2 reported by BRAHMS collaboration in j57j . 

BRAHMS data on forward suppression in d + Au from jST] is independently confirmed by other 
RHIC experiments. In Fig. |^ we show the data on hadronic nuclear modification factor R^^^^ reported 
by PHOBOS collaboration in jSHj- While PHOBOS collaboration does no have such a wide rapidity 
acceptance as BRAHMS, one can see the onset of suppression of at high already at pseudo- 

rapidities of ?7 = 1.0 -i- 1.4. 
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Figure 54: Data on nuclear modification factor for hadrons presented by PHOBOS collaboration 

in inH|. 
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Figure 55: Data on R'^^ for hadrons presented by PHENIX collaboration in |59j . 



In Figs. EHl and EHl we show the data presented by PHENIX jSH] and STAR jHO] collaborations. 
Instead of RF"^ from Eq. ()208|) one sometimes is interested in ratio of particle production rates in 
central over peripheral collisions, which is defined by 



(Pkdy 



(central) 



d^kdy 



(peripheral) 



(325) 



As once can see by comparing Eq. ()325|) to Eq. ()208p . R^^ is very similar to i?^"^, with the difference 
of using the peripheral pA collisions instead of pp as a reference in the denominator. The conclusions 
of the above analysis of R^"^ summarized in Fig. |^ would also apply to R^^: we would also expect a 
transition from Cronin enhancement to suppression in R'~^^ with increasing rapidity. 

Fig. Eni plots PHENIX data on R'~^^ for hadrons as a function of px for three different centrality bins 
both for forward rapidity (1.4 < rj < 2.2, denoted by squares) and for backward rapidity (—2.2 < t] < 
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Figure 56: Data on for neutral pions presented by STAR collaboration in jnOI superimposed with 

the BRAHMS collaboration data for negatively charged hadrons from |,57j . 
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—1.4, denoted by circles) jSH]- We see that i?*"^ is suppressed at forward pseudo-rapidities lA < r] < 2.2 
in agreement with Fig. ESI Fig. ESI presents prehminary STAR data on for vr^'s as a function 

of transverse momentum pt at rapidity rj = 4 [HU] superimposed with BRAHMS data on negatively 
charge hadrons at lower rapidities. One can see that suppression gets progressively stronger at forward 
rapidities, in agreement with Fig. ESI leading to a very strong suppression of vr'^'s at ?7 = 4. 

Of course the qualitative confirmation of the CGC-based predictions is a very important experi- 
mental evidence for saturation/CGC physics. However it is important to demonstrate the quantitative 
agreement of BRAHMS data [HZl with the CGC expectations. A saturation/CGC-based model was 
constructed in jl22j . which involved a new parameterization of the dipole- nucleus scattering amplitude 
Ng (see Eq. fl328|) below and discussion around it), with the saturation scale matching that of the 
Golec-Biernat-Wiithoff model of DIS on a proton |182j . Together with a simple model for no, the re- 
sulting Ng from Eq. (|328p was used in Eq. (|232|) to give the gluon production cross section. At forward 
rapidities high-p^- particle production is close to the kinematic limit, which is obtained by demanding 
that projectile's Bjorken x is less than 1, Xp = {kT/ \/s) < 1, leading to kx < \fse^ . The effective X\ 
is rather large and valence quark contribution becomes important. The latter is given by Eqs. ()270p 
and 1)2711) and was also taken into account in |122j . To compare the quark and gluon production cross 
sections with the data one indeed needs to convolute them with fragmentation functions. The resulting 
fit of BRAHMS data from Fig. EBlis shown in Fig. 1571 (see also ^123 ). Since BRAHMS data at rapidities 
?7 = 0, 1 is for total charged hadrons, while the data at rapidities ?7 = 2.2, 3.2 is for negatively charged 
hadrons, the charge asymmetry issues had to be taken into account |183j . One may also worry that in 
the given kinematics the saturation-based description of high-p^ hadron production in collisions may 
not be well-justified: in that sense the success of the saturation model from ||122l| in describing particle 
spectra in collisions can be considered as just a parameterization of the data. What is important 
is that saturation correctly describes particle production in ci + Au collisions, i.e., the numerator in 
Eq. ()209|) . To model non-perturbative effects a shift of the saturation scale k was introduced in |122j : 
however, this shift can be put to zero without significantly impacting the quality of the fit, as one 
can see from Fig. EZI (The non-perturbative shift appears to improve the agreement with the data at 
mid-rapidity (the dashed line on the left panel in Fig. Wl\ indicating that non-perturbative effects may 




Figure 57: Fit of nuclear modification factor RdAu of charged particles for different rapidities from |122j . 
In the top two figures, corresponding to = 0, 1, the solid line corresponds to (hr + h^)/2 contribution 
calculated with k = in the isospin-independent approximation, while the dashed line gives the same 
{h~ + h^)/2 contribution but with k = 1 GeV. In the lower two plots, corresponding to t] = 2.2, 3.2, the 
solid line gives the h~ contribution calculated in the constituent quark model with k = 0, the dashed 
line gives the same h" contribution for n = 1 GeV, while the dotted line at ?7 = 2.2, 3.2 gives the 
(/i^ + h^)/2 contribution with k = (see text). Data is from (FT] . 
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still be important at corresponding values of Bjorken x.) We, therefore, conclude that BRAHMS data 
from Fig. |3ni[57j is in both qualitative and quantitative agreement with the saturation/CGC physics. 

The model from jl22j is based on Eq. (|23ip with the unintegrated gluon distribution functions 
calculated in the saturation/CGC approach. In that sense it is consistent with the model of jSHI used 
to describe total charged particle multiplicity in d + Au collisions, as shown in Fig. |^ Both models 
are inspired by saturation physics, but reflect different aspects of it: the model from |122ij deals with 
the high-pT hadronic spectra, and is, therefore, more accurate at high-pT- At the same time the model 
in jnS] concentrates on total particle multiplicity, which is most sensitive to low-pr hadron production: 
therefore, that model is more precise at Xow-px- 




Figure 58: Data from jS7j on nuclear modification factor Rcp of charged particles {h^ + h~)/2 for 
rapidities 77 = 0, 1 and of negatively charged particles h~ for r] = 2.2, 3.2 superimposed with the results 
of the model from |122j . Full and open dots, described by the solid and dashed lines correspondingly, 
give the ratio of particle yields in 0-20% and 30-50% centrality events correspondingly divided by the 
yields from 60-80% centrality events scaled by the mean number of binary collisions [HZ] • 

In Fig. EH we show BRAHMS data on R'-'^ for the same rapidities as in Fig.|H3 Full dots correspond 
to more central collisions, while open dots correspond to more peripheral collisions (see the caption of 
Fig. EH)- One can see that, in agreement with the expectation of saturation/CGC (see Sect. 13. the 
Cronin maximum at r/ = is an increasing function of centrality. At forward rapidities the centrality 
dependence changes: R^^ becomes a decreasing function of centrality, again in agreement with expec- 
tations of saturation/CGC approach (see Sect. 13. 2|) . On top of this qualitative agreement between the 
data and the theory of CGC, in Fig. EH "we observe quantitative correspondence of the data and the 
results of the saturation- inspired model from |122j . 

There have been several attempts to explain the suppression of particle production from Fig. ESI in 
the frameworks alternative or parallel to CGC. In |183|,ll84j the suppression of jHZj was analyzed in the 
framework of standard shadowing models. The conclusion of |183j was that the observed suppression is 
much stronger than one would anticipate using the common parameterizations of nuclear shadowing.^ 
In |185j it was argued that multiple rescattering in the nuclei would lead to observed suppression: how- 
ever, this conclusion appears to be in contradiction with our Sect. 13.11 where we prove that multiple 
rescatterings lead to enhancement of the nuclear modification factor. Recently a new model was pro- 
posed in |186j in which the origin of suppression in Fig. 1531 was attributed to energy conservation which 
was imposed by Sudakov from-factors. The model of 186j includes both multiple rescatterings and 
small-x evolution: hence at the moment the conceptual differences between the CGC approach and the 

successful fit of the data from '57' by shadowing models would only produce new parameterizations of quark and 
gluon shadowing, without calculating them from first principles, as we have done here in the CGC framework. 
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model from |186j are not clear. Parton recombination models also appear to be able to accommodate 
the data by completely redefining the hadronization mechanism ^187j . 

4.2.2 Future Experimental Tests 

Despite the successes of the Color Glass Condensate in predicting the suppression of the hadron pro- 
duction cross section in proton (deuteron)-nucleus collisions as compared to proton-proton collisions in 
the forward rapidity region at RHIC and the change of centrality dependence as one goes from mid 
to forward rapidity, one would like to quantitatively understand to which degree CGC is the domi- 
nant physics at forward rapidity RHIC and LHC collisions and possibly at mid-rapidity LHC collisions. 
So far, the main evidence for CGC at RHIC comes from the pt spectrum of single inclusive hadron 
production in deuteron-gold collisions in the forward rapidity region. 

A difficulty with a quantitative understanding of the physics behind the observed suppression of 
the hadron spectra in deuteron-nucleus collisions at RHIC is that hadronization by its nature is non- 
perturbative and is not well-understood from first principles. To avoid non-perturbative physics or to, at 
least, minimize its effects, one needs to consider high pj- processes. Exactly at what pt non-perturbative 
physics can be safely neglected can not be addressed by pQCD and one needs experimental input. The 
RHIC data on transverse momentum spectra of pions produced in proton-proton collisions in forward 
rapidity suggest that NLO pQCD works well down to pr ~ 1.5 — 2 GeV |188j . Therefore one may expect 
the weak coupling methods employed by CGC to work at similar p-r's for differential cross sections, even 
though perturbative approaches usually work much better for ratios of cross sections due to possible 
(partial) cancellations of large corrections in both the numerator and the denominator. 

The experimental situation is complicated by the fact the forward rapidity data published by 
BRAHMS is for negatively charged hadrons, production of which is suppressed, as compared to neutral 
or charge averaged hadrons, in proton-proton collisions. This can bias the nuclear modification factor as 
pointed out in |183j . However, the preliminary data on neutral pion production from STAR at rapidity 
T) = 4 shown above in Fig. ESI demonstrates clear suppression up to ~ 2 GeV. If this preliminary 
STAR data is confirmed, it will help clarify the role of isospin effects in the observed suppression of the 
hadron spectra at forward rapidity. 

The situation will improve drastically at LHC where due to an order of magnitude higher center of 
mass energy, the saturation scale is expected to be higher so that the weak coupling methods of CGC 
will be much more reliable and the available phase space in transverse momenta will be much larger 
even at rapidities as high as 5 — 6 which should eliminate or minimize phase space edge effects even for 
moderately large transverse momenta. 

Nevertheless, one can further test applicability of the Color Glass Condensate at RHIC by considering 
its electromagnetic signatures, for example in photon and dilepton production considered earlier. As 
noted above, the CGC building blocks for single particle production cross sections in proton-nucleus 
collisions are the fundamental and adjoint dipole cross sections, describing the scattering of quarks 
and gluons from the target. Therefore, one can "measure" or (more precisely) constrain the dipole 
cross section in one process and then use it to make predictions for other processes. This is in direct 
analogy with particle production cross sections in pQCD where one has universal (process independent) 
ingredients in the form of distribution and fragmentation functions, which are measured in one process 
and then used to make predictions for other processes, which give rise to the predictive power of the 
theory. 

Here we consider models of the dipole cross section (we consider those which include the small x 
evolution governed by JIMWLK or BK equations only) which have been used to fit the data, either in 
DIS at HERA or pA at RHIC. We then use these models to make predictions for photon and dilepton 
production cross sections and their modification factors at RHIC and LHC. We show that measuring 
the photon production cross section at RHIC can impose significant constraints on the models of the 
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dipole cross section and further clarify the role of the Color Glass Condensate at RHIC. 

The first model of the dipole cross section we consider is that of lancu, Itakura and Munier |189j 
which was used to fit the proton structure functions at HERA. It has the correct behavior in the 
saturation region and in the geometric scaling region but does not have the correct double logarithmic 
limit built in. Furthermore, it is fit to the data on proton targets and not on nuclear targets. The IIM 
model of the dipole cross section is 



where 



and 



(326) 
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The constants a, b are determined by matching the solutions at vtQs = "2 while 7^ = 0.63 and k = 9.9 
are determined from LO BFKL. The form of the saturation scale is taken to be = {xq/x)^ GeV^ 
with xo, A, -R, A/o determined from fitting the HERA data on proton structure function F2 




Figure 59: Invariant mass dependence of the nuclear modification factor for dilepton production in 
proton(deuteron)-nucleus scattering at RHIC. 

This model of the dipole cross section was used in |121j to make predictions for the transverse 
momentum integrated nuclear modification factor RpA for dilepton production cross section in deuteron 
(proton)-nucleus colhsions at RHIC. (For the py-dependence of dilepton RpA see |lUUj .) Here we show 
the main result of |121j for the nuclear modification factor for central collisions and aX y = 2.2 in Fig. EHl 
which shows a clear suppression in analogy with the nuclear modification factor in hadron production 
at RHIC that can be understood as being due to the quantum evolution of the dipole-nucleus cross 
section, present in both hadron and dilepton production cross sections. In the case of a deuteron 
projectile, nuclear shadowing (anti-shadowing) of the deuteron wave function is included by using the 
HKM parameterization |19()j . 

We also show the nuclear modification factor for central collisions and at rapidity y = 5 for the 
integrated over transverse momentum dilepton nuclear modification factor at LHC as a function of the 
dilepton invariant mass M in Fig. 1^01 where again a similar suppression is seen. Here one can extend 
the calculations to much higher dilepton masses since the saturation scale is very large due to the large 
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Figure 60: Invariant mass dependence of the nuclear modification factor for dilepton production in 
proton-nucleus scattering at LHC 



forward rapidity and the order of magnitude increase in the center of mass energy at LHC as compared 
to RHIC. 

We note that the IIM parameterization of the dipole cross section is originally done for a proton 
target. Here we have assumed a simple scaling of the dipole profile function by A^^^ in order to use it for 
nuclear targets. There is another parameterization of the dipole cross section, due to Kharzeev, Tuchin 
and one of the authors |122j , which has been used to fit the hadron production data in deuteron-nucleus 
collisions at RHIC. This model has the advantage that it has the correct high pt limit built in, unlike 
the IIM parameterization. It is given by 



N{rT,Xg 



1 — exp 



4 l""'^ N, 



(328) 



where the anomalous dimension 'y{y,r^) is parameterized as 
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and c is a constant fitted to the data. This parameterization has the correct high-p^ behavior where 
7 — i> 1. In Fig. I^we show the profile function for the two dipole parameterizations as a function of the 
dimensionless variable rrpQs for a nuclear target for Xg = 1.6 x 10~^ (which corresponds to Pt = 1.5 
GeV at y = 3.8 at RHIC). Clearly, the two dipole profiles are quite different. 

We now use the two parameterizations of the dipole cross section to calculate the photon production 
cross section in deuteron-gold collisions at RHIC. To get the cross section, we need to convolute the 
partonic production cross section in ()280|) with the quark (and anti-quark) distribution function in a 
deuteron. We therefore use 



(l(jd.{A)A^'y{kT,yj)X 

(PhdPkdy.^ 



(2vr) 



/ 



(^dipoleiXg^^S (331) 



to calculate the nuclear modification factor for photon production in deuteron-gold collisions at RHIC 
dX y = 3.8 and for the most central collisions. This will be eventually measured at RHIC by the STAR 
detector. 
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Figure 61: Quark anti-quark dipole profile for a nuclear target. 



It is clear that the two parameterizations lead to quite different results shown in Fig.lHSl This may be 
understood to be due to the different dipole shapes in the two approaches, so that a future measurement 
of this ratio can help constrain the models of the dipole cross section. Each parameterization of the 
dipole cross section has its advantages and disadvantages. For example, the JIM parameterization is 
a fit to the proton structure function, which involves a convolution of the dipole cross section with 
the photon wave function squared. Therefore, the dipole cross section is weighed preferentially in the 
physical cross section and some dipole sizes play a bigger role than others. The KKT parameterization 
on the other hand is a fit to particle production data which is a less inclusive observable than a fully 
inclusive structure function and therefore is more constraining on the dipole cross section. However, it 
is partly complicated by hadronization and the need for a convolution with the hadron fragmentation 
functions. Electromagnetic signatures of the Color Glass Condensate therefore offer a more precise 
means by which the quantitative aspects of CGC may be probed. 
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Figure 62: Nuclear modification factor for photon production. 



Having a parameterization of the dipole-nucleus scattering amplitude should allow one to make 
predictions for particle production at LHC. In Fig. 1^21 we present a prediction from (122|| for the nuclear 
modification factor Rp^ of charged particles at mid-rapidity at LHC (dashed line). It appears that the 
amount of suppression to be seen at mid-rapidity pA collisions at LHC is roughly equal to the currently 
observed at forward rapidity d + A collisions at RHIC (solid line in Fig. l6H|) . It would be very interesting 
to verify this prediction. If suppression at mid-rapidity pA collisions at LHC is observed, then one would 
be able to appreciate how special RHIC energy range is: while the center of mass collisions energy at 
RHIC is high, it is still low enough to have Cronin enhancement at mid-rapidity in d + Au collisions. 
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Figure 63: Nuclear modification factor R^^ of charged particles {h^ + h~)/2 at LHC energies ^/s = 
5500 GeV at mid-rapidity r/ = as predicted in |I22j (dashed line) versus i?''^" of (/i+ + /i~)/2 for RHIC 
energies y/s = 200 GeV at r/ = 3.2 (sohd line). 



allowing to use d + Au as a control experiment for the suppression of particle production observed in 
Au + Au indicating creation of quark-gluon plasma. 

Another interesting experimental test of Color Glass Condensate physics is in measuring two-particle 
correlations. While the exact numerical analysis of Eq. ()291|1 leading to numerical predictions to be 
verified experimentally has not been done yet, one can use the fact that for fciT, ^ Qs all the 
higher-twist multiple rescattering effects are not important and the two-gluon production cross section 
(generalized to pA) reduces to the following expression for gluon multiplicity |176j 



d^ki dyi d'^k2 dy2 d^B 
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j d'^qid^q2(ppi(f^,Y - y) f{k2 - q^,q^,y2 - yi] 

—X —2, 



X M{ki-q,r,yi), (332) 
where f{k2 — 2/2 ^ is given by the BFKL equation (pT| with the initial conditions from Eq. (j^n|) 




Figure 64: Two-gluon production in pA collision in /cj^-factorization approximation, which includes 
linear BFKL evolution only. 
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and both (j)p and (pA have hnear evolution only. Eq. ()332|) is the expression for two-gluon production in 
/cT-factorization formalism [176 applied to pA collisions. It is illustrated in Fig. IMl 

Following Kharzeev, Levin, and McLerran |l()lj . we apply Eq. ()332|1 to the case when the rapidity 
interval between the two produced gluons is very large. For instance, in RHIC d + Au kinematics, 
we are interested in the case when one of the gluons is produced at mid-rapidity and the other one is 
produced at forward deuteron rapidity. This is the case of Mueller-Navelet jets, which were originally 
proposed in |177j by Mueller and Navelet as an observable allowing to study experimentally the BFKL 
dynamics. 
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Figure 65: Prediction from |inij of suppression of back-to-back correlations of high-p^ hadrons produced 
m d + Au collisions as compared to pp collisions. 

Since Eq. ()332|) also applies to gluon production in pp collisions, one can compare the strength of 
back-to-back correlations in pA and pp collisions with one gluon at mid-rapidity and the other one at 
forward rapidity. The authors of |l(Jlj argued that BFKL evolution between the two produced gluons 
would deplete back-to-back correlations of produced gluons due to onset of anomalous dimension ()163|) 
in pA scattering. The effect is similar to suppression in RP"^ discussed above in Sect. 13.21 The resulting 
prediction for suppression of back-to-back correlations from |lUlj is shown in Fig. IHHl where the two- 
particle correlator is plotted as a function of the azimuthal angle between the particles for d + Au and 
pp collisions (with a constant background subtracted). Once can see that the strength of back-to-back 
correlations should be reduced due to the anomalous dimension effects in d+Au collisions. Experimental 
verification of the qualitative behavior predicted in Fig.inSlby the authors of |lUlj would be an interesting 
and important test of CGC Preliminary results from STAR collaboration |178j appear to confirm the 
qualitative expectation of suppression of back-to-back correlations for particles produced with a large 
rapidity interval between them, but more experimental tests are needed to reach a conclusion. 

5 Conclusions 

In writing a review of a rapidly developing field, one is bound to omit recent developments, which 
happen in parallel to the writing of the article. The topic of pomeron loop corrections to JIMWLK 
and BK evolution equations has received a lot of attention in the recent literature |124t I125| 11261 1127| 
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I128j . Pomeron loops are corrections to the non-linear evolution equations which are not enhanced by 
powers of color charge density, or, equivalently, powers of A. A typical pomeron loop brings in an 
extra (parametric) factor of exp{(ap — 1)Y} ^Tl I24j . Such corrections become important when 
ctg exp{(ap — 1) Y} ~ 1. Therefore, JIMWLK and BK evolution equations are valid up to rapidities 

Y < Yioop = ^-In^. (333) 

Eq. (j333|) . along with Eq. 0, describes the applicability region of JIMWLK and BK evolution equations. 
For rapidities higher than V/oop pomeron loop corrections become important and have to be included. 
In the recent papers |124t I125| I126t I127t I128j corrections to JIMWLK evolution were discussed, which 
would take into account pomeron loop contributions. The consequences of pomeron loop corrections 
for observables in pA collisions are not clear at present. 

Another important class of processes which we have omitted here are the exclusive and semi-exclusive 
processes. In the quasi-classical approximation the diffractive DIS structure functions were calculated 
in |1291 1130j . In the framework of dipole model one can construct an evolution equation governing 
single diffractive dissociation amplitudes for DIS, as it was done in |131j . The results of both |129| I13U] 
and [TBlJ can easily be generalized to pA collisions. Diffractive particle production was calculated for 
both DIS and pA in the quasi-classical approximation in [132 . Unfortunately, due to limited detector 
capabilities, it appears to be very difficult, if not impossible, to measure these exclusive processes at 
RHIC. Since our review here is dedicated to RHIC physics we had to omit such processes. 

We conclude by pointing out once again that we have reviewed the state of the field at the moment 
of writing. Future developments are likely to bring in many new exciting results, possibly modifying 
some of the discussed conclusions and adding on to the material presented here. 
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